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Abstract 

Near  and  Far-Field  Acoustic  Scattering  through  and  from  Two  Dimensional 

Fluid-Fluid  Rough  Interfaces 


Recent  experimental  results  [F.  E.  Boyle  and  N.  P.  Chotiros,  J.  Acoust  Soc.  Am.  91, 
2615-2619, 1992;  N.  P.  Chotiros,  J.  Acoust  Soc.  Am.  97,  199-214,  1995]  reveal  acoustic 
penetration  from  water  into  sandy  sediments  at  grazing  angles  below  the  compressional 
critical  angle  in  relation  to  the  mean  surface.  These  authors  interpret  their  results  to 
indicate  the  excitation  of  a  Biot  slow  wave  in  the  sediment.  Another  explanation  is 
considered  here.  Modeling  the  ocean  as  a  homogenous  fluid  and  the  sediment  as  a  lossy 
homogenous  fluid,  computer  simulations  of  these  experiments  based  on  analytical 
derivations  in  this  work  show  that  roughness  of  the  water-sediment  interface  causes 
propagation  of  acoustical  energy  from  water  into  the  sediment  at  grazing  angles  below 
the  compressional  critical  grazing  angle;  these  simulations  indicate  that  the  experimental 
results  can  be  explained  in  terms  of  diffraction  of  an  ordinary  longitudinal  wave.  These 
simulations  use  an  analytical  expression  for  the  time-dependent  mean  square  incoherent 
field  scattered  through  (and  from)  a  rough  2-D  fluid-fluid  interface  that  is  derived  in  terms 
of  the  bistatic  scattering  cross  section  per  unit  area  per  unit  solid  area  (differential  cross 
section)  of  the  rough  interface.  First-order  perturbation  theory  is  used  to  derive  an 
expression  for  the  differential  cross  section.  The  coherent  field  is  calculated  using  the 
flat-surface  result  (zero-order  perturbation  theory)  and  compared  to  the  coherent 
component  of  the  second-order  perturbation  theory  result.  Effects  of  sound-speed 
gradients  on  the  field  scattered  from  the  rough  water-sediment  interface  are  also  shown 
using  the  first-order  perturbation  derivations. 
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CHAPTER  1 


Introduction 


1.1  Purpose 

A  current  area  of  considerable  interest  in  ocean  acoustics  is  the  penetration 
of  acoustical  energy  from  a  source  in  water  into  the  seafloor  (Williams, 
Satkowiak,  and  Bugler,  1989;  Chotiros,  1995).  Recent  experimental  results 
(Boyle  and  Chotiros,  1992;  Chotiros,  1995)  reveal  acoustic  penetration  from 
water  into  sandy  sediments  at  grazing  angles  below  the  cOmpressional  critical 
angle  in  relation  to  the  mean  surface.  In  addition,  energy  appears  to  propagate 
in  the  sediment  at  speeds  slower  than  the  compressional  wave  speed,  in  a 
direction  concurring  with  Snell's  law.  These  authors  Interpret  their  results  to 
indicate  the  excitation  of  a  Biot  slow  wave  in  the  sediment. 

A  primary  goal  of  this  work  is  to  demonstrate  that  their  experimental  results 
can  be  explained  in  terms  of  diffraction  of  an  ordinary  longitudinal  wave,  rather 
than  refraction  of  a  slow  wave.  The  effect  of  roughness  on  penetration  of  the 
seafloor  is  demonstrated  by  developing  a  rough  surface  scattering  model  that 
provides  a  solution  to  scattering  into  the  rough  seafloor.  In  order  to  provide  a 
tractable  model,  simplifying  assumptions  about  the  ocean  and  seafloor  are  made, 
and  an  approach  for  solving  the  corresponding  time-dependent  rough  surface 
scattering  problem  is  taken  that  leads  to  simple  analytical  results.  The  ocean 
acoustic  model  used  and  the  methods  of  solution  are  discussed  in  the  following 
sections.  Although  the  main  problem  is  acoustic  scattering  through  the  rough 
fluid-fluid  interface  for  an  incident  field  in  the  water,  analytical  results  derived  in 
this  work  also  include  scattering  from  the  rough  fluid-fluid  interface.  Furthermore, 
these  analytical  results  are  general,  and  can  be  applied  to  other  rough  surface 
scattering  problems. 
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1.2  Ocean  acoustic  model 

The  propagation  of  sound  in  the  ocean  is  very  complicated  (see  for  example, 
Urick,  1975;  Brekhovskikh  and  Lysanov,  1991;  Burdic,  1991).  Sound  speed  is 
not  constant,  but  varies  as  a  function  of  temperature,  salinity,  and  depth.  In 
addition,  acoustical  energy  propagating  in  the  ocean  is  attenuated  by  absorption 
and  by  scattering  due  to  inhomogeneities  in  the  water.  However,  for  the 
examples  considered  here,  the  propagation  distances  in  the  water  are  short 
enough  so  that  the  sound  speed  does  not  vary  appreciably.  Furthermore,  at  the 
frequencies  of  interest  (300  Hz  -  60  kHz)  attenuation  of  acoustical  energy  over 
these  short  distances  is  also  negligible.  For  this  reason,  the  water  is  modeled  as 
a  lossless  homogenous  fluid. 

The  ocean  seafloor  is  in  general  an  inhomogeneous  porous,  viscoelastic 
medium,  possibly  supporting  the  propagation  of  Biot  slow  waves  (Biot,  1956a, 
1956b:  Biot,  1962a,  1962b)  and  transverse,  or  shear,  waves  (e.g.,  Brekhovskikh 
and  Lysanov,  1991)  in  addition  to  longitudinal,  or  compressional,  waves.  The 
seafloor  has  a  roughness  spectrum  approximating  that  of  a  power  law  (Fox  and 
Hayes,  1985;  Goff  and  Jordan,  1989;  Briggs,  1989).  Acoustical  energy  incident 
on  the  seafloor  will  be  scattered  from  this  rough  interface  back  into  the  water, 
and  through  this  rough  interface  into  the  sediment. 

The  main  purpose  of  this  work  is  to  show  the  effect  of  roughness  in  the 
seafloor  on  scattering  into  the  sediment.  Although  volume  inhomogeneities  in 
the  sediment  can  be  a  significant  source  of  scattering  (Ivakin  and  Lysanov,  1981; 
Jackson,  Winebrenner,  and  Ishimaru,  1986;  Jackson  and  Briggs,  1992),  volume 
scattering  is  beyond  the  scope  of  this  work.  However,  the  effect  of  vertical 
stratification  of  the  sediment  below  the  rough  interface  on  the  scattered  field  is 
discussed.  Boyle  and  Chotiros  (1992)  and  Chotiros  (1995)  neglect  the  effects  of 
roughness  in  their  interpretation  of  sediment  penetration  experiments,  and 
suggest  that  penetration  of  acoustical  energy  into  a  sandy  sediment  at  shallow 


3 


grazing  angles  (below  the  compressional  critical  angle)  can  be  explained  by  Biot 
wave  propagation  —  provided  nonstandard  Biot  parameters  are  used  (Chotiros, 
1995).  We  show  that  their  experimental  results  can  be  reproduced  in  a  computer 
simulation  that  models  the  sandy  sediment  as  a  fluid  and  the  interface  between 
the  water  and  the  sediment  as  randomly  rough.  Ivakin  and  Jackson  (in 
preparation)  show  that  for  the  case  of  sandy  sediments,  the  component  of  energy 
propagating  as  shear  waves  is  small,  and  can  therefore  be  neglected. 
Attenuation  of  acoustical  energy  is  not  negligible  in  the  sediment  for  even  very 
short  distances,  and  the  sediment  is  therefore  modeled  as  a  lossy  fluid. 

Since  it  is  not  feasible  to  obtain  an  exact  analytical  description  of  even  a 
relatively  small  portion  of  the  seafloor,  a  probabilistic  analysis  is  most  practical. 
A  statistical  quantity  of  particular  interest  to  rough  surface  scattering  is  the 
bistatic  scattering  cross  section  per  unit  area  per  unit  solid  angle  (Ishimaru, 
1978b).  For  convenience,  this  quantity  will  be  referred  to  as  the  differential  cross 
section.  It  is  useful  to  represent  solutions  for  scattered  acoustical  energy  in 
terms  of  differential  cross  section.  However,  exact  analytical  expressions  for  the 
differential  cross  section  are  not  available. 

1.3  Rough  interface  scattering  soiution  techniques 

Determining  the  field  scattered  from  or  through  (penetrating)  a  rough 
interface  given  an  arbitrary  incident  field  is  a  problem  with  applications  in  several 
research  areas,  in  addition  to  ocean  acoustics  (e.g.,  remote  sensing,  solid  state 
physics,  astronomy,  optics).  In  the  area  of  ocean  acoustics,  rough  surface 
scattering  is  viewed  as  a  contributor  to  propagation  loss  for  sound  propagating  in 
shallow  water  (e.g.,  Kuperman  and  Ingenito,  1977)  and  a  source  of  interfering 
reverberation  in  active  sonar  systems.  In  the  acoustic  detection  of  subsurface 
scatterers,  the  rough  seafloor  can  be  considered  a  source  of  interfering 
reverberation  and  signal  loss  (Nagy  and  Rose,  1993),  and  in  medical  imaging 
applications,  rough  interfaces  in  intervening  tissue  layers  often  lead  to  distorted 
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images  (Berkhoff,  Thijssen,  and  van  den  Berg,  1996).  Although  the  rough 
interface  scattering  problem  has  been  of  interest  for  several  decades,  there  is  no 
exact  solution  method  for  the  statistical  problem. 

Solutions  of  the  rough  interface  scattering  problem  are  obtained  by  means  of 
either  numerical  or  approximate  analytical  methods.  Although  numerical  methods 
can  result  in  accurate  calculations  for  the  scattered  field,  statistical  results  require 
a  numerically  intensive  “Monte  Carlo”  approach.  In  addition,  numerical  methods 
are  very  computer  intensive  for  the  three  dimensional  problem  (Tsang,  Chan, 
and  Pak,  1993, 1994;  Pak  etal.,  1995). 

A  classical  analytical  technique  applied  to  acoustic  scattering  by 
Brekhovskikh  (1952)  and  by  Eckart  (1953)  is  the  Kirchhoff  approximation,  or 
tangent  plane  approximation.  Although  this  approximation  has  a  useful  region  of 
validity  for  Dirichlet  boundary  conditions  (Thorsos,  1988a,  1988b,  1990;  Thorsos 
and  Jackson,  1989,  1991;  Ishimaru  and  Chen,  1990),  it  is  not  as  suitable  for 
obtaining  statistics  of  the  scattered  field  for  more  general  seafloor  boundary 
conditions,  such  as  fluid-fluid,  where  the  reflection  coefficient  is  a  function  of 
angle.  However,  as  shown  by  Thorsos  (1996b),  the  Kirchhoff  approximation  can 
be  very  accurate  for  scattering  through  a  rough  fluid-fluid  interface  near  the 
Snell's  law  direction. 

More  recent  approximate  methods  include  the  phase  perturbation  method 
(Winebrenner,  1985;  Winebrenner  and  Ishimaru,  1985a,  1985b),  the  small  slope 
approximation  (Voronovich,  1985),  and  the  operator  expansion  method  (Milder, 
1991).  A  summary  of  all  of  these  approximate  methods  is  given  in  Kaczkowski 
(1994).  See  also  Thorsos  and  Broschat  (1995)  for  a  description  of  the  small 
slope  method.  Although  the  phase  perturbation  method  can  provide  an  accurate 
approximation  to  the  scattered  field  for  Dirichlet  and  Neumann  surfaces 
(Broschat  et  al.,  1987;  Broschat,  Thorsos,  and  Ishimaru,  1989),  it  does  not 
readily  provide  tractable  solutions  for  the  fluid-fluid  seafloor  scattering  problem. 
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The  operator  expansion  method  is  accurate  in  several  scattering  regimes  for 
Dirichlet  surfaces  (Kaczkowski,  1994;  Kaczkowski  and  Thorsos,  1994),  but  has 
not  been  developed  for  other  boundary  conditions  and  does  not  allovy/  for  formal 
averaging.  A  Monte  Carlo  approach  is  therefore  required  for  obtaining  field 
statistics.  The  small  slope  approximation  can  yield  statistical  solutions  to  the 
scattering  problem  with  a  fluid-solid  boundary  condition  (Yang  and  Broschat, 
1994),  using  a  separate  supporting  analytical  technique  in  a  region  of  known 
validity.  The  analytical  technique  commonly  used  with  the  small  slope 
approximation  is  perturbation  theory  (Rayleigh,  1945). 

Perturbation  theory,  first  used  by  Rice  (1951)  for  randomly  rough  surfaces,  is 
a  straightforward  technique  for  obtaining  solutions  to  the  scattering  problem  with 
complicated  boundary  conditions  (Kuo,  1964;  Dacol  and  Berman,  1988; 
Winebrenner,  Farmer,  and  Joughin,  1995;  Ivakin,  1994).  The  Rayleigh-Rice 
method  (Rice,  1951)  is  a  classical  perturbative  approach  that  uses  the  Rayleigh 
hypothesis  as  a  starting  point.  In  this  approach,  the  incident  and  scattered  fields 
are  each  expressed  in  terms  of  a  Weyl,  or  a  plane  wave,  representation.  The 
Rayleigh  hypothesis  assumes  that  the  scattered  field  can  be  represented 
everywhere  in  terms  of  a  superposition  of  plane  waves  traveling  away  from  the 
surface.  Although  this  assumption  is  generally  not  valid  close  to  the  surface, 
Jackson,  Winebrenner,  and  Ishimaru  (1988)  have  established  that  the  Rayleigh- 
Rice  method  and  another  more  rigorous  form  of  perturbation  theory  based  on  the 
extinction  theorem  (Nieto- Vesperinas  and  Garcia,  1981)  yield  identical  solutions 
for  the  special  case  of  a  Dirichlet  boundary  condition  through  at  least  fifth  order. 
They  conjecture  that  the  two  perturbation  methods  yield  identical  results  for  all 
orders  and  all  boundary  conditions,  and  that  the  Rayleigh  hypothesis  is  irrelevant 
to  perturbation  theory. 

The  conditions  for  validity  of  the  perturbation  solution  depend  on  the  root 
mean  square  (RMS)  height  h  of  the  rough  surface  as  well  as  the  wave  number 
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ki  of  the  incident  wave.  It  is  generally  accepted  that  for  values  of  k\h  much 
smaller  than  one,  a  first-order  perturbation  solution  will  usually,  but  not  always, 
represent  the  scattered  field  with  sufficient  accuracy  (Thorsos  and  Jackson, 
1989).  For  this  reason,  first-order  perturbation  theory  is  widely  used  for  low 
frequency  scattering  from  rough  surfaces.  Kuo  (1964)  has  presented  an 
expression  for  the  differential  cross  section  of  a  randomly  rough  interface 
separating  two  fluids  based  on  first-order  Rayleigh-Rice  perturbation  theory.  This 
result  has  been  frequently  used  to  model  the  seafloor  (Crowther,  1983,  Jackson, 
Winebrenner,  and  Ishimaru,  1986;  Stanic  ef  a/.,  1988,  1989;  Mourad  and 
Jackson,  1989,  1993;  Jackson  and  Briggs,  1992).  Dacol  and  Berman  (1988), 
Kuperman  and  Schmidt  (1986),  Liu,  Schmidt,  and  Kuperman  (1993),  and  Essen 
(1994)  have  applied  perturbation  theory  to  a  rough  shear-supporting  sea-floor. 

Comparisons  of  first-order  perturbation  approximations  with  exact  numerical 
results  for  the  2-D  problem  (1-D  surface)  have  been  made  by  Chen  and  Fung 
(1988)  and  Thorsos  and  Jackson  (1989)  for  surfaces  with  a  Gaussian  roughness 
spectrum.  These  comparisons  indicate  that  perturbation  theory  can  be  valid  for 
values  of  kxh  as  high  as  025.  Thorsos  (1990)  shows  by  comparisons  with  exact 
integral  equation  results  for  scattering  from  a  1-D  rough  fluid-fluid  interface,  that 
first-order  perturbation  theory  accurately  determines  the  differential  cross  section 
for  scattering  from  a  "Pierson-Moskowitz"  roughness  spectrum  with  k\h  =  0.45. 
Using  a  truncated  power-law  roughness  spectrum  and  sound-speed  profile 
appropriate  to  the  examples  in  this  dissertation,  Thorsos  (1996a)  shows  by 
comparisons  with  exact  integral  equation  results  for  scattering  through  a  1-D 
rough  fluid-fluid  interface,  that  first-order  perturbation  theory  accurately 
determines  the  differential  cross  section  over  a  wide  range  of  scattering  angles 
with  k\h  =  l.  For  scattering  from  a  1-D  rough  fluid-fluid  interface,  Thorsos  shows 
that  perturbation  theory  is  valid  for  k\h  =  035.  Since  a  goal  of  this  work  is  to 
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provide  a  tractable  analytical  model  for  seafloor  scattering  that  is  particularly 
accurate  for  scattering  through  a  rough  interface,  perturbation  theory  is  used. 

1.4  Overview 

Representing  the  ocean  as  a  homogenous  fluid,  and  the  sediment  as  a  lossy 
homogenous  fluid  separated  by  a  rough  interface,  an  analytical  model  based  on 
perturbation  theory  is  used  to  create  computer  simulations  of  the  penetration 
experiments  conducted  by  Boyle  and  Chotiros  (1992)  and  Chotiros  (1995)  — 
even  though  these  authors  suggest  that  their  results  can  be  explained  in  terms  of 
slow  wave  propagation  in  sediment  that  is  a  Biot  medium.  The  general  formalism 
for  scattering  of  a  pulse  through  a  rough  interface  as  well  as  from  a  rough 
interface  is  derived  in  chapter  2  of  this  dissertation.  First-order  fluid-fluid 
perturbation  theory  results  used  to  approximate  the  rough  interface  scattering 
problem  in  the  simulations  are  derived  in  chapter  3.  In  chapter  4,  the  results  of 
the  corresponding  simulations  of  the  experiments  are  given  and  compared  with 
the  experimental  results  of  Boyle  and  Chotiros  (1992)  and  Chotiros  (1995).  The 
simulations  use  the  first-order  perturbation  model  developed  in  chapters  2  and  3 
to  calculate  the  incoherent  scattered  intensity  at  a  receiver  in  the  sediment,  and 
use  the  flat  surface,  zero-order  field  to  calculate  the  coherent  intensity.  Using  the 
zero-order  field  to  represent  the  coherent  intensity  is  valid  for  small  amounts  of 
surface  roughness;  however,  the  surface  roughness  used  in  the  ocean 
experiment  simulations  is  kih  =  l.  The  second-order  component  of  the  coherent 
intensity  is  derived  in  chapter  5,  and  compared  with  the  zero-order  component  of 
the  coherent  intensity  for  the  sediment  roughness  parameters  used  in  the  ocean 
experiment  simulation. 

The  simulation  model  also  neglects  gradients  in  sound  speed  in  the 
sediment.  In  the  approach  taken  in  chapter  5,  the  medium  below  the  interface  is 
constrained  to  be  a  fluid  only  down  to  the  lowest  point  of  the  rough  interface. 
Below  this  point,  the  sediment  can  be  an  arbitrary  vertically  stratified  medium. 
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characterized  by  its  flat  surface  reflection  coefficient.  An  application  of  the 
derivations  in  this  chapter  shows  the  effect  of  an  upward  refracting  sound-speed 
gradient  on  the  field  scattered  back  into  the  water.  This  result  is  further 
generalized  in  chapter  6,  where  the  sediment  parameters  below  the  lowest  point 
on  the  rough  interface  are  taken  to  be  random.  Although  results  showing  the 
effect  of  gradients  on  scattering  through  a  rough  interface  are  not  presented  in 
chapter  6,  this  case  can  also  be  examined  using  the  analytical  approach  in  this 
chapter  and  chapter  5.  Overviews  of  chapters  2  through  6  are  presented  in 
sections  1.4.1  through  1.4.5. 

1.4,1  Time-dependent  mean  square  scattered  field 

The  total  propagation  time  from  source  to  receiver  can  be  experimentally 
determined  when  the  source  signal  is  a  pulse.  For  a  known  propagation  path 
and  known  propagation  speed,  the  distance  of  an  object  can  be  determined  from 
the  propagation  time  of  a  scattered  pulse.  Likewise,  if  the  propagation  speed  is 
unknown,  but  the  propagation  path  and  locations  of  the  source  and  receiver  are 
known,  the  propagation  speed  can  be  determined  from  the  pulse  propagation 
time.  In  Boyle  and  Chotiros  (1992),  the  propagation  speed  in  the  sediment  is 
inferred  from  the  arrival  of  signal  energy.  The  main  purpose  of  chapter  2  is  to 
derive  a  near-field  expression  for  the  time-dependent  mean  square  scattered 
field  (average  of  square  of  field  magnitude)  due  to  a  narrowband  point  source, 
including  dispersion  due  to  scattering  and  frequency  dependence  of  wave  speed. 
These  results  (Moe  and  Jackson  1996a,  1996b)  are  applied  to  the  simulation  of 
tank  experiments  (Boyle  and  Chotiros,  1992)  and  ocean  experiments  (Chotiros, 
1995):  the  results  of  these  simulations  are  discussed  in  chapter  4. 

The  time-dependent  mean  square  scattered  field  can  be  expressed  in  terms 
of  a  two-dimensional  integral  over  the  two-frequency  mutual  coherence  function, 
r(coi,(02)  (see  Ishimaru,  1978a,  or  Ishimaru  et  al.,  1994a).  Scattering  of  a  pulse 
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from  the  ocean  surface  is  addressed  in  Ziomek  (1982a,  1982b),  where  r(coi,(02) 
is  presented  in  terms  of  a  Kirchhoff  representation.  Ishimaru  et  al.  (1994a)  derive 
an  analytical  expression  for  r(t0i,Q)2)  for  scattering  from  rough  surfaces  based 
on  the  Kirchhoff  approximation,  and  Ishimaru  etal.  (1994b)  use  a  second-order 
Kirchhoff  approximation  with  shadowing  corrections  (Ishimaru  and  Chen,  1990, 
1991)  to  obtain  an  analytical  expression  for  r(coi,co2)  that  includes 
backscattering  enhancement  effects.  In  chapter  2,  we  obtain  a  general 
expression  for  r(®i,o)2)  appropriate  for  rough  interface  scattering  for  a 
narrowband  incident  plane  wave  that  does  not  depend  on  a  particular  scattering 
approximation. 

In  order  to  obtain  simple  analytic  expressions  for  the  time-dependent  mean 
square  scattered  field,  Ishimaru  etal.  (1994a,  1994b)  assume  the  scattering 
channel  is  a  WSSUS  (wide-sense  stationary  uncorrelated  scattering)  channel 
(Ishimaru,  1978a)  —  the  two-frequency  mutual  coherence  function,  r(o)i,o)2)  is 
only  dependent  on  the  frequency  difference,  coj  =  co2-a)i,  and  not  on 
(coi  +  02)/2-  Because  of  frequency  dependent  attenuation  in  the  sediment,  we 
do  not  assume  that  r(coi,o)2)  is  only  dependent  on  co^.  However,  in  chapter  2, 
we  assume  that  the  incoherent  T-matrix  correlation  function  varies  slowly  over 
the  source  frequency  range,  and  expand  phase  terms  in  a  power  series  in 
frequency.  Both  the  incoherent  T-matrix  correlation  function  approximation  and 
the  power  series  expansion  require  that  the  source  signal  be  narrowband.  This 
approach  is  used  to  derive  a  general  analytical  expression  for  the  time- 
dependent  mean  square  incoherent  field  scattered  from  and  through  a  2-D  fluid- 
fluid  rough  interface  due  to  a  narrowband  incident  plane  wave  which  is 
expressed  in  terms  of  the  second  moment  of  the  T-matrix.  At  distances 
sufficiently  far  from  the  rough  interface  to  neglect  the  evanescent  waves,  this 
expression  for  the  time-dependent  mean  square  incoherent  field  is  given  in  terms 
of  the  differential  cross  section  of  the  rough  interface.  This  derivation  leads  to 
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the  usual  sonar  equation  in  the  limit  as  the  narrowband  signal  approaches  the 
CW  (continuous  wave)  case.  For  the  special  case  of  a  CW  incident  plane  wave, 
this  result  becomes  essentially  exact  for  distances  greater  than  about  a 
wavelength  from  the  interface.  This  result  is  surprising  because  the  differential 
cross  section  is  defined  as  a  far-field  entity  (Ishimaru,  1978b:  Winebrenner  and 
Ishimaru,  1986).  The  expression  for  the  time-dependent  mean  square 
incoherent  field  due  to  an  incident  plane  wave  is  used  to  obtain  the  expression 
for  the  scattered  field  due  to  a  narrowband  point  source  heuristically.  A  simple 
analytical  expression  is  obtained  for  the  case  of  a  Gaussian  shaped  source 
pulse,  also  in  terms  of  the  differential  cross  section. 

First-order  perturbation  calculations  (first-order  perturbation  results  are 
derived  in  chapter  3)  for  the  case  of  a  baseband  Gaussian  shaped  source  pulse 
illustrate  narrowband  pulse  dispersion  effects  of  the  incoherent  field  due  to 
forward  scattering  into  a  lossy  sediment.  For  the  case  of  incidence  below  the 
critical  grazing  angle,  first-order  perturbation  computations  also  show  that  the 
incoherent  field  scattered  through  a  rough  surface  can  be  much  greater  than  the 
zero-order  field  transmitted  below  the  corresponding  flat  surface  depending  on 
loss  and  receiver  depth.  These  computations  for  the  mean  square  incoherent 
field  penetrating  the  rough  interface  are  compared  to  the  flat-surface  result,  for 
both  plane  wave  and  point  sources. 

1.4.2  First-order  perturbation  theory 

Derivations  for  the  first-order  perturbation  expressions  used  for  the  numerical 
computations  in  chapter  2  and  in  the  computer  simulations  (chapter  4)  are 
presented  in  chapter  3.  Rayleigh-Rice  perturbation  theory  is  derived  in  chapter  3 
for  the  simple  case  of  a  rough  interface  separating  two  homogenous  fluids.  The 
water  above  the  rough  interface  is  modeled  as  a  homogenous  lossless  fluid,  and 
the  sediment  below  the  interface  is  modeled  as  a  homogenous  lossy  fluid.  The 
results  of  chapter  3  are  equivalent  to  those  published  by  Kuo  (1964)  and  Mourad 


11 


and  Jackson  (1989),  except  that  the  field  scattered  into  the  second  medium  is 
also  considered  here.  Although  these  results  are  derived  for  the  two 
homogenous  fluid  case,  this  form  of  the  first-order  scattered  field  solution  is 
shown  to  be  more  general  in  chapter  5.  First-order  perturbation  theory  solutions 
for  this  case  are  also  used  in  the  simulations  of  sediment  penetration 
experiments  in  chapter  4.  As  in  Moe  and  Jackson  (1994c),  a  truncated  or  filtered 
power  law  is  used  for  the  roughness  spectrum.  This  filter  is  shown  to  be 
equivalent  to  subtracting  a  moving  average  (Papoulis,  1984)  from  a  2-D  surface 
with  power-law  roughness  spectrum. 

1.4.3  Computer  simulations  of  recent  experiments 

A  recent  tank  experiment  (Boyle  and  Chotiros,  1992)  and  a  recent  ocean 
experiment  (Chotiros,  1995)  reveal  acoustic  penetration  from  water  into  sandy 
sediments  at  grazing  angles  below  the  compressional  critical  angle  in  relation  to 
the  mean  surface.  In  their  tank  experiment  (Boyle  and  Chotiros,  1992),  a 
receiver  buried  in  sand  measures  the  intensity  of  a  signal  transmitted  in  water 
above  the  sand.  The  propagation  time  of  the  transmitted  pulse  is  the  same  as 
that  of  slow  refracted  wave.  Since  only  a  single  receiver  is  used  in  Boyle  and 
Chotiros'  tank  experiments,  no  information  about  direction  of  the  propagating 
energy  in  the  sediment  is  available.  In  an  attempt  to  determine  both  direction 
and  speed  of  the  energy  propagating  in  the  sediment,  a  buried  array  of  receivers 
was  used  in  an  ocean  experiment  (Chotiros,  1995)  that  measures  the  intensity  at 
an  array  of  receivers  in  the  sediment  due  to  an  acoustical  source  in  the  water. 
They  use  a  simple  algorithm  to  determine  the  speed  and  direction  of  the  intensity 
front  propagating  in  the  sediment.  The  results  of  their  experiments  seem  to 
suggest  propagation  of  energy  at  speeds  slower  than  the  compressional  speed 
of  the  sediment,  in  a  direction  concurring  with  Snell's  law. 

Assuming  a  small  amount  of  roughness  at  the  water-sediment  interface,  and 
using  the  3-D  scattering  model  developed  in  chapter  2  with  the  first-order 
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perturbation  theory  results  for  2-D  fluid-fluid  rough  interfaces  derived  in 
chapter  3,  we  reproduce  the  results  of  the  tank  experiments,  indicating  that  the 
acoustic  penetration  of  the  surface  may  be  due  to  scattering  (diffraction)  from  low 
levels  of  roughness  (Moe  et  a!.,  1995)  rather  than  slow-wave  refraction. 
Assuming  roughness  parameters  appropriate  to  the  ocean  experiment  reported 
by  Chotiros  (1995),  our  goal  is  to  reproduce  his  results  with  the  model  developed 
in  this  dissertation.  An  explanation  by  Thorsos  (1995)  of  how  diffraction  can 
mimic  a  slow  wave  is  given  in  chapter  4,  along  with  discussions  of  simulation 
results  (Moe  et  al.,  1995).  Some  of  the  results  presented  in  this  chapter  have 
been  presented  in  Moe  et  al.  (1995)  and  again  in  Jackson  et  al.  (1996). 

1.4.4  Perturbation  theory  results  for  a  fluid-fluid  rough  interface  including 
gradients 

Moe  and  Jackson  (1994a)  follow  the  approach  used  by  Winebrenner, 
Farmer,  and  Joughin  (1995)  in  their  treatment  of  electromagnetic  scattering  from 
rough  dielectric  surfaces,  and  constrain  the  sediment  below  the  rough  interface, 
down  to  the  lowest  point  on  the  interface,  to  be  a  fluid.  However,  the  medium 
below  the  lowest  point  on  the  interface  is  only  constrained  to  be  vertically 
stratified  (o.g.,  sound  speed  can  vary  with  depth,  but  not  horizontal  position)  and 
is  represented  by  its  reflection  coefficient.  The  derivation  in  chapter  5  takes  this 
approach,  allowing  the  sediment  below  the  lowest  point  on  the  surface  to  be 
viscoelastic,  or  porous,  supporting  shear  or  Biot  slow  waves  as  well  as  gradients 
in  physical  properties.  A  first-order  solution  for  scattering  from  a  rough  interface 
is  also  given  in  Ivakin  (1994).  Multiple  layered  rough  interfaces  in  layered 
sediments  are  also  considered  in  his  presentation. 

First-order  perturbation  theory  is  valid  for  obtaining  the  incoherent  field 
scattered  into  the  sediment  for  the  roughness  parameters  used  in  the  simulations 
(Thorsos,  1996a),  but  generally,  there  is  a  need  for  higher  order  results  to  test 
validity  and  to  correct  lower  orders.  Also,  sediments  tend  to  be  stratified,  so  the 
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assumption  of  spatial  homogeneity  of  the  second  medium  is  sometimes 
questionable.  Using  Rayleigh-Rice  perturbation  theory,  Nth-order  recursion 
equations  are  derived  in  chapter  5  that  provide  a  starting  point  for  evaluating  the 
Nth-order  scattered  field  in  terms  of  lower  orders.  The  validity  of  using  the  zero- 
order  result  to  represent  the  coherent  intensity  is  shown  in  chapter  5,  where  an 
expression  for  the  second-order  coherent  field  is  derived. 

1.4.5  Formal  average  with  random  sediment  parameters 

In  the  simulation  discussed  in  chapter  4,  the  sediment  is  taken  to  be  a 
homogenous  fluid.  The  perturbation  derivations  in  chapter  5  allow  for  the  effects 
of  sound-speed  gradients  on  the  field  scattered  through  and  from  the  rough 
water-sediment  interface.  Strong  gradients  and  discontinuities  in  density  and 
sound  speed  on  scales  of  several  cm  are  often  observed  in  shallow  water 
environments  (Briggs,  Richardson,  and  Young,  1985;  Richardson,  1986;  Stanic 
eta!.,  1988, 1989;  Jackson  and  Briggs,  1992;  Lyons,  Anderson,  and  Dwan,  1994; 
Jackson  et  al.,  1996),  and  gradients  on  scales  of  hundreds  of  meters  are 
observed  in  deep  ocean  environments  (Hamilton,  1979,  1980).  Although  the 
approach  taken  in  chapter  5  includes  the  effects  of  layering  and  gradients  in 
seafloor  parameters  below  the  rough  interface,  it  does  not  allow  for  the  inherent 
randomness  of  these  parameters,  which  will  sometimes  "wash  out"  rapid  angular 
variations  in  scattering  caused  by  constructive  and  destructive  interference. 
Perfect  vertical  stratification  leads  to  interference  effects  which  may  not  always 
appear.  By  allowing  stratification  parameters  to  be  random,  these  interference 
effects  can  be  averaged  to  reveal  more  robust  features  of  the  differential  cross 
section.  Note  that  this  approach  does  not  model  volume  scattering,  since  there 
is  no  horizontal  dependence  of  sediment  properties. 

In  chapter  6,  the  first-order  expression  for  the  differential  cross  section 
derived  in  chapters  2  and  5  is  generalized  to  account  for  randomness  in  the 
medium  below  the  lowest  point  on  the  rough  interface  by  expressing  the 
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differential  cross  section  in  terms  of  moments  of  the  stochastic  flat  surface 
reflection  coefficient.  Upward  refraction  or  reflection  due  to  large  sound-speed 
gradients  in  the  sediment  results  in  increased  scattering  strength  over  the  non¬ 
gradient  case.  This  scattering  gain  occurs  for  various  sediment  types  over  a 
wide  range  of  bistatic  angles,  but  has  a  relative  maximum  when  the  incident  and 
scattered  grazing  angles  are  equal  (backscattering  enhancement). 
Backscattering  enhancement  is  observed  and  predicted  in  scattering  from  rough 
interfaces  (e.g.,  Maradudin  and  Mendez,  1993;  Ishimaruef  a/.,  1994b),  but  it  is 
not  manifest  in  first-order  perturbation  theory  calculations  that  involve  only 
scattering  from  a  single  rough  interface  with  no  additional  specular  or  refracted 
path.  Hanson  and  Zavorotny  (1995)  show  that  backscattering  enhancement  is 
apparent  using  first-order  perturbation  theory  in  a  two-scale  model  for 
electromagnetic  scattering  from  a  rough  ocean  surface.  In  Hanson  and 
Zavorotny,  the  energy  path  consists  of  a  single  scattering  event  and  specular 
reflection:  in  the  sound-speed  gradient  problem  addressed  in  chapter  6,  the 
energy  path  consists  of  a  single  scattering  event  and  refraction. 


CHAPTER  2 


Near  field  scattering  through  and  from  a  2-D  fluid-fluid  rough 

interface 


2.1  Introduction  to  chapter  2 

In  this  chapter,  a  general  analytical  expression  for  the  time-dependent  mean 
square  scattered  field  (average  of  square  of  incoherent  field  magnitude)  due  to 
an  incident  narrowband  plane  wave  is  derived  and  expressed  in  terms  of  the 
second  order  T-matrix  statistics,  as  well  as  in  terms  of  the  bistatic  scattering 
cross  section  per  unit  area  per  unit  solid  angle  (differential  cross  section)  of  the 
rough  interface.  This  novel  approach  yields  an  analytical  expression  that  is  not 
dependent  on  a  particular  scattering  solution  technique.  In  addition,  the  effect  of 
loss  in  the  medium  below  the  interface  is  included  in  the  result.  Even  though  the 
differential  cross  section  is  typically  considered  a  far-field  entity  (Ishimaru, 
1978a),  the  form  of  this  result  in  terms  of  the  differential  cross  section  is  valid  for 
distances  on  the  order  of  a  wavelength  from  the  interface.  This  rigorously 
derived  result  is  used  to  derive  the  scattered  field  due  to  a  narrowband  point 
source  heuristically.  This  derivation  leads  to  the  usual  sonar  equation  in  the  limit 
as  the  narrowband  signal  approaches  the  CW  (continuous  wave)  case.  The 
results  of  this  chapter,  along  with  the  two-fluid  perturbation  theory  results  of 
chapter  3,  are  used  in  the  simulation  of  sediment  penetration  experiments  (Boyle 
and  Chotiros,  1992;  Chotiros,  1995),  discussed  in  chapter  4.  The  results  of  this 
chapter  are  also  presented  in  Moe  and  Jackson  (1996a,  1996b). 

Figure  2.1  shows  a  2-D  slice  of  the  3-D  ocean  model.  A  time-dependent 
incident  field  in  a  homogenous  fluid  (water)  above  a  rough  water-sediment 
interface  results  in  energy  scattered  from  the  interface  back  into  the  water  and 
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through  the  rough  interface  into  the  sediment.  The  sediment  is  represented  by  a 
homogenous  lossy  fluid.  The  effect  of  depth-dependent  sediment  parameters  on 
the  time-dependent  scattered  field  is  not  considered  in  this  formalism.  However, 
in  the  general  perturbation  theory  derivation  in  chapter  5,  as  well  as  in  the 
formally  averaged  sediment  parameters  derivation  in  chapter  6,  the  effect  of 
sound-speed  gradients  on  the  intensity  scattered  back  into  the  water  for  the  CW 
case  is  discussed.  General  background  for  the  scattering  problem  is  given  in  the 
next  section. 

The  solution  for  a  time-dependent  scattered  field  requires  obtaining  the 
solution  for  the  two-frequency  mutual  coherence  function,  r(o)i,co2)  (Ishimaru, 
1978a).  A  general  expression  for  r(«)i,f02)  for  an  incident  plane  wave  is  derived 
by  assuming  the  T-matrix  correlation  function  is  a  slowly  varying  function  of 
frequency.  Dispersion  effects  due  to  frequency-dependent  phase  speed  are 
treated  by  expanding  phase  terms  in  a  power  series  about  the  center  frequency. 
The  validity  of  this  expansion  and  the  frequency-independence  assumption  both 
depend  on  the  source  signal  bandwidth.  The  validity  of  the  frequency- 
independence  assumption  for  the  T-matrix  correlation  function  is  shown  in 
chapter  3  to  also  depend  on  the  distance  from  the  interface.  Using  this 
approximate  expression  for  r(a)i,a)2),  we  obtain  a  near-field  expression  for  the 
time-dependent  mean  square  incoherent  field  (average  of  incoherent  field 
magnitude  squared)  for  narrowband  source  signals  in  terms  of  the  second 
moment  of  the  T-matrix.  Neglecting  the  evanescent  field  components  in  this 
result,  we  obtain  a  useful  expression  for  the  time-dependent  mean  square 
incoherent  field  in  terms  of  the  differential  cross  section.  This  result  is 
remarkable  in  that  a  near  field  result  is  obtained  in  terms  of  what  is  generally 
considered  to  be  a  far  field  quantity.  First-order  perturbation  results  (derived  in 
the  following  chapter)  depict  the  contribution  of  the  evanescent  component  to  the 
incoherent  scattered  field  for  a  CW  plane  wave  source.  The  resulting  profiles  of 
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mean  square  incoherent  field  versus  depth  are  compared  with  the  flat-surface 
case,  also  illustrating  the  effect  of  roughness  on  the  field  penetrating  a  fluid-fluid 
interface. 

Using  the  narrowband  time-dependent  mean  square  incoherent  field  result 
for  an  incident  plane  wave,  we  heuristically  obtain  an  analytical  expression  for 
the  time-dependent  mean  square  incoherent  field  due  to  a  point  source.  First- 
order  perturbation  results  illustrate  dispersion  of  a  pulse  with  Gaussian  envelope 
due  to  forward  scattering  into  a  lossy  sediment.  These  computations  are  carried 
out  using  a  low-frequency  cutoff  for  the  bottom  relief  spectrum,  and  are 
compared  to  the  flat-surface  case. 

2.2  Scattering  problem  background 

The  general  scattering  geometry  is  shown  in  Fig.  2.1.  An  arbitrary  pressure 
field,  \|r/(r),  is  incident  on  a  2-D  rough  interface  separating  the  homogenous 
lossless  fluid  in  medium  1  (water)  from  the  lossy  fluid  in  medium  2  (sediment). 
The  zero-mean  rough  interface  is  defined  by 

z  =  hf{R),  (2.1) 

where  R  is  the  transverse  component  of  the  three-dimensional  position  vector  r, 
and  h  is  the  root  mean  square  (RMS)  height  of  the  rough  interface.  The  total 
field,  xi/i(r),  above  the  rough  interface  is  a  sum  of  an  arbitrary  incident  field, 
\|r,(r),  and  the  resulting  scattered  pressure  field,  \^f  (r); 

¥i(r)  =  Vj(r) +  ¥/(*•)•  (2-2) 

The  boundary  conditions  on  the  interface  consist  of  the  continuity  of  pressure, 

Vl(r)|2=y(R)  =  V2(*’)l2=;,/(R)> 


and  the  continuity  of  normal  velocity. 
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•  V\l/i (r)|^^ "  p2  ” 

where  \|/2  is  the  pressure  field  in  the  sediment,  h  is  the  unit  vector  normal  to  the 
interface,  pi  is  the  density  of  the  water,  and  p2  is  the  density  of  the  sediment  The 
scattered  field,  v|r/(r),  in  the  region  above  the  highest  point  on  the  interface  and 

the  incident  field  are  expressed  in  terms  of  the  following  (Weyl)  plane  wave 
expansions  (Devaney  and  Sherman,  1973): 

(r)  =  J  (2.5) 

and 

x|/^.(r)  =  J  ,  (2.6) 

where  K  is  the  two-dimensional  transverse  wave  vector  with  magnitude  K,ki  is 
the  wave  number  in  the  medium  above  the  interface,  and 

Pi(K)  =  Vr^^^.  (2.7) 

The  square  root  in  Eq.  (2.7)  is  chosen  so  that  p]  is  either  positive  or  positive 

imaginary.  For  the  CW  (continuous  wave)  case,  a  factor  of  is  suppressed, 
but  in  general  the  above  equations  are  the  Fourier  transforms  of  the  time  domain 
field, 

¥!(**)  =  Vi(*’,(o)  =  Jffr\iri(r,r)e'“'.  (2.8) 

The  field  below  the  rough  interface  is  expressed  as  a  plane  wave  expansion  of 
waves  traveling  in  the  negative  z  direction, 

Vi2(r)  =  V2-(r)^JA>P2-(K)e-'^l*2<*'>V*’',  z<0,  (2.9) 

where 

^2(K)  =  ^l-K^/ki 


(2.10) 


is  the  complex  wave  number  in  the  lossy  medium  below  the  interface  with 
frequency-independent  loss  parameter,  5.  The  speed  ratio,  v,  is  the  ratio  of  the 
sound  speed  in  the  sediment,  C2,  to  the  sound  speed  in  the  water,  q.  Although 
the  definition  of  ^2  'n  Eq-  (2.11)  violates  causality  (Aki  and  Richards,  1980; 
Wingham,  1985),  we  show  these  noncausal  effects  are  negligible  and  consistent 
with  the  approximations  made  in  deriving  the  analytical  result.  The  square  root  in 
Eq.  (2.10)  is  chosen  such  that  P2  'S  in  the  first  quadrant  of  the  complex  plane. 
The  field  scattered  back  into  medium  1  can  be  expressed  in  terms  of  the  incident 


field  and  the  T-matrix,  T\i: 

^f{K)  =  J  A'7ii(K,K')'Fj(K').  (2.12) 

Likewise,  the  downward  component  of  the  scattered  field  penetrating  below  the 
interface  is  expressed  in  terms  of  the  T-matrix,  ri2: 

'i'2-(K)  =  I  ATi2(K,K')'Pj(K').  (2.13) 

The  scattered  field  for  the  case  of  an  incident  plane  wave  is  found  by  substituting 

%{K)  =  b{K-Ki)  (2.14) 

into  Eq.  (2.12)  and  the  result  into  Eq.  (2.5),  yielding 

\l/y(r)  =  (2.15) 


The  scattered  field  penetrating  the  interface  for  the  incident  plane  wave  case  is 
found  in  the  same  way  from  Eqs.  (2.9),  (2.13),  and  (2.14): 

xj/2_(r)  =  I  A7i2(K,K,)e"'^2p2(K)Zg/KR 


(2.16) 
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For  convenience,  let  V(/(r)  represent  both  \\ff{r)  and  \K2-(r)-  Define 

\jjf(r)  =  \|/(r)  -  (\ir(r)) ,  (2-17) 

where  (\|/(r)^  is  the  mean  of  the  field  averaged  over  the  ensemble  of  rough 
interfaces.  Similarly,  define 

f{K,Ki)  =  T{K,Ki)-{T{K,Ki)) 

to  be  the  incoherent  T-matrix.  Using  the  above  definitions,  \j/(r)  is  expressed  as 
\j/(r)  =  J  (K,K,)e"'*2p2(K)Ze**^ .  (2.19) 


2.3  General  result  for  incident  plane  wave 

The  purpose  of  this  section  is  to  derive  an  expression  for  the  time-dependent 
mean  square  incoherent  field  at  the  receiver  position,  r,  due  to  a  narrowband 
plane  wave  source.  Although  the  following  derivation  considers  the  more  general 
case  of  scattering  through  a  rough  interface,  it  also  applies  to  scattering  from  a 
rough  interface. 

2.3.1  Time  domain  received  signal 

A  narrowband  source  signal,  s{t),  with  units  of  pressure  and  with  Fourier 
transform 

5(co)=  (2-20) 

— oo 

can  be  represented  as 
5(r)  =  Re|M(r)e“'“‘'*|, 


(2.21) 
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where  u{t)  is  the  corresponding  complex  baseband  input  signal  with  Fourier 
transform  t/(o)),  and  is  the  carrier  frequency,  or  center  frequency  of  the  pulse 
(see,  for  example,  Ishimaru,  1978a:  Proakis,  1989).  An  incident  plane  wave 
pulse  with  direction  denoted  by  the  transverse  unit  vector,  a,-,  and  incident 
grazing  angle,  6,  ,  is  expressed  in  terms  of  its  2-D  Fourier  transform  as 

'F,(K,(o)  =  5(a))5(K-K,)e'®^^®‘"®'^,  a)>0,  (2.22) 

where 


K,=  — cose,6c/  (2.23) 

is  the  incident  transverse  wave  vector,  and  is  taken  to  be  an  arbitrary  source 
height.  Using  Eqs.  (2.9),  (2.13),  (2.17),  (2.18),  (2.20),  (2.21),  and  (2.22)  and 
chapter  5  of  Ishimaru  (1978a),  we  express  the  incoherent  field  %fe(r,r)  at  a 
position  r  =  (R,z)  below  the  interface  as 

^bb  M  =  .  (2.24) 

where 


+  (dc)  =  \d\[  fi2(K,K,.,Q)'  +  0)  Je' 


xe 


-iKP2  (K,(0'+C0c  )2(a)'+(Dc  )/q  JKR 


(2.25) 


is  the  Fourier  transfer  of  the  impulse  response  of  the  incoherent  field,  generally 
known  as  the  transfer  function.  The  complex  wave  number  ratio  is  given  by 


K  =  ^2/^1 


(1  +  /5) 

V 


(2.26) 


and 


We  have  explicitly  shown  the  dependence  of  P2(K)  on  frequency  and  K. 


2.3.2  Time-dependent  mean  square  incoherent  field 

The  time-dependent  mean  square  incoherent  field  at  position  r  can  be 
expressed  as  (Ishimaru  1978a,  chapter  4) 

/2(r,O  =  (|%fc(r>0f)- 

Combining  Eqs.  (2.24)  and  (2.28)  yields  Ishimaru’s  (1978a)  Eq.  (5-17): 

/2(r,0  =  [^\ d(ii"[u{(d')uXGi") 

xr(G)'  +  C0c.tt>"  +  ®c)}>  (2.29) 

where 

r(co'  -I-  (Oc ,  CO"  +  COc  )  =  (^12  (®'  +  )‘^12  +  ®c  )) 

is  defined  by  Ishimaru  (1978a)  as  the  two-frequency  mutual  coherence  function. 
Using  the  expression  for  ??^i2  in  Eq.  (2.25),  Eq.  (2.30)  is  equivalent  to 

r(a)' -I- 0)^,0)"  +  (Oc) 

=  jd^K"  J  A'  (f,2(K',K;,(0'  +  co,)f,2(K",K",0)'’  +  «,)) 

^  g/K'Rg-iK"-Rg/  sin  0,z,  co'/q  ^-i  sin  0,z,  (a"/ci 
^g-/Kp2  (K',(0'+C0c  )z(co'+co,  )/c,  ^/(Kp2  (K",{o"+(Dc  ))*  z(a)"+(o,  )/c, 

where  the  transverse  wave  vectors  for  the  incident  field  at  the  angular 
frequencies  co"  and  co'  are 


(2.31) 


23 


cos9/a/  and  =  cosG/a,-. 


(2.32) 


Define  ^]2(K',K",K-,K",0Dc +  a)',(o^ +  ©")  (see,  for  example,  Zipfel  and  DeSanto, 
1972;  Voronovich,  1995)  such  that 

(f,2  (K',  Ki.o)' + <0,  )f;2  (k",k;'  CO" + ffl, )} 

=  ^i2(K;'’K",K-,K-',C0c  +  co'.cOc  +  a)")6(K"  -  K'  +  K-  -  K-O.  (2.33) 

We  refer  to  as  the  incoherent  r-matrix  correlation  function.  In  the  above 
expression,  the  subscript  is  used  to  represent  scattering  into  the  sediment  and 
the  dependence  of  fx2  on  cd  is  included  in  the  argument.  Substituting  Eq.  (2.33) 
and 


into  Eq.  (2.31),  and  changing  the  integration  variable  K'  to 
K  =  K'  +  K^/2, 
yields 

r(co'-i-roc’®''  +  ®c) 

=  J  {^i2(K  -  Kj/2,K  +  Kj/2,K;,K;;cOe  +  co',®^  +  ®") 

^  g-/K^  Rgi  sin  e,Zj  co'/cj  sin  e,z^  ©"/ci 

^^-iKp2(K-K^/2,©<,+®')z©7c, 

^^/(K32(K-K^/2,®,+®")r^K+®")Ai| 


(2.34) 


(2.35) 


(2.36) 


Equation  (2.36)  is  now  evaluated  by  making  a  few  approximations.  We  assume 
^12  is  a  slowly  varying  function  of  frequency,  and  for  the  case  of  a  narrowband 
source  signal, 
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^12(K-Kj/2,K  +  K^/2,k;,k;',co^  +  o)",0)c+co') 

=  t^i2(K,K,K,-,K;,a)c,(0c) 

^^12(K,K,-,coJ,  (2.37) 

where  the  relation  (Thorsos  and  Jackson,  1989;  Berman,  1992) 

^,2(K,K,.,»,)8(K-K'>(7i2(K,K,.,(0,)7i*2(K",K, ■,«,))  (2.38) 

is  a  special  case  of  Eq.  (2.33),  and  K,  is  evaluated  at  the  center  frequency: 

K.=^os0A-.  (2.39) 

q 

It  is  easy  to  show  using  first-order  perturbation  theory  (see  chapter  3)  that 
the  approximation  in  Eq.  (2.37)  is  valid  for  isr  +  |K^|/2  <  ,  or 

K<ki-  Tcfd  cos  0  j-  /ci ,  (2  -40) 

where  Inf^  =  |(o"-(o'|  can  be  taken  to  be  the  source  signal  bandwidth.  For  larger 
K,  the  phase  of  the  quantity  on  the  left  side  of  Eq.  (2.37)  cannot  be  considered 
independent  of  frequency.  The  first-order  T-matrix  representing  scattering  from 
the  water  into  the  sediment,  ri2,  is  derived  in  chapter  3  and  given  in  Eqs.  (3.41) 
and  (3.45),  and  the  first-order  approximation  of  the  incoherent  T-matrix 
correlation  function  is  given  in  Eq.  (3.50). 

Define  the  vertical  component  of  the  wave  vector  in  medium  2  as 

co^  +  co',  cOc  +  ® ")  =  ((cOc  +  0)')/  Cj  )Kp2  (K  -  /2 ,  ©c  +  to') .  (2.41) 

Analytical  expressions  for  dispersion  in  a  waveguide  can  be  found  by  expanding 
the  propagation  constant  of  the  time-dependent  field  in  a  power  series  (for 
example,  see  Ishimaru,  1991,  chap.  4).  Similarly,  with  reference  to  Eq.  (2.36), 
((o'/ci)kP2  is  expanded  in  a  power  series  in  co'  and  to".  For  the  case  of  a 
narrowband  signal. 
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B2{K,(iic  +  CO',CO^  +  CO")  =  B2c{K)  +  ^(ii'  +  ^(0" 

^  a 


9  ^  ,2  ^  ^  Bi  ,  „  9^52  »2 

3(0'^  9  0)9(0  3(0"^ 


(2.42) 


where  the  derivatives  are  evaluated  at  co'  =  o)"  =  0,  and 


B2,(/:)^K^ip2(K,®c).  (2-43) 

where  the  wave  number,  ki ,  in  the  above  equation  and  all  following  equations  is 
evaluated  at  the  center  frequency  (ifei  =o)c/ci).  Solving  for  the  derivatives  in  Eq. 
(2.42)  (see  Appendix  A),  combining  like  powers  of  co'  and  co"  in  Eq.  (2.36),  and 
substituting  the  result  into  Eq.  (2.29)  yields 


;2(r,r)  =  J/Ji:  |^i2(K,Ki,ra,)e"‘' 

x±ld<0"  £/*(0)")e''“"('-'‘  )■  e'”"' 


where 


R  n  Zs  .  a  d/  Kcos0,-  I 

h  = - ^cos0j  +  -^sm0j  +  Z— - ;; - ^Re  „  - r 

*  q  q  '  q/:i  l^Kp2(K,(0c) 


(2.44) 


(2.45) 


and 


t  =_z£]S/£l_ 
2  P2(K,®c)’ 


(2.46) 


The  physical  significance  of  these  parameters  will  be  discussed  later.  Note  that 
?!  is  real  and  ti  is  complex.  The  complex  coefficient. 


P  =  |(2k^/ )(«!•  •  K  (Oc  cos  0j  /q^  -  {Kjc^  - 1  Im(l/52e  )(cos  0,-  /q  f  _ 

-iTm(l/ 5|c)(di  •  Kcos0,/qf  |z 


(2.47) 
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and  the  real  coefficient, 

^^  =  Im|(d,'  •Kcose;/B2c)(K^®cA?  -«/  •*^cosejy(2c?))-cos^ei/(2cfB2c)|z,  (2.48) 

can  represent  the  contribution  to  pulse  dispersion  of  the  frequency-dependent 
phase  term.  This  result  is  simplified  by  setting 

2'7t) 

(2-49) 

2k  ^ 

The  above  definition  leads  to  a  simple  expression  for  the  time-dependent  mean 
square  incoherent  field  in  medium  2: 

In  addition  to  being  a  function  of  the  complex  argument,  r,  g(t)  is  a  function  of 
K,  the  incident  field  direction,  and  fluid-sediment  parameters.  Note  that 

+  (2.51) 

is  also  a  complex  function  of  K.  For  an  arbitrarily  small  pulse  bandv\/idth  (CW), 
F,  P,  Im(r2)  are  arbitrarily  small,  and  g(t)-^u(t). 

The  definition  of  k  in  Eq.  (2.25)  violates  causality  (see  Aki  and  Richards, 
1980,  chap  5,  and  Wingham,  1985).  When  a  corrected  version  of  Eq.  (8)  from 
Wingham  (1985)  is  used,  causality  will  not  be  violated  if 
K-26ln|(co^  +  co')/a)c|/(7rv)  is  substituted  for  K  in  (A2).  Using  this  substitution, 

and  following  the  derivation  in  Appendix  A,  results  in  causal  versions  of  the 
parameters  P,P,t\,  and  t2-  Since  6«1  for  ocean  sediment,  the  difference 
between  the  causal  representation  of?,^^/i,  and  t2  and  Eqs.  (2.45)  -  (2.48)  is 
negligible. 
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While  Eq.  (2.50)  involves  an  integral  over  the  transverse  wave  vector,  it  can 
be  converted  to  an  integral  over  the  scattering  interface  through  the  change  of 
variable, 

R'  =  R-vrjK/*l.  (2.52) 

As  shown  in  Fig.  2.2,  R'  is  the  location  of  a  small  scattering  patch  on  the 
interface,  and  rj  =  |r^|  is  the  distance  from  this  patch  to  the  field  point,  r . 


(2.53) 


We  constrain  K  to  the  range  in  Eq.  (2.52),  that  is,  we  neglect 

evanescent  waves.  This  constraint  is  justified  later  when  the  contribution  of 
evanescent  waves  on  the  mean  square  scattered  incoherent  field  is  illustrated 
using  first-order  perturbation  computations  for  an  incident  plane  wave.  This 
constraint  is  also  consistent  with  the  frequency-independence  assumption  of  the 
T-matrix  correlation  function.  This  assumption  is  invalid  for  arbitrarily  large 
integration  variable  K  in  Eq.  (2.50)  (see  Eq.  2.40).  The  Jacobian  follows  from 
Eq.  (2.52). 


|_ 

yd+z^  xjyd 

[dR'J 

Xdyd 

(2.54) 


where  and  are  the  transverse  coordinates  of  r^.  The  geometric 
significance  of  these  definitions  is  apparent  when  one  assumes  that  the  loss 
parameter  is  small  (5  « 1) .  Then 


«:ikP2(K.®c)  =  a/K?^ 

=  (itj  sin  02/v)(l  +  iS/  sin^  62). 


(2.55) 


where 


sin  02  =  v'^/ki  =  \z\/r^ 


(2.56) 
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is  the  sine  of  the  scattered  grazing  angle  defined  in  Fig.  2.2.  When  terms  of 
second  order  and  higher  in  the  loss  parameter  are  dropped,  the  significance  of  ti 
and  t2  becomes  apparent,  as 


?!  =  R'  •  OLi  COS  9,-  /q  +  Zs  sin  0,-  /  q 

(2.57) 

is  the  time  required  for  the  incident  plane  wave  front 

with  direction 

(a/cosGjssine,)  to  travel  from  the  source  reference  point, 
scattering  patch,  (R',0),  and 

(0,0,zj,  to  the 

(2.58) 

is  the  time  required  for  a  spherical  wave  scattered  from  the  patch  at  (R',0)  to 
travel  through  the  sediment  to  the  field  point  (R,z),  z<0.  According  to  Eq. 

(2.55),  the  absorption  exponent  in  medium  2  is  given  by 

=  -^kl/sin  02  =  ^^{hYd  ■ 

(2.59) 

Substituting  Eqs.  (2.54)  and  (2.59)  into  Eq.  (2.50)  results  in 

/2(r,()=  frfV  /(r-tj). 

•'  V  ra 

(2.60) 

Substituting  the  following  quantities. 

a,2(a.a,.)  ^  *?™^fe)^|2(K.K,.(i),), 

(2.61) 

a  =  K/A:,  and  d, 

(2.62) 

into  Eq.  (2.60)  yields 

•’  r-rf 

(2.63) 
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The  interpretation  of  012  as  a  differential  cross  section  will  be  justified  later  in  this 
section.  Since  K  =  *i(R-R')/(vr^).  t\,  tj,  g,  and  therefore  F  and  F  are 
functions  of  |R-R'|.  Note  that  the  CW  result  follows  from  Eq.  (2.63)  by  simply 
setting  g(r  -  =  1  •  For  this  case  of  an  arbitrarily  narrowband  source  signal,  the 
result  in  Eq.  (2.50)  for  the  mean  square  incoherent  field  is  exact,  and  Eq.  (2.63) 
approaches  the  exact  result  at  distances  at  least  a  few  wavelengths  from  the 
interface,  where  the  evanescent  waves  contribute  very  little  to  the  mean  square 
incoherent  field. 

Although  the  mean  square  incoherent  field  scattered  into  the  water,  /i(r,r),  is 
similarly  derived,  the  corresponding  expression  for  is  most  easily  obtained 
as  a  special  case  of  Eq.  (2.63)  by  substituting  medium  1  parameters  for  medium 
2  parameters.  For  example,  kj  is  substituted  for  ^2.  Pi  for  KP2,  etc.  For  this 
special  case,  Eq.  (2.63)  simplifies  to 

(2.64) 

|r-r1 


where  from  Eq.  (2.61) 

o„(d,aj)  =  *?sm2(e2)|<?,i(K,K,.,m,)f.  (2.65) 

In  this  case,  62  is  also  the  scattered  field  grazing  angle  measured  from  the  mean 
horizontal  plane,  |z|/r^  =sin(e2).  Again,  the  quantity  is  found  from  the  relation 
(Eq.  2.38) 

^l,(K,K,..0),)5(K-K'')  =  (f„(K,K„ffl,)7i*,(K",K,.,(0,)).  (2.66) 

There  is  no  loss  in  this  case:  F=  0,  and  therefore. 


gfiiO- 


—  (d(o'U((o')( 
2n^ 


Plll4 


(2.67) 
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Pj  1  =  (2/pfie)(co^^:cose,'  cos  V cf  -  (/sT/qf ), 
and 


(2.68) 


(2.69) 


Note  that  tj  is  still  given  by  Eq.  (2.51),  but  the  propagation  time  from  the  surface 
patch  to  the  receiver/observation  point  above  the  interface  is  a  real  quantity. 


t 

2  pi(^,co,)‘ 


(2.70) 


As  given  in  Voronovich  (1994),  the  quantity  an(d,d,)  in  Eq.  (2.65)  is  actually  the 
differential  cross  section  defined  in  Ishimaru  (1978b). 


aii(d,di)  = 


|r-rf  (1^/0 

M  |,^,f  ' 


(2.71) 


where  \|r,  is  a  plane-wave  field  incident  on  a  surface  patch  of  area  M,  and  d,  is 
the  unit  vector  in  the  direction  of  propagation  of  the  incident  field;  vjiy  is  the 


scattered  coherent  field  at  a  long  distance  |r-r'p  from  the  interface,  in  the 
direction  denoted  by  the  unit  vector,  d.  Similarly,  012  ^^e  differential  cross 

section  relating  the  mean  square  incident  field  to  the  mean  square  scattered 
incoherent  field  in  medium  2  (the  sediment),  defined  as 


ai2(d,d,)  = 


[r-rf  (fr2-P)  2lm(^,)|r-r1 
M  |y.p 


(2.72) 


where  ^2-  's  the  incoherent  field  at  the  receiver  in  the  sediment,  a  distance 
|r-r'|  from  the  surface  patch,  and  the  unit  vector  a  represents  the  transverse 
direction  of  the  scattered  field  in  the  lower  medium.  The  result  given  in  Eq.  (2.63) 
(or  Eq.  2.64)  is  remarkable,  as  the  differential  cross  section  is  a  far-field  entity. 
Close  inspection  of  Eq.  (2.63)  shows  that  it  is  the  sonar  equation  when  u{t)  is 
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substituted  for  g{t),  in  which  the  mean  square  scattered  incoherent  field  is 
obtained  by  integrating  the  differential  cross  section  over  the  interface  with 
appropriate  attenuation  due  to  spreading  and  sediment  loss.  Equations  (2.38) 
and  (2.61)  are  a  convenient  way  to  find  0)2  in  theoretical  developments;  Eqs. 
(2.65)  and  (2.66)  are  a  convenient  way  to  find  an. 

2.3.3  Dispersion  of  a  Gaussian  input  pulse 

Two  types  of  pulse  dispersion  occur  in  this  problem.  One  type  is  simply  due 
to  scattering  and  is  treated  by  the  integral  in  Eq.  (2.63)  over  the  scattering 
surface.  The  other  type  of  dispersion  is  due  to  the  frequency  dependence  of  the 
propagation  constant  and  is  embodied  in  Eq.  (2.49).  This  subsection  focuses  on 
the  latter  form  of  dispersion.  Consider  the  case  of  a  pulse  with  a  peak  pressure, 
/7peak,  and  Gaussian  envelope, 


2/2 

u{t)  =  Ppeak  ^  ^  * 

(2.73) 

Its  Fourier  transform  is  given  by 

U{Gy)  —  ^  ^  . 

(2.74) 

The  parameter  is  chosen  to  vary  the  width  of  this  input  pulse,  and  therefore  the 

pulse  bandwidth  and  energy.  Substituting  Eq.  (2.74)  into  Eq. 
Appendix  B)  and  integrating  over  co'  and  co"  yields 

(2.49)  (see 

(2.75) 

where 

q'  =  ^tj+iP, 

(2.76) 

(2.77) 
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and 


P"  =  P*  -i(2Pjq'f,  (2.78) 

Substituting  Eqs.  (2.76)  and  (2.77)  into  Eq.  (2.75),  and  simplifying,  yields  an 
expression  equivalent  to  Eq.  (2.75): 


g  (t)  —  /?peak 


^-exp 


tj-ip* 


-{ipf 


-2Re{(r,^-/P*)r^} 


-{IPf 


(2.79) 


As  mentioned,  Re(r^)  =  Re(r2)  +  fi  represents  the  propagation  delay  of  the  pulse. 
The  imaginary  part  of  is  a  result  of  modeling  the  sediment  as  lossy.  Since 
higher  frequency  signals  are  attenuated  more  than  lower  frequency  signals,  the 
lower  frequency  components  of  a  narrowband  signal  will  be  less  attenuated  than 
the  center  frequency.  The  imaginary  component  of  t2,  along  with  P,  and  the 
imaginary  component  of  P  compensate  for  typically  excessive  loss  in  the  CW 
loss  term  —  the  exponential  term  in  Eq.  (2.50)  or  Eq.  (2.63).  In  addition,  P  is  a 
frequency  dispersion  term  that  results  in  pulse  broadening  of  the  signal  below  the 
rough  interface.  Note  that  for  the  zero  loss  case,  the  parameter  P  =  0.  The 
effect  of  loss  and  the  parameters  P  and  P  is  seen  in  Fig.  2.3  for  from 

Eq.  (2.79)  evaluated  at 


lr=Re(r^) 


tii  =  Im(r^)  =  Im(r2) 


(2.80) 


and  plotted  as  a  function  of  scattered  grazing  angle.  In  Fig.  2.3,  the  pulse  length 
parameter,  t^,  is  equal  to  two  periods  of  the  center,  or  carrier,  frequency  for  both 

cases  shown.  When  the  loss  and  dispersion  are  small,  |g(-^2/)|  approaches 
m(0)  =  1 .  The  smaller  the  scattered  grazing  angle,  the  greater  the  dispersion,  and 
the  smaller  one  would  expect  the  pulse  peak  amplitude.  However,  the  loss  of  the 

CW  signal  represented  by  the  term  jg  greater  than  the  loss  of  the 
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narrowband  signal,  and  the  propagation  distance  becomes  large  for  small 
scattered  grazing  angles  (z  is  fixed).  Even  though  there  is  dispersion,  the  loss 
for  the  narrowband  case  is  significantly  less  than  the  CW  case,  and  |g(-r2,)|>  1  . 

In  the  case  of  Fig.  2.3,  with  depth  of  0.3  m,  |g(-^2t)| 's  greater  than  one  for  very 
small  grazing  angles  since  the  propagation  distance  becomes  very  large.  An 
example  with  zero  loss  always  has  P"=  P*  and  -  0,  and  it  is  always  true  that 

k(o)|^i. 

2.4  Applications  and  calculations 

Although  the  theoretical  results  in  this  dissertation  are  general  and  include 
scattering  back  into  the  water,  the  examples  deal  with  scattering  through  a  rough 
fluid-fluid  interface  (z<0).  Both  plane-wave  and  point  sources  are  considered, 
and  the  mean  square  scattered  incoherent  field  is  computed  for  both  the  flat 
surface  case  and  the  rough  surface  case. 

In  the  examples  of  the  following  section,  the  mean  square  scattered  field  in 
medium  2  is  divided  into  coherent  and  incoherent  parts.  The  incoherent 
component  is  treated  using  first-order  perturbation  theory,  and  the  coherent 
component  is  treated  in  zeroth  order  —  the  flat  surface  solution  is  used.  There  is 
a  slight  inconsistency  in  this  approach,  in  that  the  mean  square  incoherent  field  is 
proportional  to  the  second  power  of  kih,  while  the  square  of  the  coherent  field 
magnitude  is  only  computed  to  an  accuracy  of  zeroth  order.  As  Rice  (1951)  has 
shown,  the  square  of  the  coherent  field  magnitude  to  second  order  in  kih  is 
obtained  by  subtracting  the  power  carried  by  the  mean  square  incoherent  field. 
For  small  kih,  this  correction  is  necessarily  slight.  In  chapter  5,  the  second-order 
perturbation  scattered  field  and  second-order  coherent  field  are  derived. 


34 


2.4.1  Continuous  plane  wave  source 

The  contribution  of  the  evanescent  waves  to  the  field  scattering  through  a 
rough  interface  is  shown  in  the  following  CW  incident  plane  wave  examples.  In 
Fig.  2.4a,  the  first-order  field  penetrating  through  a  rough  interface,  with 
roughness  spectrum  given  in  Eq.  (3.66)  and  =  0.2545,  is  calculated  using  Eqs. 
(3.70)  in  Eq.  (2.50)  and  setting  g(r-r^)  =  lPa.  The  solid  line  in  the  plot 

represents  the  first-order  result  including  evanescent  waves,  and  the  dashed  line 
is  the  first-order  result  excluding  evanescent  waves,  which  is  the  result  in  Eq. 
(2.63).  Since  the  grazing  angle  in  this  example  is  below  critical,  the  zero-order 
field  (dotted  line)  is  evanescent.  Near  the  interface,  the  mean  square  incoherent 
field  decays  with  increasing  depth  at  a  rate  comparable  to  the  square  of  the 
magnitude  of  the  zero-order  field,  but  decays  at  a  rate  dictated  by  the  sediment 
loss  farther  from  the  interface.  Including  the  evanescent  waves  ( K  >  fcj/v)  does 
not  make  a  significant  difference  at  a  distance  greater  than  about  a  wavelength 
(7.5  cm)  below  the  interface.  The  sediment  parameters  used  in  the  plot  in  Fig. 
2.4b  are  the  same  as  Fig.  2.4a,  except  the  roughness  parameters  are  chosen 
such  that  kih  =  \.0.  Perturbation  theory  is  shown  valid  for  kih-\  by  Thorsos 
(1996a)  for  scattering  through  a  rough  interface  with  truncated  power-law 
roughness  spectrum.  In  this  case,  including  the  evanescent  waves  of  the  integral 
( Ar>  iti/v)  does  not  make  a  significant  difference  at  a  distance  greater  than  about 
12cm  below  the  interface. 

These  examples  show  that  including  the  evanescent  waves  in  the  calculation 
of  the  first-order  mean  square  incoherent  field  is  only  necessary  close  to  the 
interface,  and  the  contribution  of  evanescent  waves  to  the  mean  square 
incoherent  field  is  dependent  on  the  interface  roughness.  It  is  important  to  note 
that  in  the  region  below  the  interface  where  the  incoherent  evanescent 
contribution  is  significant,  the  zero-order  coherent  contribution  to  the  mean 
square  field  is  typically  much  greater. 
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2.4.2  Mean  square  incoherent  field  due  to  a  point  source 


Starting  with  the  plane  wave  result  in  Eq.  (2.63)  and  expressing  the  2-D 
integral  in  terms  of  a  summation  over  the  interface  yields 


/2(r,0  = 

n 


hn  ^2n) 


-2Im(fc2)r2„ 

C  » 


(2.81) 


where 


r2n  =  ^-^h 


(2.82) 


is  the  distance  from  surface  patch  n  at  position  of  size  AA  to  the  observed 
field  point,  r,  and 


R-R^ 

R-R' 


(2.83) 


The  propagation  time  from  this  surface  patch  to  the  receiver,  Re(r2n).  is  also 
given  by 


Re(^2n)  =  ^2nlc2  ■ 


(2.84) 


As  the  area  of  each  surface  patch,  AA,  approaches  zero,  the  summation  in  Eq. 
(2.81)  approaches  the  exact  integral  in  Eq.  (2.63). 

An  expression  for  the  first-order  mean  square  field  below  the  interface  due  to 
a  point  source  follows  heuristically  from  Eq.  (2.81).  Here,  the  surface  patch  (area 
AA)  is  chosen  sufficiently  small  in  relation  to  the  distance,  ri„,  from  the  point 
source  to  surface  patch  n,  to  ensure  that  the  field  incident  on  the  interface  is  a 
plane  wave.  However,  the  size  of  this  patch  cannot  be  chosen  arbitrarily  small, 
and  the  distance  between  the  source  and  interface  is  therefore  constrained  to  be 
greater  than  some  minimum  value.  Limitations  of  this  approximation  remain  an 
issue.  The  scattered  or  penetrating  mean  square  incoherent  field  is  given  by 
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/2(r,0  =  AA2— 4^ - -ro8  (t-tdn)^  '' 

n  nnrin 


(2.85) 


where 


a 


in 


K-^s 

R;-R.r 


(2.86) 


R_5  is  the  transverse  coordinate  vector  of  the  source,  g(r)  is  the  pressure  at  1  m 
from  the  source,  and  tq  =  1  m.  The  total  propagation  time  is  given  by  the  real  part 
of 


(2-87) 

with 

n„  =  >;-R/+z7.  (2.88) 

When  the  approximation  g{t)  =  u{t)  is  made,  Eq.  (2.85)  is  in  the  form  of  a 
convolution  with  reduced  computation  time. 

Using  the  expression  for  the  first-order  differential  cross  section  given  in  Eq. 
(3.77),  the  approximate  mean  square  incoherent  field  pulse  in  lossy  sediment 
below  a  rough  interface  due  to  a  Gaussian  narrowband  point  source  above  the 
interface  is  found  for  two  cases.  In  Fig.  2.5a,  the  time-dependent  mean  square 
incoherent  field  at  a  position  below  a  rough  interface  in  a  lossy  medium  due  to  a 
point  source  is  found  using  Eq.  (2.79)  for  g^(r)  in  Eq.  (2.85),  and  is  compared  to 
the  result  obtained  setting  P  and  P  to  zero  —  g{t)  =  u{t).  Although  the  resulting 
pulse  shapes  are  close  in  magnitude,  as  well  as  peak  arrival  time,  using  g{t) 
does  result  in  an  incoherent  pulse  that  has  a  smaller  peak  magnitude.  In  this 
example,  the  parameter  in  Eq.  (2.73)  is  equal  to  two  cycles  of  the  center 
frequency  of  =  20  kHz.  This  frequency  results  in  fej/i  =  0.25,  which  is  within  the 
region  of  accuracy  for  first-order  perturbation  theory  (Thorsos  and  Jackson, 
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1989;  Thorsos,  1990;  Thorsos,  1996a).  In  Fig.  2.5b,  the  point  source  is  farther 
from  the  interface,  but  the  incident  angle  is  the  same.  The  incident  grazing  angle 
is  defined  here  to  be  the  angle  between  the  mean  surface  and  the  line  containing 
the  source  and  receiver  points. 


2.4.3  Zero-order  calculations 


Setting 

'P2_(K)  =  T(K)'Pi(K) 
and 


’  2jcfeiPi(K) 


(2.89) 


(2.90) 


in  Eq.  (2.9)  results  in  an  expression  for  the  field  penetrating  a  flat  fluid-fluid 
interface  due  to  a  point  source  at  with  unit  pressure  magnitude  at  Im 

from  the  source.  The  quantity  U  has  unit  magnitude  and  dimensions  of  pressure 
X  time.  The  resulting  integral  simplifies  to 


X|/2_(r)  =  / J 2^2p2(K))  ^ 

J  ^iPi(K) 


(2.91) 


where  Jq  is  the  Bessel  function  of  order  zero.  For  comparison  purposes,  the 
square  of  the  time-dependent  zero-order,  or  flat-surface  coherent  field 
magnitude, 

|2 


/(r,r)  = 


[ day [/(m)ro \ dK  z*2P2(K)) 

J  J  ^iPi(K) 


(2.92) 


is  plotted  in  Fig.  2.6a  along  with  the  time-dependent  first-order  mean  square 
incoherent  field  using  the  same  parameters  as  in  Fig  2.5a.  In  this  example,  the 
incident  field  is  below  the  critical  angle.  Although  the  zero-order  component  can 
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be  approximated  from  Eq.  (2.91)  using  the  method  of  stationary  phase  to  solve 
the  resulting  integral  (see  for  example  Brekhovskikh,  1980;  Ishimaru,  1991)  when 
the  incident  field  is  above  the  critical  angle,  a  numerical  method  is  more  suitable 
when  the  incident  grazing  angle  is  close  to  critical  (Westwood,  1989).  Figure  2.6a 
shows  an  “exact”  solution  of  the  square  of  the  flat-surface  transmitted  coherent 
field  magnitude,  along  with  the  first-order  mean  square  incoherent  field  from  Fig. 
2.5a.  For  this  example,  the  zero-order  coherent  pulse  arrives  earlier  and  with 
greater  magnitude  than  the  first-order  incoherent  pulse.  In  Fig.  2.6b,  the  source 
is  moved  farther  from  the  surface  than  the  previous  two  examples.  Even  though 
the  sediment  parameters,  receiver  point,  and  incident  angle  (the  angle  between 
the  mean  surface  and  the  line  containing  the  source  and  receiver  points)  are 
unchanged  from  Fig.  2.6a,  the  incoherent  first-order  pulse  is  greater  in  magnitude 
than  the  zero-order  pulse.  As  the  source  is  moved  further  from  the  surface,  the 
stationary  phase  path  for  the  zero-order  pulse  is  longer,  and  the  zero-order  pulse 
amplitude  is  more  attenuated  in  the  lossy  sediment.  Since  Fig.  2.4  is  a  plane 
wave  CW  version  of  this  same  example,  one  would  expect  that  the  first-order 
mean  square  incoherent  field  would  be  greater  than  the  zero-order  field  squared 
for  the  source  point  far  from  the  rough  interface. 

2.5  Summary 

A  general  analytic  expression  for  the  time-dependent  narrowband  mean 
square  incoherent  field  scattered  from  and  through  a  fluid-fluid  rough  interface  is 
expressed  in  terms  of  the  second  moment  of  the  r-matrix,  and  in  terms  of  the 
differential  cross  section.  This  result  is  valid  for  narrowband  source  signals  when 
the  T’-matrix  correlation  function  is  a  slowly  varying  function  of  frequency  over  the 
source  signal  frequencies.  Using  the  exact  CW  incident  plane  wave  result,  along 
with  first-order  perturbation  theory,  and  sediment  parameters  characteristic  of  a 
sandy  seafloor,  we  show  that  the  evanescent  component  of  the  scattered 
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incoherent  field  in  the  sediment  becomes  insignificant  about  a  wavelength  from 
the  surface,  depending  on  the  surface  roughness.  These  examples  show  that 
the  expression  in  terms  of  the  differential  cross  section  is  valid  except  within  one 
or  two  (depending  on  the  roughness  spectrum)  wavelengths  of  the  surface.  This 
result  is  surprising  since  the  differential  cross  section  is  a  far  field  entity. 

The  above  formalism  applied  to  the  scattering  through  a  rough  surface 
shows  the  effect  of  roughness  on  penetration  through  a  surface.  When  the 
grazing  angle  of  a  plane  wave  incident  field  is  below  the  critical  angle  in  relation 
to  the  mean  surface,  the  zero-order  component  of  the  transmitted  field  is 
evanescent,  and  does  not  penetrate  deeply  into  the  seafloor.  Higher  order 
components  contain  downward  traveling  waves,  which  can  increase  the  depth  of 
penetration  of  sound  relative  to  the  flat-surface  case  (Moe  and  Jackson,  1994b; 
Moe  etal.,  1995). 


Figure  2.1  :  Scattering  problem  geometry  for  rough  interface  separating  a 
lossless  fluid  in  medium  1  (z  >  hf{R))  from  a  lossy  fluid  in  medium  2  (z  <  ^/(R))- 
This  diagram  can  be  viewed  as  a  slice  through  a  two-dimensional  surface. 


41 


Figure  2.2  ;  Diagram  describing  transverse  position  variables. 


42 


scattered  grazing  angle  (degrees) 


Figure  2.3  ;  Illustration  of  effect  of  dispersion  on  peak  mean  square  incoherent 
field.  Peak  of  g^(t},  normalized  with  respect  to  max(i/^(r)j  as  a  function  of 

grazing  angle.  5  =  0.0163, /  =  20kHz,  ci  =  1500m  s,  v  =  1.13,  r^  =  100p^. 
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First  order  and  zero  order  field  strength 


Figure  2.4a  :  First-order  and  zero-order  field  strength  as  a  function  of  depth. 
Field  strength  is  given  in  dB  with  respect  to  incident  pressure.  Solid  line 
includes  evanescent  component  in  first-order  calculation,  dashed  line  does  not 
inciude  evanescent  component  in  first-order  calculation.  Incident  angle  is  20°. 
Critical  angle  is  27.75°,  /  =  20kHz,  a  =  0.1m,  y  =  3,  H'2=2xl0"^m,  p  =  2, 
5  =  0.019,  C]  =  1500  m/s,  v  =  1.13,  =  0.254. 
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First  order  and  zero  order  field  strength 


Figure  2.4b  :  First-order  and  zero-order  field  strength  as  a  function  of  depth. 
Field  strength  is  given  in  dB  with  respect  to  incident  pressure.  Solid  line 
includes  evanescent  component  in  first-order  calculation,  dashed  line  does  not 
include  evanescent  component  in  first-order  calculation.  Incident  angle  is  20  . 
Critical  angle  is  27.75°,  /  =  20kHz,  a=0.5m,  y=3,  W2=6.2x10"V,  p  =  1.96, 

6  =  0.0163,  cj  =  1500  m/s,  V  =  1.126, /Ji/i  =  1.00. 
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First-order  mean  square  incoherent  field 


Figure  2.5a  ;  Time-dependent  mean  square  incoherent  field  at  depth  0.2  m  below 
the  mean  interface.  Point  source  height  above  interface  is  0.528  m;  transverse 
distance  between  source  and  receiver:  2m.  Source  frequency  =  20kHz, 
/j  =  100p^.  Solid  line  uses  Eq.  2.85;  dotted  line  is  obtained  from  setting 
g{i)  =  u{t).  Same  sediment  parameters  and  incident  angle  as  Fig.  2.4a. 
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First-order  mean  square  Incoherent  field 


Figure  2.5b  :  Time-dependent  mean  square  incoherent  field  at  depth  0.2m  below 
the  mean  interface.  Point  source  height  above  interface  is  5.5  m;  transverse 
distance  between  source  and  receiver:  2m.  Source  frequency  =  20kHz, 
=  100|xs.  Solid  line  uses  Eq.  2.85;  dotted  line  is  obtained  from  setting 
g(/)  =  u{t).  Same  sediment  parameters  and  incident  angle  as  Figs.  2.4a,  2.5a. 
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First-order  mean  square  incoherent  field  and 
square  of  zero-order  coherent  field  magnitude 


Figure  2.6a  :  Square  of  zero-order  coherent  field  magnitude  shown  with  dotted 
line,  first-order  mean  square  incoherent  field  shown  with  solid  line.  All 
parameters  identical  to  Fig.  2.5a. 
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First-order  mean  square  incoherent  field  and 
square  of  zero-order  coherent  field  magnitude 


x10'® 


Figure  2.6b  :  Square  of  zero-order  coherent  field  magnitude  shown  with  dotted 
line,  first-order  mean  square  incoherent  field  shown  with  solid  line.  All 
parameters  identical  to  Fig.  2.5b, 


CHAPTER  3 


First-order  perturbation  solution  for  rough  fluid-fluid  interface 


3.1  Introduction  to  chapter  3 

In  the  numerical  computations  in  chapter  2,  and  in  the  simulation  in  chapter 
4,  the  water  above  the  seafloor  is  modeled  as  a  homogenous  fluid,  and  the 
sediment  is  modeled  as  a  lossy  homogenous  fluid.  The  rough  interface  between 
the  water  and  sediment  is  modeled  as  a  2-D  random  process  with  a  truncated 
power-law  roughness  spectrum.  The  first-order  fluid-fluid  perturbation  results 
used  in  chapters  2  and  4  are  derived  in  this  chapter.  The  power-law  filter  used  in 
chapters  2  and  4  is  also  discussed  in  this  chapter. 

Perturbation  derivations  are  also  given  in  chapter  5.  However,  the 
derivations  in  chapter  5  allow  the  sediment  bulk  properties  below  the  lowest  point 
on  the  interface  to  be  arbitrarily  stratified,  and  the  approach  allows  evaluation  of 
higher  orders.  The  goal  of  this  chapter  is  to  obtain  the  zero-order  coherent  field 
and  the  first-order  incoherent  field  for  the  fluid-fluid  rough  interface  problem  with 
homogenous  half  spaces.  In  addition  to  derivation  of  an  expression  for  the  field 
scattered  back  into  the  water  as  in  Kuo  (1964)  and  Mourad  and  Jackson  (1989), 
a  solution  for  the  field  scattered  into  the  sediment  is  also  derived.  Although  the 
results  in  this  chapter  are  derived  for  the  two  homogenous  fluid  case,  the  form  of 
the  first-order  scattered  field  solution  derived  in  this  chapter  is  shown  to  be  more 
general  in  chapter  5. 

3.2  Scattering  problem  description 

A  2-D  slice  of  the  3-D  scattering  problem  is  shown  in  Fig.  2.1.  An  arbitrary 
pressure  field,  \ir,(r),  is  incident  on  a  2-D  rough  interface,  z  =  /i/(R),  separating 


the  homogenous  lossless  fluid  in  medium  1  (water)  from  a  homogenous  lossy 
fluid  in  medium  2  (sediment),  resulting  in  the  pressure  field,  \|/y(r),  scattered 

back  into  the  water  and  the  pressure  field,  x|/2-(r),  scattered  into  the  sediment. 
The  mean  of  /(R)  is  taken  to  be  zero  for  convenience  and  without  loss  of 
generality.  Substituting  the  plane  wave  representations  (Devaney  and  Sherman, 
1973)  for  the  field  above  the  interface  (Eqs.  2.5  and  2.6)  and  the  field  below  the 
interface  (Eq.  2.9)  into  the  equation  for  the  continuity  of  pressure  (Eq.  2.3)  yields 


e^KR 


(3.1) 


Substituting  the  plane  wave  representations  for  incident  and  scattered  fields  into 
the  equation  for  continuity  of  normal  velocity  (Eq.  2.4)  yields 


— J  ¥/(K)e'' 

J  d^KK  ■  V(/?/’(R))['F,(K)e"‘^'Pi^*^^^-^^^4  ^f{Ky 

=  _±fA*P2(K)'J^2-(K)e 

J 


^fc,p,(K)/i/(R) 


,iKR 


,ifc,Pi(K)/j/(R) 


jKR 


-/K)t,p2(K)/!f(R)^/KR 


_±  f  j2^K-V(;z/(R))'F2_(K)e 
P2-’ 


(3.2) 


where  the  unit  vector  normal  to  the  interface  in  Eq.  (2.4)  is  given  by 

z-^hfjR) 

"  |£-V;?^(R)|’ 


(3.3) 


and  V  represents  the  gradient  over  the  transverse  coordinates. 
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Since  8«1,  the  expressions  in  Eqs.  3.1  and  3.2  are  valid  (Landau  and  Lifshits, 
1987). 

The  plane  wave  representation  for  the  field  scattered  from  the  interface, 
given  in  Eq.  (2.5),  represents  the  scattered  field  above  the  highest  point  of  the 
interface.  Likewise,  the  plane  wave  representation  for  the  field  scattered  through 
the  interface,  given  in  Eq.  (2.9),  represents  the  scattered  field  below  the  lowest 
point  of  the  interface.  Using  Eqs.  (2.5)  and  (2.9)  to  represent  the  field  on  the 
interface,  z  =  hf{R),  in  Eqs.  (3.1)  and  (3.2)  invokes  the  Rayleigh  hypothesis  (e.g., 
Jackson,  Winebrenner,  and  Ishimaru,  1988). 

The  goal  of  this  chapter  is  to  find  the  scattered  field  in  terms  of  the  incident 
field.  Since 

'?y(K)  =  JATn(K,K')'P/(K').  (3.5) 

and 

'F2_(K)  =  J  ATi2(K,K')'T,(K')  ,  (3.6) 

the  scattering  problem  is  solved  when  the  T-matrices  TJj  and  7j2  are  found  from 
Eqs  (3.1)  and  (3.2). 

3.3  First-order  and  zero-order  perturbation  theory  solution  for  scattered 
field 

The  perturbation  theory  approach  assumes  that  the  product  of  the  wave 
number  in  the  water  and  the  RMS  surface  height,  kih,  is  small.  The  fields  and 
r-matrices  are  expanded  in  a  power  series  in  kih,  and  the  problem  is  solved  one 
order  at  a  time.  In  this  chapter,  we  are  only  interested  in  the  first-order  solution 
used  in  chapters  2  and  4,  and  therefore  keep  only  the  zero-order  terms  and  the 
first-order  terms  of  the  scattered  fields. 
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'F(K)  =  +  (3.7) 

and  the  zero-order  and  first-order  r-matrix  components, 

nK,Ki)  =  T^^\K,Ki)  +  {kih)T^^\K,Ki).  (3.8) 

Likewise,  exponential  terms  In  Eqs.  (3.1)  and  (3.2)  are  also  expanded  in  the  first 
two  terms  of  a  power  series, 

QikM^WW  s  i  +  /jtjPi(K)/i/(R).  (3.9) 

The  Fourier  transform  of  /(R)  is  given  by 

F(K)  =  ^-^fj^^/(R)e“*-*^.  (3-10) 

(27c)^  •’ 

Therefore  /(R)  and  V/(R)  can  be  expressed  as 

/(R)  =  JA'F(K')e*'*^,  (3-11) 

V/(R)  =  J  A'jK'F(K')e*'*^ .  (3.12) 


Using  Eqs.  (3.9)  and  (3.11)  in  Eq.  (3.1)  yields  an  approximate  expression  for  the 
continuity  of  pressure  boundary  condition  appropriate  for  first-order  and  zero- 
order  fields: 

JA[vi/.(K)  +  'Fy(K)]e*-^ 

+J  d^Ki{kih%  (K)[^,(K)(-1)  +  ^/(K)]J 

=  j  A4^2-(K)e'*^’^-  j  •  (3.13) 

Multiplying  both  sides  of  the  above  expression  by  e"'*^  and  integrating  this 
equation  over  R  by  means  of  the  operator  ^d^R,  yields 
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J  A['F,(K)  +  ^y.(K)]6(K  -  K") 

+J  d^Ki(kih%(K)[-^iiK)  +  'P/(K)]J  d^K'F(K')b{K  +  K'-  K") 
=  'A'F2-(K)5(K-K") 

-  ■^/2/s:/(Kiti;i)P2(K)4^2-(K)J d^K'FiK')d{K  +  K'  -  K") , 


(3.14) 


where  6(K)  is  the  2-D  Dirac  delta  function,  integrating  the  above  expression 
over  K  and  over  K'  results  in  the  following  simplified  expression  for  the 
continuity  of  pressure  boundary  condition  appropriate  for  first-order  and  zero- 
order  fields: 


'P,(K")  +  'I7(K") 

+J  A/(A:,/i)Pi(K)[-^,(K)  +  'Fj.(K)]f(K"  -  K) 

=  'F2_(K")  -  J  ci2i:i(Kfci/i)p2(K)'P2-(K)F(K"  -  K) 


Following  the  above  procedure,  we  now  use  Eqs.  (3.9),  (3.11),  and  (3.12)  in  Eq. 
(3.2)  and  obtain  an  approximate  expression  for  the  continuity  of  pressure 
boundary  condition  appropriate  for  first-order  and  zero-order  fields: 

J  Api  (K)[-^,(K)  +  'Fy.(K)]e**^ 

+/J  d^K{kih)^^(Kf¥i{K)  +  'i'/(K)]J  d^K'F(K'W^^'^^''>'^ 

-ij d^K{kih)[¥i(K)  +  'F/(K)]J d^K' 

ki 

=  --  f  i^2/s:Kp2(K)'P2-(K)e'*^ 

.P 

+- J  d^K{kih)K^fil(K)'¥2.(K)j  d^K'FiK'W^^-^^'^'^ 

--  f d^K(kihy¥2-0^)  f  . 

pJ  ■'  ki 

Multiplying  both  sides  of  the  above  expression  by  and  integrating  this 

equation  over  jd^R  yields 
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j  (K)[-^j(K)  +  'Fj.(K)]5(K"  -  K) 

+/ f  d^K  (feift)pf  (K)[t,(K)  +  'Fy  (K)]J  d'^K'FiK'M^  +  K'  -  K") 

-,j  A(t,A)[>P,(K)  +  >P/(K)]j  A'  ■^F(K')5{K  +  K'  -  K”) 

=  --{  AkP2(K)'F2_(K)5(K"  -  K) 

p-* 

+- J  A  (jtj/i)K2p^(K)'F2-(K)j  A'F(K')6(K  +  K'  -  K") 

-^jd^K  {khyi>2- (K)|  d^K'  F(K')5(K  +  K'  -  K")  • 

Integrating  the  above  expression  over  K  and  K'  results  in  the  following 
simplified  expression  for  the  continuity  of  normal  velocity  boundary  condition, 
appropriate  for  first-order  and  zero-order  fields: 

Pl(K")[-'i'/(K")  +  'P/(K")] 

+/[  A(^j/j)p2(K)['F,(K)  +  'Fjc(K)]F(K"-K^ 

-ijd^K{k,h)[¥i(K)  +  -  K) 

Ki 

=  -!-kP2(K")'F2_(K") 

.p 

+i  f d^K (yti/i)K2pi(K)'F2_(K)F(K" -  K) 

p-* 

-i  J  A  (M)^2-(K)  ^  ‘  ^^2  ~ 

3.3.1  Zero-order  field  (flat  surface  reflected  and  transmitted  field) 

Substituting  for  'P2-(K)  and  ^f\K)  for  into  Eqs  (3.15)  and 

(3.18),  and  keeping  only  zero-order  terms,  yields  the  zero-order  equations  for 
continuity  of  pressure  and  continuity  of  normal  velocity: 
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Solving  Eqs.  (3.19)  and  (3.20)  for  and  yields 


(3.19) 

(3.20) 


'F)®\K)  =  r(K)^,(K) 
and 

'F^?(K)  =  T(K)'F,(K), 
where 

l-Kp2(K)/pPi(K) 
^  ’  l  +  Kp2(K)/pPi(K) 


(3.21) 

(3.22) 

(3.23) 


is  the  flat  surface  reflection  coefficient,  and 


T(K)  =  l  +  r(K)  = 


_ 2 _ 

1  +  kP2(K)/pPi(K) 


(3.24) 


is  the  flat  surface  transmission  coefficient.  The  zero-order  field  is  therefore  the 
field  transmitted  and  reflected  from  the  flat  zero-mean  surface.  Comparing  Eqs. 
(3.21)  and  (3.22)  with  Eqs.  (3.5)  and  (3.6)  shows  that  the  reflection  coefficient 

can  be  expressed  in  terms  of  the  zero-order  T'-matrix, 


r(K)  =  r/f  ^(K',  K)5(K'  -  K) ,  (3.25) 

and  the  transmission  coefficient  can  be  expressed  in  terms  of  the  zero-order 
r-matrix, 

T(K)  =  r/2®)(K',K)5(K'  -  K) . 


(3.26) 
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3.3.2  First-order  scattered  field 


Substituting  Eq.  (3.7)  into  Eqs  (3.15)  and  (3.18),  and  keeping  only  first-order 
terms  (terms  with  a  kih  factor),  yields  the  first-order  equations  for  continuity  of 
pressure, 

-h  j  ^2i^/Pi(K)[-4',-(K)  -h  'p)®\k)]f(K"  -  K) 

=  'F^l.^K")- j  A/kP2(K)'p(®\k)F(K"-K),  (3  27) 


and  continuity  of  normal  velocity, 

Kp2(K")  f 

+  (K)['Pi(K)  -h  'Ff  \k)]f(K"  -  K) 

- —  f  d^KWiiK)  +  'f(®)(K)1^^^^^^F(K"  -  K) 

kP2(K")-’  l  '  f  ^  kf 

=  -'f|!.Vk") 

+  Kp2^j^l  ^^^k2p|(K)'f(?(K)F(K"  -  K) 

- i - f  (K)  F(K"  -  K)  . 

kP2(K")-’  " 


3.3.2.1  First-order  field  scattered  from  interface  back  into  water 

Equations  (3.27)  and  (3.28)  are  added  together  to  eliminate  Eqs. 
(3.21),  (3.22),  and  (3.24)  are  used  to  express  4^!°^  and  'F^®^  in  terms  of  'F,-. 
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f  \  Kp2(K"), 


Ip 


Kp2(K") 


I  d^K^j(K)^i(K){l  +  r(K))F(K"  -  K) 


ip 


f  d^K^i(K)(l  +  r(K))^-^^^-^F(K"  -  K) 
J  ki 

J  ^/2/s:k2p^(K)t,(K)(i + r(K))F(K"  -  k) 
f  d^K%(K)(l  +  r(K))^'^^”~-^F(K"  -  K) 

^  hr" 


lcp2(K") 

i 

Kp2(K'') 

i 

kP2(K") 


f  j2/i:/pj(K)^,.(K)(l  -  r(K))F(K"  -  K) 
AzKp2(K)'i'i  (K)(l  +  r(K))F(K"  -  K). 


(3.29) 


Noting  Eq.  (3.5),  and  multiplying  the  above  equation  by  kP2(K")/p.  we  obtain  an 
equivalent  expression  in  terms  of  the  first-order  F-matrix: 


Pl<K")f  1  +  =  |i(l  -  r(K))^22^P,(K) 


pPl(K'') 


1  ,.2r2, 


Kp2(K") 


-i(l  +  r(K))  pf(K)— k^P^(K)  + 

l  p  p 


kP2(K) 


+i(i+r(K)) 


K-(K"-K) 


^1 


1-i 


^F(K"-K). 


For  convenience,  define  H\  \  such  that 
r/J)  (K",  K)  =  //„  (K",  K)F(K"  -  K) . 


Substituting 

kf  kf 

Pf  (K)  -  (KP2(K)f  =  =  1  -  k2 

ki  ki 


(3.30) 


(3.31) 


(3.32) 

(3.33) 
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into  Eq.  (3.30),  and  using  Eq.  (3.31),  yields 


Pi  (K"(l  +  1  (K".  K)  =  i{l  - 

V  PPiv^  ))  P 


+/(i+r(K)) 


K  K 


/A 


(3.34) 


With  the  use  of  Eq.  (3.24),  the  left  side  of  Eq.  (3.34)  is  expressed  in  terms  of  one 
plus  the  reflection  coefficient: 

f/„  (K”,  K) = -(1 + r(K"))(i  - 


(H-r(K"))(l  +  r(K))rKP,(K--)K:P2(K)  ^  ^  1  ^  Ir  » 

iP,(K'')  I  p  P  p'  ^ 


iP,(K'') 

(i+r(K"))(i+r(K)) 

2iP,(K") 


K-K""! 

Ap  ^ 

*^1  J 

(3.35) 


The  first  factor  on  the  right  side  of  the  above  equation  and  the  first  two  factors  in 
the  second  term  are  converted  from  one  plus  the  reflection  coefficient  to  one 
minus  the  reflection  coefficient  through  the  following  equality,  obtained  from  Eq. 
(3.23): 

l  +  r(K)  ppi(K) 


l-r(K)  Kp2(K)' 

The  expression  for  Hu  follows  from  Eq.  (3.35): 

= iipiin 

+b(K",K)(l  -  r(K"))(l  -  r(K))} , 


(3.36) 


(3.37) 
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where 

^  1  ^  Mi"  ■  1? 

a(K",K)  =  l-  — +  --1  (3.38) 

p  Ip  J  K 

and 

b{K",K)  =  Pi(K")Pi(K)(p  - 1) .  (3.39) 

The  first-order  incoherent  field  is  given  by  k\h  multiplied  by  the  2-D  inverse 
Fourier  transform  of  Eq.  (3.5),  where  is  given  by  Eq.  (3.31)  and  Hu  is  given 
by  Eq.  (3.37). 

3.3.2.2  First-order  field  scattered  into  sediment 

The  field  scattered  into  the  sediment  follows  from  rearranging  the  terms  in 
Eq.  (3.27),  and  expressing  the  reflected  and  transmitted  fields  in  terms  of  the 
incident  field  by  using  Eqs.  (3.21)  and  (3.22), 

'F^!](K")  =  t)^)(K")  -  J  d^Ki^i{K){l  -  r(K))F(K"  -  K) 

+J  d^KiK^2iH^  +  r(K))F(K"  -  K).  (3  40) 

As  in  Eq.  (3.31),  define  H12  in  terms  of  the  first-order  T-matrix,  and  2-D 
Fourier  transform  of  the  interface  profile  function  F, 

r/2*\K",K)  =  Hi2(K",K)F(K"  -  K),  (3.41 ) 

where 

Hf^)(K",K)  =  //{|)(K",K)  -  iPi  (K)(l  -  r(K))  +  /kP2(K)(1  +  r(K))  (3.42) 

follows  from  Eqs.  (3.40)  and  (3.41).  Substituting  Eq.  (3.37)  for  Hu  into  the 
above  equation  yields 
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+fc(K",K)(i-r(K"))(i-r(K))} 

-iP,(KXl-  r(K))+  iK|!2(K)(l  +  r(K)).  (3.43) 

The  following  identity  is  shown  by  using  Eq.  (3.36): 

2!pi(K")[-iPi(K)(l  -  r(K))  +  iKP2(K)(l  +  r(K))] 

+i.(K",K)(i-r(K"))(i-r(K)) 

=  -P,(K)P,  (K")(p  - 1)(1  -  r(K))(l  +  r(K")).  (3,44) 

Substituting  Eq.  (3.44)  into  Eq.  (3.43)  yields 

Hi2(K".K)  =  -|^[a(K'',K)(l  +  r(K))-  f.(K",K)(l  -  r(K))],  (3.45) 

2iPi(K  ) 

where  7  is  the  flat  surface  transmission  coefficient  given  in  Eq.  (3.24).  The  first- 
order  approximation  to  the  T-matrix  for  fonward  scattering  into  the  sediment  is 
obtained  from  Eqs.  (3.6),  (3.41),  and  (3.45). 

3.4  Time-dependent  mean  square  incoherent  field  approximation 

Since  the  first-order  T-matrices  are  proportional  to  the  2-D  Fourier  transform 
of  the  zero-mean  interface,  their  expected  value  over  the  rough  interface 
ensemble  is  zero.  In  addition,  the  zero-order  component  of  each  T-matrix  is  only 
dependent  on  the  zeroth  moment  of  the  rough  interface  ensemble.  Therefore 

f,2(K,K,)  =  (/l,*)r»(K,K,)=(*:i/.)«,2(K,K,)F(K-K,)  (3.46) 

and 

T„(K,K,.)  =  (;t,A)7i<;>(K,K,.)  =  {*,'!)Wil(K,K,)f(K  -  K,)  (3.47) 

are  the  first-order  approximations  for  the  incoherent  T-matrices  given  in  Eq. 
(2.18).  Using  the  first-order  approximation  to  the  T-matrix  in  Eq.  (3.41), 
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=  (i,/.fH,2(K',K<,0i)Hr2(K".Kf,O)2XFl(K'-K;)F,(K''-K;')}  (3.48) 

is  the  first-order  approximation  to  the  incoherent  T-matrix  correlation  function, 
given  in  Eq.  (2.33),  for  forward  scattering  through  a  fluid-fluid  interface.  The 
dependence  on  angular  frequency  is  shown  explicitly  in  the  T-matrix  and  Hi2- 
Using  the  identity 

h^{F(K"  -  Ki)F(K'  -  K,))  =  W(K'  -  K,)5(K"  -  K')  (3.49) 

together  with  Eq.  (2.33)  and  Eq.  (3.48)  yields 

^12(K',K",K;,Kf,Q)^  -hQ)',Q)c  +(0") 

=  +a)0H;2(K^K;>,  +co'0w(K-(k;  +K<V2)^  (3.50) 

where  Wis  the  roughness  spectrum.  For  the  case  of  a  wide  sense  stationary 
process  (WSS), 

W(K)  =  f  d^R  C(R)e“*  ^ ,  (3-51 ) 

where 

aRj)  =  (/(R  +  Ri)/(R))  (3.52) 

is  the  surface  autocorrelation  function.  For  this  special  case,  the  identity  in  Eq. 
(3.49)  is  simply  derived  by  making  the  substitutions 

Rj  =R,  -Rj,  R„  ={R,  +R2)/2  (3.53) 

in 

(F(Ki)F(K2))  = -^ 

(271) 


(3.54) 
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and  using  Eqs.  (3.51)  and  (3.52). 

3.4.1  Seafloor  roughness  spectrum 

The  two  dimensional  Gaussian  random  process  /(R)  describing  the  seafloor 
surface  is  assumed  to  be  an  isotropic  power-law  (Fox  and  Hayes,  1985;  Goff, 
and  Jordan,  1989;  Briggs,  1989)  with  a  roughness  spectrum 

W(K)  =  ^.  (3-55) 

Kl 

In  the  notation  of  Mourad  and  Jackson  (1989,  1993),  the  spectral  strength  and 
spectral  exponent  are  represented  by  the  parameters  W2  and  y,  respectively.  For 
2<y<4.  this  surface  random  process  is  an  example  of  a  fractal  (Mandelbrot, 
1983),  and  the  perturbation  theory  requirement  that  k^h  be  small  is  not  satisfied. 
In  fact,  the  RMS  height  is  infinite,  and  the  random  process  is  not  WSS.  However, 
since  this  random  process  is  stationary  in  increments,  Eq.  (3.49)  is  valid 
(Ishimaru,  1978b),  and  the  large  scale  roughness  can  be  removed  by  a  filtering 
operation  (Jackson,  Winebrenner,  Ishimaru,  1986). 

The  power-law  spectrum  is  an  estimate  of  the  roughness  spectrum  over  all 
spatial  frequencies.  The  seafloor  roughness  spectrum  on  scales  over  one  meter 
is  poorly  known,  and  therefore  assuming  a  roughness  spectrum  that  is  a  power- 
law  over  all  spatial  frequencies  (Eq.  3.55)  may  not  be  valid.  Large  scale 
components  of  the  roughness  spectrum  are  set  by  the  spatial  scale  of  the 
experiment.  In  the  case  of  a  tank  experiment,  moderate  to  large  scale  roughness 
will  not  be  present.  Removal  of  the  large  scale  roughness  in  Eq.  (3.55)  is 
therefore  appropriate. 

Let  /(R)  be  the  surface  with  the  power-law  roughness  spectrum  given  in  Eq. 
(3.55),  with  2-D  Fourier  transform,  F(K).  A  filtered  surface,  /(R),  is  obtained  by 
subtracting  a  moving  average  (Papoulis,  1984)  of  /(R)  from  itself: 
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271  a 

/(R)  =  /(R)  -  f J dR'f{^  -  R').  (3.56) 

^00 

Equation  (3.56)  can  be  expressed  as  a  2-D  convolution  of  /(R)  with  a  filter  ^(R), 

/(R)  =  ^|/(R-R>(R'),  (3.57) 

(27C) 


or 


/(R)  =  /(R)05(R),  (3.58) 

where  j(R)  is  of  the  form 

s(R)  =  (27c)^5(R)  -  m(R),  (3.59) 


and 


/r,  x2circ(/?/a) 

mi(R)  =  (27r)^ - 

Tca 


(3.60) 


is  a  moving  average  filter.  The  circ  function  is  defined  as  (Goodman,  1968) 


circ(/?)  = 


fl  R<1 

1 0  otherwise 


(3.61) 


Filtering  the  original  power-law  random  process  with  roughness  spectrum  given 
in  Eq.  (3.57)  with  the  filter  ^(R)  corresponds  to  subtracting  from  each  point  on  the 

random  surface,  /(R),  an  average  of  the  surface  over  a  radius  a  about  each 
point.  The  filter's  Fourier  transform  follows  from  Eq.  (3.59), 

5(K)  =  1-MK)  (3.62) 

where 


M(K)  = 


IJxiKa) 


(3.63) 


Ka 


64 


The  resulting  roughness  spectrum  obtained  by  filtering  the  surface  with  the 
power-law  roughness  spectrum  given  in  Eq.  (3.55)  using  the  filter  in  Eqs.  (3.62) 
and  (3.63)  results  in  the  2-D  surface  random  process  given  by 


Ka  >  ■ 


(3.64) 


Instead  of  subtracting  a  uniformly  averaged  area  around  each  surface  point 
from  the  surface,  subtracting  a  Gaussian  weighted  average  around  each  surface 
point  from  the  surface  results  in  a  roughness  spectrum  without  the  Bessel 
function  oscillations.  In  this  case 


Here,  T  refers  to  the  gamma  function  —  not  the  reflection  coefficient,  and  not  the 
two-frequency  mutual  coherence  function. 


The  roughness  spectrum  given  in  Eq.  (3.66)  is  the  roughness  spectrum  used 
in  this  dissertation.  The  choice  of  the  filter  parameter,  a,  is  discussed  for  each 
numerical  calculation.  The  larger  this  parameter  is,  the  closer  the  filtered 
roughness  spectrum  is  to  the  original  power-law  (Eq.  3.55).  Subtracting  a 
moving  average  of  the  power-law  random  process  from  itself  results  in  removal 
of  low  frequency  components  of  the  original  power-law  roughness  spectrum. 
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Plots  of  this  roughness  spectrum  as  a  function  of  K  for  two  different  values  of  the 
filter  parameter,  a,  are  given  in  Fig.  (3.1). 

3.4.2  Incoherent  T-matrix  correlation  function  approximation  in  Eq.  (2.37) 


In  order  to  evaluate  the  approximation  in  Eq.  (2.37)  using  first-order 
perturbation  theory,  K-K^/2  is  substituted  for  K',  and  K+Kj/2  is  substituted 
for  K"  in  Eq.  (3.50) 

^12(K  -  K^/2,K  +  Kj/2,K<,K;',a)c  +  +  cd") 

=  iti2/fi2(K-K^/2,K;,C0e  +a)'K2(K  +  K^/2,K<',o)c  +co") 


xW(K-(K<  +  K<V2). 

where  is  given  in  Eq.  (2.34)  and  is  repeated  here  for  convenience. 


(3.68) 


k^  =  k;'-k;- 


COc+CO  ^  . 

-  — - -cosO/a,-  - 


c\ 

(o"-co' 

c\ 


(O^  +  CO  -  . 
— = - cos0,a/ 


cos0,d/ , 


(3.69) 


where  0,-  is  the  incident  grazing  angle  and  d,-  is  the  unit  transverse  direction 
vector  of  the  incident  field.  From  Eq.  (2.37),  ^^2  Eq-  (3.68)  is  approximated  by 

^12(K  -  K^/2,K  +  K^/2,K;,Kf,0)^  +  0)',(D^  +  to") 

s^^i2(K,Kj,C0c) 

=  *f|//l2(K,K,-,co,)|V(K-Ki),  (3.70) 


where 


K,-  =  (a)c/q)cos6jd,-. 


(3.71) 
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The  roughness  spectrum,  w(K-(KJ  +  K;^/2),  is  not  a  function  of 
0^  =  0"-©'. 

For  a  narrowband  signal,  it  is  a  slowly  varying  function  of 
0'+0" 

=®C  + - 2 - ’ 


(3.72) 

(3.73) 


except  for  K=  kx.  However,  scattered  energy  with  K>k2  is  evanescent.  From 
inspection  of  Eqs.  (3.38),  (3.39),  and  (3.45),  the  quantity 

//i2(K-Kj/2,K;,0c  +  0')7fi*2(K  +  Kj/2,K;',0,  +  0")  depends  on  0^  only  in  the 
last  term  of  Eq.  (3.38).  In  addition,  the  dependency  on  0^  is  primarily  in  the 
factor  7/(3  outside  the  square  brackets  in  Eq.  (3.45). 


7(K  +  K^/2)  _ _ 2p _ 

Pl(K  +  K^/2)  k(32(K  +  K^/2)  +  pPi(K  +  K^/2) 


Although  the  magnitude  of  Eq.  (3.68)  varies  slowly  with  frequency,  the  phase 
becomes  very  dependent  on  0^  when  either  |K  +  K^/2|  or  |K-K^/2|  is  greater 
than  both  idtj  and  ki-  For  a  narrowband  signal,  this  situation  occurs  for 
evanescent  scattered  waves  and  possibly  scattered  waves  at  shallow  grazing 
angles.  However,  as  shown  in  chapter  2,  including  the  evanescent  scattered 
field  is  only  necessary  within  a  wavelength  or  two  of  the  interface,  depending  on 
the  roughness. 

Figure  3.2a  is  a  plot  of  the  magnitude  of  Eq.  (3.68)  (excluding  the  roughness 
spectrum,  W)  and  Fig.  3.2b  is  a  plot  of  its  phase  as  a  function  of  A"  =  |K|  and 


The  roughness  spectrum  is  not  included  in  these  plots,  because  it  is  not  a 
function  of  fd.  The  direction  of  K  is  the  same  as  the  incident  transverse  wave 
vector  in  this  example.  The  center  frequency  is  20  kHz,  and  the  sediment 


67 


parameters  are  the  same  as  Fig  2.4b.  Figure  3.2a  shows  that  the  magnitude 
becomes  slightly  dependent  on  only  for  shallow  grazing  angles.  Figure  3.2b 
shows  that  the  phase  of  Eq.  (3.68)  becomes  dependent  on  large  fj  only  for  the 
evanescent  scattered  field  and  for  scattering  at  shallow  grazing  angles  —  even 
for  a  relatively  wideband  signal.  The  simplifying  approximation  for  e  in  Eq. 
(2.37)  is  valid  for  all  values  of  K,  except  those  values  corresponding  to 
evanescent  waves  and  those  corresponding  to  scattering  at  small  grazing 
angles.  Because  of  loss  in  the  sediment,  energy  scattered  at  these  small  grazing 
angles  is  insignificant  compared  to  energy  scattered  at  larger  grazing  angles. 
Therefore,  Eq.  (3.68)  can  be  considered  to  be  independent  of  f^. 


3.4.3  First-order  time-dependent  mean  square  incoherent  field 


Substituting  ei2{^^^h<^c)  Eq.  (3.70)  into  Eq.  (2.50)  yields  the  first-order 
expression  for  the  time-dependent  mean  square  incoherent  field  scattered 
through  a  rough  fluid-fluid  interface  due  to  an  incident  plane  wave: 

The  differential  cross  section  from  first-order  perturbation  theory  is  found  by 
substituting  Eq.  (3.70)  into  Eq.  (2.61). 


ai2(a,a,)  = 


_  ^^^5111^(92) 


/fi2(K,K,-,co,)|V(K-K,). 


(3.77) 


Although  the  differential  cross  section  in  this  approximation  is  actually  second 
order  in  kih  (Thorsos  and  Jackson,  1989),  we  refer  to  it  as  the  first-order 
differential  cross  section,  and  will  refer  to  the  corresponding  mean  square 
incoherent  field  in  Eq.  (3.76)  (hirj))  as  the  first-order  mean  square  incoherent 
field.  Another  first-order  expression  for  /2(r,r)  for  an  incident  plane  wave  that 
neglects  evanescent  waves  is  obtained  using  Eq.  (3.77)  in  Eq.  (2.63).  The  effect 
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of  neglecting  evanescent  waves  is  discussed  in  chapter  2.  The  first-order 
differential  cross  section  is  also  substituted  into  Eq.  (2.85)  to  obtain  the  first-order 
expression  for  /2(r,r)  due  to  an  incident  field  from  a  point  source.  This  resulting 
expression  is  used  in  the  simulations  of  experiments  that  are  discussed  in 
chapter  4. 

3.5  Summary 

First-order  perturbation  solutions  of  the  scattered  field  from  and  through  a 
rough  interface  between  two  homogenous  fluids  are  derived  in  this  chapter. 
These  results  are  used  to  show  the  validity  of  approximating  the  incoherent 
T-matrix  correlation  function  as  only  dependent  on  the  center  frequency  of  a 
narrowband  source  signal  (Eq.  2.37).  The  results  derived  in  this  chapter  are 
primarily  used  in  chapters  2  and  4.  The  Gaussian  weighted  filtered  power  law 
(Eq.  3.66)  is  used  for  all  numerical  computations  in  this  dissertation.  The  filter 
parameter,  a,  is  chosen  to  create  a  spectrum  appropriate  for  each  particular 
example.  For  example,  in  the  simulation  of  scattering  into  the  sediment  in  a  tank, 
discussed  in  chapter  4,  a  small  value  is  used  for  this  filter  parameter.  In  ocean 
experiments,  the  value  of  this  parameter  is  considerably  larger. 


fd  (kHz) 


Magnitude 


K  (1/m) 


80 


Figure  3.2a  :  Magnitude  of  Eq.  (3.68),  incoherent  T-matrix  correlation  function  (divided 
by  the  roughness  spectrum),  incident  angle  =  20  ,  fg  =  20kHz,  p  =  1.9608,  5  =  0.0163, 
v=  1.12565,  Cl  =1536  m/s. 
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Figure  3.2b  :  Phase  of  Eq.  (3.68),  incoherent  T-matrix  correlation  function  (divided  by 
the  roughness  spectrum),  same  parameters  as  Figure  3.2a. 


CHAPTER  4 


The  effect  of  roughness  on  acoustic  penetration  of  the  seafioor  as 
given  by  a  fiuid-fiuid  perturbation  model  and  comparison  with  recent 
sediment  penetration  experiments 

4.1  Introduction  to  chapter  4 

Simulation  results  of  a  sediment  penetration  tank  experiment  by  Boyle  and 
Chotiros  (1992)  and  an  ocean  penetration  experiment  by  Chotiros  (1995)  are 
presented  in  this  chapter.  For  the  incoherent  intensity  calculations,  the  first-order 
perturbation  theory  results  derived  in  chapter  3  are  used  in  Eq.  (2.85);  the 
incoherent  intensity  in  the  sediment  is  obtained  by  dividing  Eq.  (2.85)  by  the 
characteristic  impedance  of  the  sediment  (P2C2).  In  the  ocean  experiment 
simulation,  the  time  domain  zero-order  coherent  intensity  is  calculated  using  Eq. 
(2.92)  (divided  by  P2C2).  Note  that  the  incoherent  intensity  calculated  using  first- 
order  perturbation  theory  is  proportional  to  the  second  power  of  kih,  and  the 
coherent  intensity  is  computed  to  an  accuracy  of  only  zero  order.  The 
contribution  of  the  second-order  coherent  intensity  contribution  is  derived  and 
demonstrated  in  chapter  5.  The  effect  of  sound  speed  gradients  in  the  vertical 
direction  on  the  field  scattered  back  into  the  water  is  considered  in  chapters  5 
and  6.  Although  the  influence  of  sound  speed  gradients  on  the  field  scattered 
through  the  interface  is  beyond  the  scope  of  this  dissertation,  possible  effects  are 
mentioned  in  chapter  6. 

4.2  Background 

Figure  4.1  depicts  the  scattering  model  derived  in  chapter  2.  The  source 
pulse  originates  at  a  point  in  the  water  (medium  1)  and  propagates  to  the  «th 
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surface  patch,  its  mean  square  incoherent  field  magnitude  is  multiplied  by  the 
differential  cross  section  of  the  patch,  and  then  it  propagates  to  the  receiver  in 
the  sediment  (medium  2).  In  addition  to  spherical  spreading  from  the  source  to 
the  patch,  and  from  the  patch  to  the  receiver,  there  is  a  sediment  loss  term.  The 
contributions  of  all  patches  are  summed  up  and  divided  by  the  characteristic 
impedance  of  the  sediment,  yielding  the  total  incoherent  intensity  below  the 
rough  interface.  Figure  4.1  includes  the  term  g(t)  derived  in  chapter  2.  We 
argue  later  that  the  baseband  source  signal,  «(r),  can  be  substituted  for  g(r)  for 
both  simulations.  Substituting  u(t)  for  g(t)  in  Eq.  (2.85)  yields:*-'- 


/2(r,0  =  AA^r^ 


2  2 
'in^2n 


(4.1) 


where  is  the  sum  of  the  propagation  time  from  the  source  to  the  «th  patch 
and  the  propagation  time  from  the  nth  patch  to  the  receiver.  Note  that  the 
incoherent  intensity  is  given  by  /2(*’’0/(P2‘^2)>  where  P2C2  is  the  characteristic 
impedance  of  medium  2.  The  first-order  perturbation  expression  for  scattering 
into  the  sediment,  o„,  is  derived  in  chapter  3  and  given  in  Eq.  (3.77).  The  exact 
expression  for  the  time  domain  coherent  intensity  is  given  by  Eq.  (2.92)  divided 
by  the  characteristic  impedance  of  medium  2  (also  see  Fig.  4.1).  Since  we  use 
the  zero-order  field  to  represent  the  coherent  field,  the  transmission  coefficient,  7 
is  substituted  for  ;5’i2  in  this  figure  —  which  then  becomes  Eq.  (2.92).  Following 
the  definition  of  the  coherent  reflection  coefficient  in  Thorsos  (1990),  ^12  is 
defined  in  chapter  5  to  be  the  coherent  transmission  coefficient.  The  power 
spectrum  used  for  both  simulations  is  that  of  a  filtered  power  law  (Eq.  3.66). 
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4.3  Tank  experiment  simulation 

This  simulation  determines  the  incoherent  field  scattered  into  the  sand, 
assuming  a  small  amount  of  roughness  at  the  interface.  A  diagram  of  the  tank 
experiment  is  shown  in  Fig  4.2. 

4.3.1  Simulation  parameters 

A  source  in  water  just  under  one  half  meter  above  sand  transmits  a  signal 
with  Gaussian  envelope  and  60  kHz  carrier  frequency.  The  baseband  source 
signal  used  in  the  simulation  is 

2/2 

u{t)  =  constant  x  ,  (4-2) 

with  r5  =  3a3|xs.  This  signal  approximates  the  time  resolution  achieved  in  the 
tank  experiment.  The  water  sound  speed,  ci  =  1489m/s,  and  signal  center 
frequency  of  60  kHz  correspond  to  a  wavelength  of  X  =  2.5cm  in  the  water.  In 
order  to  agree  with  the  experiment  described  in  Boyle  and  Chotiros  (1992),  the 
value  for  constant  is  chosen  such  that  the  field  Im  from  the  source  is  191pPa. 
The  sound  speed  in  sand  is  given  by  Boyle  and  Chotiros  (1992)  to  be 
C2  =  1675 m/s.  The  loss  parameter  in  the  sand  is  not  given,  and  is  assumed  to  be 
a  value  appropriate  for  sand,  6  =  0.015,  or  0.49dB/m/kHz.  No  roughness 
parameters  for  the  experiment  characterizing  the  rough  interface  for  the  sand  are 
given.  We  assume  only  a  small  amount  of  roughness,  and  set  the  filter 
parameter  in  Eq.  (3.66)  to  fl  =  0.01m.  The  spectral  exponent  and  spectral 
strength  in  Eq.  (3.66)  are  chosen  to  be  y  =  3  and  vv2  =2.17x10“^,  respectively, 
resulting  in  an  RMS  height  (Eq.  3.67)  oih-  1mm  and  k^]^  =  0.25.  This  roughness 
spectrum  is  well  within  the  range  of  validity  for  first-order  perturbation  theory  for 
scattering  through  a  rough  interface  (Thorsos,  1996a).  Of  the  three  roughness 
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parameters  {a,y,w2),  y  =  3  was  chosen  arbitrarily,  but  it  is  consistent  with 
observed  seafloor  roughness  spectra.  Values  for  W2  and  a  in  Eq.  (3.66)  are 
chosen  to  provide  only  a  small  amount  of  roughness  while  yielding  scattering 
intensities  close  to  those  obtained  by  Boyle  and  Chotiros.  The  roughness 
spectrum  is  shown  in  Fig.  4.3. 

4.3.2  Simulation  results 

As  mentioned,  the  baseband  source  signal  u{t)  is  substituted  for  g{t)  in  Eq. 
(2.85)  for  the  simulations  (Eq.  4.1).  The  justification  for  using  this  simpler  model 
in  the  tank  experiment  simulation  is  shown  in  Figs  4.4a  and  4.4b,  where  the  time- 
dependent  mean  square  incoherent  fields  using  u(t)  and  g{t)  are  compared  for 
two  receiver  depths,  with  2  m  transverse  distance  between  source  and  receiver. 
At  the  lower  depth,  the  two  cases  are  almost  indistinguishable:  at  the  depth  of 
0.08cm ,  the  differences  become  more  noticeable,  but  are  still  negligible. 

The  incoherent  intensity  (first-order)  as  a  function  of  time  and  horizontal 
distance  between  source  and  receiver  is  shown  in  Figs  4.5a,  4.5b,  and  4.5c  for 
receiver  depths  of  0.05  m,  0,33  m,  and  0.614  m.  Since  the  baseband  source 
signal  is  Gaussian,  the  arrival  time  is  the  peak  of  the  intensity  pulse.  The  source 
level  is  given  by  Boyle  and  Chotiros  (1992)  to  be  at  191  dB  referenced  tolpPa  at 
1  m  (in  the  sediment).  Figures  4.4a  and  4.4b  are  vertical  slices  of  Figs  4.5a  and 

4.5b,  respectively  (in  Pa^  instead  of  dB).  The  simulation  results  have  been 
presented  in  Figs.  4.5a,  4.5b,  and  4.5c  for  the  purpose  of  direct  comparison  with 
Figs.  5,  6,  and  7  of  Boyle  and  Chotiros  (1992).  By  choosing  the  appropriate 
roughness  parameters,  these  simulations  match  measured  intensity  space  -  time 
behavior  as  presented  by  Boyle  and  Chotiros  (1992),  not  only  in  terms  of  pulse 
arrival  times,  but  also  in  terms  of  the  resulting  field  strength.  Although  the 
simulation  only  contains  the  incoherent  intensity,  the  coherent  contribution, 
dominates  only  for  incident  grazing  angles  above  critical.  The  curve  plotted  in 
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Fig  4.5  is  the  propagation  time  for  the  pulse,  assuming  the  existence  of  a 
1200m/s  refracted  "slow  wave."  The  arrival  time  of  the  incoherent  intensity  pulse 
coincides  with  this  curve.  This  agreement  occurs  for  all  three  depths,  as  well  as 
over  a  wide  range  of  incident  angles.  This  model  does  not  allow  for  slow  wave 
propagation,  and  the  "apparent  slow  wave"  propagation  can  be  accounted  for  by 
considering  the  diffraction  path  for  the  signal.  A  result  due  to  Thorsos  (1995) 
showing  the  relationship  between  an  apparent  slow  wave  and  downward 
diffracted  energy  is  presented  in  the  next  section. 

4.4  How  roughness  can  mimic  a  slow  refracted  wave 

The  purpose  of  this  section  is  to  determine  the  apparent  slow  wave  speed  for 
energy  diffracted  in  a  given  direction.  This  derivation  follows  Thorsos  (1995), 
and  allows  for  the  diffracted  energy  in  the  sediment  to  be  in  directions  ranging 
from  straight  down  to  the  compressional  refracted  direction  (determined  by 
Snell's  law).  Figure  4.6  shows  a  plane  wave  pulse  propagating  with  a  sound 
speed  q  and  at  a  sub-critical  grazing  angle  Gj ,  incident  on  a  water  sediment 
interface.  Two  rays  strike  the  interface;  the  second  intersects  the  interface  at  a 
distance  d  from  the  first,  and  at  a  time  x  later  than  the  first.  Both  rays  are 
scattered  by  the  rough  interface,  and  propagate  into  the  sediment  at  a  grazing 
angle  02  and  at  a  sediment  sound  speed  C2-  The  path  AB,  at  a  grazing  angle 

0^=02 -6',  (4.3) 

represents  the  apparent  slow  wave  refracted  direction  corresponding  to  the  slow 
wave  speed  4.  All  points  perpendicular  to  this  path  are  therefore  taken  to  have 
the  same  total  propagation  time.  From  Fig.  4.6, 

cos  01  =-^  (4-4) 

d 


and 
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cos  02 


(4.5) 


Solving  Eqs.  (4.4)  and  (4.5)  for  x,  and  equating,  results  in  Snell's  law  for  the 
apparent  refracted  slow  wave  direction: 

COSG]  _  €0502  (4  0) 

q  Cl  ' 

Using  Eq.  (4.3)  and 


Cl  —  C2  COS0'  (4.7) 

in  Eq.  (4.6),  and  solving  for  0',  yields 

(4.8) 

The  apparent  refracted  grazing  angle,  02,  follows  from  the  above  result  and  Eq. 
(4.3).  The  apparent  slow  wave  refracted  speed  is  given  by  Eq.  (4.7). 

For  the  same  sound  speed  ratio  used  in  the  simulation.  Fig.  4.7  shows  the 
dependence  of  the  apparent  slow  wave  velocity  on  the  scattered  grazing  angle, 
02.  Figure  4.7  indicates  that  the  apparent  slow  wave  propagation  speed  of 
1200 m/s  in  the  tank  experiment  simulations  (see  Fig.  4.5)  corresponds  to  an 

"average"  scattered  grazing  angle  of  about  80°.  Because  of  loss  in  the  sediment, 
the  average  diffraction  path  is  at  large  grazing  angles  in  a  downward  direction,  as 
other  directions  suffer  larger  absorption  loss. 


/ 

0'  =  tan"^ 

—  COS  01  -00802 

/sin  02 

vq 

/  J 

Taking  the  scattered  grazing  angle  to  be  90°  (straight  down),  we  plot  the 
effect  of  sediment  sound  speed  and  incident  grazing  angle  on  the  apparent  slow 
wave  speed  and  apparent  refracted  angle  in  Fig.  4.8a  and  Fig.  4.8b,  respectively. 

For  this  case,  02  =  90°,  the  apparent  refracted  speed  is  equivalently  expressed 


as 
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^1  +  (C2/Cif  COS^01 


(4.9) 


It  is  evident  from  this  example  that  the  apparent  refracted  slow  wave  speed  and 
the  apparent  refracted  grazing  angle  do  not  depend  strongly  on  the  sediment 
sound  speed,  even  for  these  extreme  sediment  compressional  sound  speed 
examples.  Although  the  apparent  slow  wave  speed  increases  with  increasing 
incident  grazing  angle,  it  is  relatively  constant  for  incidence  below  critical  grazing. 
Below  critical  grazing  is  the  region  of  interest,  since  the  refracted  component 
dominates  above  this  angle. 


Consequently,  roughness  can  mimic  a  slow  refracted  wave  (=1200m/s), 
even  though  the  acoustic  sound  speed  is  considerably  larger  (1675m/s  in  this 
example).  This  tank  experiment  simulation  shows  that  the  effects  of  roughness 
must  be  considered  in  interpretation  of  penetration  experiments. 


4.5  Ocean  experiment  simulation 

The  tank  experiment  determines  the  arrival  time  of  the  pulse,  but  not  the 
direction  that  the  wave  propagates.  In  considering  the  propagation  of  the 
penetrating  acoustic  field,  an  ocean  experiment  using  a  buried  array  of  receivers 
(Panama  City  experiment,  Chotiros,  1991,  1995)  measures  the  direction  and 
speed  of  the  intensity  front.  The  results  of  computer  simulations  of  this 
experiment  using  Eq.  (4.1)  for  the  incoherent  intensity  and  Eq.  (2.92)  for  the 
coherent  intensity  are  presented  in  this  section. 

4.5.1  Experiment  simulation  layout 

The  receiver  locations  used  in  the  ocean  simulation  are  taken  from  Fig.  3.3.1 
of  Chotiros  (1991).  This  array  layout  is  also  shown  in  Fig.  3.4  of  Chotiros  (1995). 
As  in  their  analysis  of  the  experimental  data,  the  receiver  shown  on  the  y  axis 
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and  the  receiver  farthest  out  on  the  x  axis  (see  Fig  3.3.1  of  Chotiros,  1991)  are 
not  used  in  the  simulation.  The  layout  used  in  the  simulation  is  shown  in  Fig.  4.9. 
The  source  is  at  a  height  of  4in  above  the  sediment,  positioned  along  the  x  axis 
at  several  different  distances  from  the  array  corresponding  to  incident  angles  with 
respect  to  the  array  center  of  0,  =  19°,  40°,  24.9°,  117°,  and  9.7°. 

As  in  Chotiros  (1995),  the  center  frequency  for  the  ocean  experiment 
simulations  is  20  kHz.  The  water  sound  speed,  sediment  sound  speed,  and 
sediment  loss  parameter  are  given  by  Chotiros  (1995)  to  be  ci=  1536 m/s, 
C2  =  1729m/s,  and  5  =  0.0163  (.514dB/m/kHz),  respectively,  resulting  in  a  critical 
incident  grazing  angle  of  0,-  =  27.3°.  The  spectral  strength  and  spectral  exponent 
are  consistent  with  Stanic  et  al.  (1988),  which  are  from  measurements  in  the 
vicinity  of  the  experiment  (^2  =6.2x10"^  m  and  y  =  3).  The  filter  truncation 
parameter  (Eq.  3.66),  a  =  0im,  is  chosen  to  be  consistent  with  the  scale  of  the 
experiment,  resulting  in  kih=0.9^  from  Eq.  (3.67).  Note  that  the  transverse 
distance  between  all  receivers  is  within  2a  =  lm.  Thorsos  (1996a)  shows  that 
perturbation  theory  is  valid  for  scattering  through  the  water-sediment  interface 
with  k]h  as  high  as  2.0. 

4.5.2  First-order  intensity  model 

As  in  the  tank  simulation,  the  baseband  source  signal  u{t)  is  also  substituted 
for  g(r)  in  Eq.  (2.85)  for  the  ocean  simulations.  The  baseband  source  signal  is 
taken  to  be 

M(f)  =  e-''/'^  (4.10) 

with  =  lOOpsec.  This  pulse  width  is  comparable  to  that  achieved  in  the  actual 
experiment.  As  mentioned  later,  the  signal  at  each  receiver  is  normalized  such 
that  the  signal  peak  is  the  inverse  of  the  number  of  array  elements.  Therefore, 
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the  magnitude  of  the  source  signal  is  irrelevant  to  the  simulation.  The 
justification  for  using  this  simpler  model  (Eq.  4.1)  for  the  incoherent  intensity  is 
shown  in  Figs  4.10a  and  4.10b,  where  the  mean  square  incoherent  field  on  the 
most  shallow  receiver  (0.05  m  depth)  is  plotted  using  both  the  expression  in  Eq. 
(2.85)  and  Eq.  (4.1),  that  is,  u{t)  is  substituted  for  g(r)  in  Eq.  (2.85).  Although  the 
mean  square  incoherent  field  magnitude  using  the  simplification  g{t)  =  u{t)  is  22% 
greater  than  the  result  obtained  from  Eq.  (2.85)  at  a  depth  of  5cm,  and  17% 
greater  than  the  result  obtained  from  Eq.  (2.85)  at  a  depth  of  18cm,  the  arrival 
times  are  not  affected  by  the  approximation,  g{t)  =  u{t).  Furthermore,  the  error  in 
incoherent  intensity  resulting  from  this  simplification  does  not  affect  the  results 
presented  in  this  chapter. 

4.5.3  Coherent  intensity  model 

Analytical  approximations  for  the  integral  in  Eq.  (2.91)  can  be  found  (for 
example,  Brekhovskikh,  1980;  Ishimaru,  1991).  When  the  incident  field  is  above 
critical,  steepest  descent  solutions  agree  closely  with  the  exact  numerical 
solution.  Furthermore,  for  narrowband  signals,  the  received  zero-order  signal 
amplitude  is  determined  by  the  magnitude  of  the  Green’s  function  evaluated  at 
the  pulse  center  frequency.  Therefore,  the  time  domain  zero-order  received 
pulse  can  be  determined  by  evaluating  the  steepest  descent  solution  at  the  pulse 
center  frequency.  Although  the  approximate  steepest  descent  analytical 
solutions  agree  with  the  numerical  solution  of  Eq.  (2.91)  very  closely  for 
incidence  above  the  critical  angle,  the  solution  near  critical  can  differ  significantly 
(Westwood,  1989).  In  addition,  the  Green’s  function  is  very  frequency  dependent 
when  the  incident  field  is  at  or  below  the  critical  angle.  For  this  reason,  the  zero- 
order  transfer  function  or  Green’s  function  is  found  by  evaluating  Eq.  (2.91) 
numerically  for  a  discrete  range  of  angular  frequencies,  co.  In  order  to  obtain  the 
zero-order  output  pulse  at  a  receiver,  the  discrete  Fourier  transform  of  the  source 
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signal  in  Eq.  (3.21)  is  multiplied  by  the  corresponding  Green’s  function,  and  the 
result  is  transformed  back  to  the  time  domain.  The  resulting  exact  zero-order 
time-dependent  field  magnitude  is  squared  and  added  to  the  time-dependent 
first-order  mean  square  incoherent  field.  The  simpler  steepest  descent  solutions 
have  been  used  for  comparison  purposes. 

Although  the  zero-order  field  will  be  seen  to  be  somewhat  greater  than  the 
actual  coherent  signal  (Rice,  1951;  Thorsos,  1990)  for  the  values  of  k\h  used  in 
these  simulations,  Chotiros'  intensity  velocity  algorithm  is  affected  only  if 
coherent  and  incoherent  intensities  are  close  in  magnitude.  As  shown  in  chapter 
5,  using  only  the  zero-order  field  to  represent  the  coherent  intensity  results  in 
coherent  intensity  calculations  about  10%  greater  than  calculations  using  both  the 
zero-order  and  second-order  contribution.  Including  the  second-order  component 
of  the  coherent  intensity  would  not  produce  noticeably  different  simulation 
results. 

4.5.4  Algorithm  for  estimating  intensity  front  velocity 

Since  the  receiver  locations  in  Chotiros'  ocean  penetration  experiment  are 
not  known  precisely  enough  to  use  a  coherent  spatial  spectral  estimation 
technique,  Chotiros  (1995)  proposed  and  used  the  following  intensity  speed  and 
direction  finding  algorithm.  This  velocity  finding  algorithm  assumes  a  direction 
and  speed  for  a  plane  front  of  propagating  energy,  determines  the  corresponding 
delays  for  each  receiver  with  respect  to  a  reference  receiver,  and  then  adds 
normalized  and  delayed  versions  of  all  receiver  intensities.  If  the  position  of  the 
wth  receiver  is  denoted  by  and  the  reference  receiver  by  rj,  the  relative  signal 
delay  at  the  wth  receiver  with  respect  to  the  reference  receiver  is  given  by 


Ar^  —  ^2  ■  1^2 


(4.11) 
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for  a  plane  wave  traveling  in  a  direction  ^2  speed  C2.  For  a  given  speed  and 
propagation  direction,  the  normalized  received  signals  are  delayed  by  Eq.  (4.11), 
and  then  summed.  The  resulting  magnitude  is  plotted  as  a  function  of  grazing 
angle  (depression  angle)  and  wave  speed  values.  A  peak  in  the  resulting  plot 
occurs  when  the  assumed  direction  and  speed  coincide  with  the  actual  intensity 
front.  As  can  be  deduced  from  section  4.4,  this  type  of  incoherent  processing 
cannot  unambiguously  resolve  wave  velocity.  In  addition,  the  intensity  front 
measurement  will  be  off  if  the  flat  surface  coherent  component  of  the  signal 
dominates  at  the  reference  receiver  but  incoherent  diffracted  energy  dominates 
at  the  other  receivers. 

In  order  to  determine  the  direction  of  the  intensity  front,  an  array  is  generally 
required  to  span  3-D  space  (see  Appendix  D).  However,  the  direction  of  the 
intensity  front  is  assumed  to  have  no  y  component,  and  the  receivers  are 
therefore  only  required  to  span  the  xz  plane. 

4.5.5  Results 

Speed-angle  distributions  for  the  five  source  locations  for  e,  =  79°,40°, 

24.9°,  127°,  and  9.7°  are  shown  in  Figs.  4.11a,  4.11b,  4.11c,  4.1  Id,  and  4.1 1e, 
respectively.  Chotiros  (1991)  states  that  the  receiver  at  z=-0.05m  is  usually 
taken  to  be  the  reference  receiver.  These  plots  are  therefore  obtained  using  the 
receiver  at  z  =  -0.05m  for  the  reference.  Figures  4.11a  and  4.11b  correspond  to 
incidence  above  the  critical  grazing  angle,  and  show  energy  propagating 
approximately  at  the  sediment  sound  speed  and  at  angles  given  approximately 
by  Snell’s  Law.  The  curve  in  the  figures  obeys  Snell's  Law.  Since  the  incident 
grazing  angle  is  above  critical  in  Figs.  4.11a  and  4.11b,  almost  all  of  the  energy 
penetrating  into  the  sediment  is  due  to  the  zero-order  refracted  energy.  The 
coherent  intensity  due  to  the  flat  surface  solution  is  also  greater  than  the 
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incoherent  intensity  in  Fig.  4.11c,  but  the  first-order  incoherent  intensity  is  of 
comparable  magnitude.  At  sub-critical  angles  (Figs  4.1  Id,  4.1 1e),  the  energy 
appears  to  be  propagating  at  about  1400m/s,  at  angles  that  also  agree  with 
Snell’s  Law  for  this  apparent  slow  wave  speed.  This  model  does  not  include 
slow  wave  propagation;  this  apparent  slow  wave  propagation  is  explained  by 
diffractive  scattering.  Equation  (4.6)  shows  that  the  diffracted  intensity  front  will 
always  appear  to  obey  Snell's  Law. 

Although  Chotiros  reports  apparent  slow  wave  propagation  at  speeds 
between  llOOm/s  and  about  1350m/s  (see  Fig.  6,  Chotiros,  1995),  we  observe 
an  apparent  slow  wave  speed  of  1400m/s  in  Figs  4.11c,  4.1  Id,  and  4.1 1e. 
However,  different  choices  of  roughness  spectrum  will  result  in  different  apparent 
slow  wave  speeds.  Although  the  roughness  spectrum  used  in  the  simulation  is 
suitable  for  the  ocean  experiment  location,  burying  receivers  into  the  sediment 
will  likely  alter  the  roughness  spectrum  at  the  interface  above  the  receiver. 
Therefore,  it  is  possible  that  the  roughness  spectrum  used  in  the  simulation  is  not 
an  accurate  representation  of  the  seafloor  at  the  experiment  site.  The  effect  of 
decreasing  the  roughness  by  changing  the  filter  parameter  in  Eq.  3.66  is  shown 
in  section  4.5.7.  In  the  next  section,  another  factor  that  can  affect  the  apparent 
slow  wave  speed  is  discussed. 

4.5.6  Significance  of  reference  receiver  choice 

The  apparent  slow  wave  speed  as  determined  by  Chotiros'  sound  speed  and 
direction  finding  algorithm  can  be  affected  by  reference  receiver  choice. 
Refraction  is  dominant  when  the  incident  angle  is  above  critical,  and  diffraction 
can  be  dominant  at  subcritical  incidence,  depending  on  the  receiver  depth  and 
incident  grazing  angle.  Since  the  refracted  path  is  much  shorter  for  the  receivers 
closest  to  the  interface,  the  refracted  coherent  intensity  will  sometimes  dominate 
over  the  incoherent  intensity  only  at  shallow  receivers.  In  fact,  the  zero-order 
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field  is  negligible  for  the  transverse  distance  of  19.19m  between  source  and 

receiver  (incident  grazing  =  12.7°)  at  all  receiver  depths  except  for  the  depth  of 
0.05  m  —  the  reference  receiver.  Figures  4.12a,  and  4.12b  show  the  mean 
square  incoherent  field  and  square  of  coherent  field  magnitude  (1  Pa  at  1  m  from 
source  in  water)  for  all  receivers  in  the  simulation  as  a  function  of  time  for  the 

12.7°  incidence  grazing  example.  For  all  receivers  except  the  reference  receiver, 
the  mean  square  incoherent  field  is  much  greater  than  the  square  of  the  coherent 
field  magnitude;  at  the  reference  receiver,  the  coherent  field  is  much  greater  than 
the  incoherent  field.  Since  the  refracted  path  is  the  fastest,  the  peak  in  the 
incoherent  intensity  pulse  will  always  arrive  after  the  zero-order  component  of  the 
coherent  pulse.  As  a  result,  the  relative  time  delay  is  greater  than  would  be  the 
case  if  another  receiver  were  chosen  as  the  reference  receiver  and  the  receiver 
signal  at  z  = -0.05m  not  included  in  the  intensity  velocity  algorithm.  The  apparent 
slow  wave  speed  using  this  reference  receiver  is  therefore  slightly  slower  than 
the  apparent  slow  wave  speed  using  any  other  receiver.  For  both  incidence  at 

12.7°  and  9.7°,  the  first-order  incoherent  intensity  is  less  than  the  zero-order 
coherent  intensity  for  the  receiver  at  depth  0.05m:  for  all  other  receivers,  the 

opposite  is  true.  The  outputs  for  12.7°  incidence  and  9.7°  incidence  are  plotted  in 
Figs.  4.13a  and  4.13b  using  the  receiver  at  z=-0.18m  as  the  reference  receiver 
and  removing  the  receiver  at  z  = -0.05m.  These  simulation  results  show  slightly 
faster  apparent  slow  wave  propagation  (1500m/s)  than  the  results  in  Figs.  4.1  Id 
and  4.1 1e.  Note  that  this  speed  is  significantly  smaller  than  the  sediment  sound 
speed  of  1729 m/s. 

Note  that  if  the  apparent  slow  wave  in  the  experiments  is  also  due  to 
diffraction,  then  using  any  other  receiver  as  the  reference  than  the  one  at 
z=-0.05m  in  the  experiment  would  also  result  in  observation  of  an  apparent  slow 
wave  that  is  faster  than  would  otherwise  be  observed.  It  is  important  to  note  that 
the  slower  the  observed  apparent  slow  wave,  the  greater  is  this  discrepancy. 
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4.5.7  Effect  of  the  filter  parameter 

The  purpose  of  this  section  is  to  show  the  effect  of  the  filter  parameter  a  (see 
Eq.  3.66)  on  the  apparent  slow  wave  speed  using  Chotiros'  speed  and  direction 
finding  algorithm.  In  Figs  4.11  and  4.13,  the  filter  parameter  a  =  Oim,  resulting 
in  kih  =  05m,  was  chosen.  Figures  4.14a  and  4.14b  show  speed-angle 

distributions  for  the  two  source  locations  corresponding  to  0,  =12.7°  and  9.7°, 
respectively,  for  a  filter  parameter  value  of  a  =  2.0  m  —  kih  =  \.96.  As  shown  by 
Thorsos  (1996a),  perturbation  theory  is  valid  for  this  case.  As  in  Fig.  4.13,  the 
receiver  at  the  depth  of  0.18  m  is  used  as  the  reference  receiver,  and  the  receiver 
at  a  depth  of  0.05m  is  discarded.  The  differences  between  Figs.  4.14a  and  4.13a 
and  Figs.  4.14b  and  4.13b  are  negligible.  Increasing  the  filter  parameter  beyond 
a  =  Oi  m  to  a  value  more  than  twice  as  large  as  the  scale  of  the  receiving  array 
has  negligible  effect  on  the  speed-angle  distribution  output. 

Figs.  4.15a  and  4.15b  use  a  filter  parameter  a  =  0.02  m,  resulting  in  an 
apparent  slow  wave  speed  of  about  1200  m/s  —  the  same  as  the  tank  experiment 
simulation  results.  Since  Chotiros'  sound  speed  and  direction  finding  algorithm 
measures  the  propagation  of  the  intensity  front,  it  does  not  determine  the  actual 
direction  of  propagation.  The  reason  that  this  algorithm  yields  a  direction  and 
speed  corresponding  to  an  apparent  refracted  slow  wave  follows  from  the 
discussion  in  section  4.4  (Thorsos,  1995). 

4.6  Summary 

By  modeling  sand  as  a  fluid,  and  including  a  small  amount  of  roughness 
/i^  =  l,  the  tank  simulation  results  match  the  acoustical  penetration 
experimental  results  obtained  by  Boyle  and  Chotiros  (1992),  both  in  magnitude 
and  in  arrival  time  of  the  intensity  pulse.  The  simulation  results  of  this  tank 
experiment  show  that  roughness  must  be  included  in  a  sediment  penetration 
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model,  and  that  even  a  small  amount  of  roughness  can  account  for  energy 
penetrating  the  interface  measured  in  their  experiments. 

Although  the  measured  results  in  the  ocean  array  experiments  (Chotiros, 
1995)  yield  an  apparent  slow  wave  that  is  between  50m/s  and  300m/s  slower 
than  the  results  obtained  from  the  simulation,  the  assumed  roughness  spectrum 
may  not  be  an  appropriate  representation  of  the  experimental  site  after  the 
receivers  were  buried.  However,  the  results  of  this  chapter  show  that  diffraction 
results  in  apparent  slow  wave  propagation  that  appears  to  obey  Snell's  Law. 


Intensity  time  series  simulation 


Mean  square  incoherent  field 


n  mn 


7  2/,  ,  \  -2M*2K 
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Square  of  coherent  field  magnitude 
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Figure  4.1  :  Scattering  model  diagram 
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2-D  Isotropic  Roughness  Spectrum 


x10'® 


Figure  4.3  :  Roughness  spectrum  used  in  tank  simulation,  a  =  0.01m,  y  =  3.0, 
h'2  =  2.17  X 10”^  m,h  =  1mm. 
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First-order  mean  square  incoherent  field 


XlO"^ 


arrival  time  (ms) 


use  g(t)  ->  u(t) 
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Figure  4.4a  :  Time-dependent  mean  square  incoherent  field  determined  in  tank 
experiment  simulation;  solid  line  uses  Eq.  (2.85),  dotted  line  uses  Eq.  (4.1). 
Roughness  parameters  are  given  in  Fig.  4.3,  frequency  =  60kHz,  p  =  2.0,  6  =  0.015, 
v=  1.125,  q=  1489 m/s,  4=  33.33 ps.  Source  height  =  0.49m,  receiver  depth 
=  0.08  m ,  transverse  distance  between  source,  receiver  =  2  m . 


91 


First-order  mean  square  incoherent  field 
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Figure  4.4b  :  Time-dependent  mean  square  incoherent  field  determined  in  tank 
experiment  simulation;  solid  line  uses  Eq.  (2.85),  dotted  line  uses  Eq.  (4.1).  Receiver 
depth  =  0.33  m ,  othenvise  same  parameters  as  Fig.  4.4a. 


Receiver  Depth  =  0.08  m 
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Figure  4.5a  :  First-order  incoherent  intensity  as  a  function  of  time  and  transverse 
distance  between  source  and  receiver.  Arrival  time  is  at  pulse  peak;  black  curve 
corresponds  to  1200  m/s  arrival  time  of  apparent  slow  wave.  Same  parameters  as  Figs 
4.3,  and  4.4  -  frequency  =  60kHz,  p  =  2.0,  5  =  0.015,  v  =  1.125,  Cj=1489m/s, 
tg  =  33.33  ps.  Source  height  =  0.49  m ,  receiver  depth  =  0.08  m . 
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Figure  4.5b  ;  First-order  incoherent  intensity  as  a  function  of  time  and  transverse 
distance  between  source  and  receiver.  Receiver  depth  -0.33  m,  otherwise  same 
parameters  as  Fig,  4.5a. 


elapsed  time  (ms) 


Receiver  Depth  =  0.614  m 


Figure  4.5c  :  First-order  incoherent  intensity  as  a  function  of  time  and  transverse 
distance  between  source  and  receiver.  Receiver  depth  =  0,614  m,  otherwise  same 
parameters  as  Fig'  4.5a. 
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Apparent  slow  wave  speed  as  a 
function  of  scattered  grazing  angle 


Figure  4.7  :  Dependence  of  the  apparent  slow  wave  velocity  on  scattered  grazing  angle 
for  incident  grazing  angle,  0,  =  10°. 
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Apparent  refracted  slow  wave 
speed  vs.  incident  grazing  angle 


Incident  grazing  angle  (degrees) 

=  compressional  wave  speed  in  sediment 


Figure  4.8a  :  Dependence  of  the  apparent  slow  wave  velocity  on  incident  grazing  angle 
assuming  average  scattered  field  is  straight  down. 
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Apparent  refracted  grazing 
angle  vs.  incident  grazing  angle 


Incident  grazing  angle  (degrees) 


Figure  4.8b  :  Dependence  of  the  apparent  refracted  grazing  angle  on  incident  grazing 
angle  assuming  average  scattered  field  is  straight  down. 
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Figure  4.9  :  Ocean  experiment  simulation  layout. 
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First-order  mean  square  incoherent 
pressure  at  5  cm  below  interface 


x10'^ 


Figure  4.10a  ;  Time-dependent  mean  square  incoherent  pressure  in  ocean  experiment 
simulation  for  receiver  depth  =  0.05  m.  Solid  line  uses  Eq.  (2.85),  dotted  line  uses  Eq. 
(4.1).  Transverse  distance  between  source  and  receiver  =  19.19m,  a  =  0.5  m,  y  =  3.0, 

W2  =  6.2  X 10“^  m,  /i  =  1.2  cm,  frequency  =  20  kHz,  p  =  1.9608,  6  =  0.0163,  v=  1.1257, 
q  =1536 m/s,  =100|xs. 
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First-order  mean  square  incoherent 
pressure  at  1 8  cm  below  interface 
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Figure  4.10b  ;  Time-dependent  mean  square  incoherent  pressure  in  ocean  experiment 
simulation  for  receiver  depth  =  0.18  m,  othenwise  same  parameters  as  Fig.  4.10a. 
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Incident  Grazing  Angle  =  79.0  (degrees) 
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Figure  4.11a  ;  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Velocity  distribution  for  6^  =79°  (0.43  m 
transverse  distance  between  source  and  reference  receiver),  same  parameters  as  Fig. 
4.10.  White  curve  corresponds  to  energy  propagating  at  speeds  and  directions 
coinciding  with  Snell's  law. 
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Incident  Grazing  Angle  =  40.0  (degrees) 


sound  speed  (m/s) 


Figure  4.11b  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Velocity  distribution  for  =40°  (4.91m 
transverse  distance  between  source  and  reference  receiver),  otherwise  same 
parameters  as  Fig.  4.11a.  White  curve  corresponds  to  energy  propagating  at  speeds 
and  directions  coinciding  with  Snell's  law. 
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Figure  4.11c  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Velocity  distribution  for  0,  =24.9°  (9.22m 
transverse  distance  bet\A/een  source  and  reference  receiver),  otherwise  same 
parameters  as  Fig.  4.11a.  White  curve  corresponds  to  energy  propagating  at  speeds 
and  directions  coinciding  with  Snell's  law. 
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Figure  4.11d  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Velocity  distribution  for  0/  =12.7°  (19.19 
transverse  distance  between  source  and  reference  receiver),  otherwise  same 
parameters  as  Fig.  4.11a.  White  curve  corresponds  to  energy  propagating  at  speeds 
and  directions  coinciding  with  Snell's  law. 


sound  speed  (tos) 


Figure  4.11e  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Velocity  distribution  for  9,-  =9.7°  (25.83  m 
transverse  distance  between  source  and  reference  receiver),  otherwise  same 
parameters  as  Fig.  4.11a.  White  curve  corresponds  to  energy  propagating  at  speeds 
and  directions  coinciding  with  Snell's  law. 


Figure  4.12a  :  Time  dependent  mean  square  incoherent  pressure  in  ocean  simulation 
at  each  receiver  for  the  12.7°  incidence  grazing  (Fig.  4.11d  example). 
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Figure  4.12b  :  Square  of  coherent  pressure  magnitude  at  each  receiver  in  ocean 
simulation  for  the  2.7°  incidence  grazing  (Fig.  4.1 1e  example). 
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Figure  4.13a  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Same  parameters  as  Fig  4.1  Id  except 
receiver  at  z  =  -0.05m  is  removed,  and  receiver  at  0.18  m  depth  (see  Fig.  4.9)  is  used 
as  reference  receiver. 


no 


Figure  4.13b  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Same  parameters  as  Fig  4.1 1e  except 
receiver  at  z  =  -0.05m  is  removed,  and  receiver  at  0.18  m  depth  (see  Fig.  4.9)  is  used 
as  reference  receiver. 
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Figure  4.14a  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Same  parameters  as  Fig  4.13a  except 
roughness  spectrum  filter  parameter  a  =  2.0  m ,  rather  than  a  =  0.5  m. 
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Figure  4.14b  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Same  parameters  as  Fig  4.13b  except 
roughness  spectrum  filter  parameter  a  =  2.0  m ,  rather  than  a  =  0.5  m . 
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Figure  4.15a  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Same  parameters  as  Fig  4.13a  except 
roughness  spectrum  filter  parameter  a  -  0.02  m ,  rather  than  a  =  0.5  m . 
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Figure  4.15b  :  Normalized  output  intensity  of  ocean  simulation  using  intensity  velocity 
finding  algorithm  given  in  Chotiros  (1995).  Same  parameters  as  Fig  4.13b  except 
roughness  spectrum  filter  parameter  a  =  0.02  m,  rather  than  a  =  0.5m. 


CHAPTER  5 


Perturbation  derivation  including  the  effect  of  gradients 


5.1  Introduction  to  chapter  5 

Solution  to  the  first-order  scattered  field  from  and  through  a  rough  interface 
separating  two  homogenous  fluids  is  derived  in  chapter  3.  The  simulations  of  the 
tank  and  ocean  experiments  in  chapter  4  use  first-order  perturbation  theory  to 
calculate  the  incoherent  scattered  intensity,  and  the  flat-surface  (zero-order) 
results  are  used  to  calculate  the  coherent  intensity.  Although  first-order 
perturbation  theory  results  allow  solution  for  the  incoherent  intensity  proportional 
to  the  second  power  of  kih,  the  zero-order  solution  results  in  computation  of  the 

coherent  intensity  only  to  {kihf.  In  order  to  calculate  the  coherent  intensity 

component  proportional  to  [kih)  ,  the  second-order  scattered  field  is  required. 
The  second-order  field  and  second-order  coherent  field  are  calculated  in  this 
chapter,  and  compared  with  the  zero-order  component  of  the  coherent  field. 

Since  the  sediment  in  both  the  tank  and  ocean  experiment  simulations  is 
modeled  as  a  homogenous  fluid,  the  effect  of  sound-speed  gradients  on  the 
scattered  field  is  not  included.  In  this  chapter,  the  sound-speed  gradients  are 
included  in  the  model.  Moe  and  Jackson  (1994a)  take  an  approach  suggested 
by  Winebrenner  etal.  (1995)  that  allows  for  layering  below  the  rough  interface. 
Although  the  medium  just  below  the  rough  interface  is  constrained  to  be  a 
homogenous  lossy  fluid,  the  medium  below  the  lowest  point  on  the  interface  is 
constrained  only  to  be  vertically  stratified  (see  Fig.  5.1),  and  its  complexity  is 
contained  in  the  reflection  coefficient  for  the  corresponding  flat-surface  problem. 
This  approach  is  used  to  derive  recursion  equations  that  allow  for  calculation  of 
the  Nth-order  scattered  field.  The  zero-order  field  and  first-order  scattered  field 
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are  derived  and  compared  with  the  results  in  chapter  3.  The  second-order  field  is 
derived  in  order  to  obtain  an  expression  for  the  second-order  component  of  the 
coherent  field.  The  second-order  coherent  field  is  compared  to  the  flat  surface 
transmitted  field.  The  effect  of  including  sound-speed  gradients  as  well  as  a 
shear  supporting  medium  is  also  discussed. 


5.2  Scattering  problem 


The  seafloor  is  represented  as  a  two  dimensional  random  rough  fluid-fluid 
interface.  The  geometry  of  the  scattering  problem  is  shown  in  Fig.  5.1.  A 
pressure  field,  \i/;(r),  is  incident  on  a  2-D  rough  interface  separating  the 
homogenous  lossless  fluid  in  medium  1  (water)  from  a  homogenous  lossy  fluid  in 
medium  2  (sediment).  Unlike  the  problem  in  chapter  3,  the  region  below  the 
lowest  point  on  the  rough  interface  (z  =-d)  is  arbitrary.  As  in  chapter  3,  the  total 
field  above  the  highest  point  of  the  rough  interface  in  medium  1  is  a  sum  of  an 
arbitrary  incident  field,  \|/i(r),  and  the  resulting  scattered  pressure  field,  \}r/(r), 

and  is  expressed  as 


^iir)  =  jd^K 


(K)  ^  (K)  ^ 


e/KR 


(5.1) 


A  factor  of  is  suppressed.  Because  of  the  flat  interface  between  medium  2 
and  medium  3,  the  field  in  medium  2,  \y2,  is  expressed  as  a  sum  of  plane  waves 
traveling  in  the  positive  z  direction,  as  well  as  a  sum  of  a  plane  waves  traveling 
in  the  negative  z  direction: 


x|/2(r)  =  J  j2^['F2-(K)e"'^2P2(K)z+  'F2+(K)e'^2p2(K)z 


jKR 


(5.2) 


Note  that  in  the  absence  of  upward  refracted  or  reflected  energy,  4^2+  =  0 .  The 
derivation  proceeds  similarly  to  the  derivation  in  chapter  3.  In  Eq.  (5.1),  4^  is 

unknown;  in  Eq.  (5.2),  both  4'2_  and  4^2+  are  unknown.  Three  equations  are 
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required  to  solve  for  these  three  unknowns.  Two  of  the  boundary  conditions  are 
the  continuity  of  pressure  on  the  interface,  given  in  Eq.  (2.2),  and  the  continuity 
of  normal  velocity  on  the  interface,  given  in  Eq.  (2.4).  These  boundary  conditions 
are  used  in  the  scattered  field  derivation  for  the  fluid-fluid  model  in  chapter  3.  An 
additional  equation  is  supplied  by  a  technique  originated  by  Winebrenner  et  al. 
(1995)  in  their  treatment  of  electromagnetic  scattering  from  rough  dielectric 
surfaces.  This  third  equation  uses  the  general  boundary  condition  that  relates 
4^2+  to  'F2-  at  z in  terms  of  the  reflection  coefficient  r23  a\z  =  -d  for  a  wave 
incident  from  medium  2  to  medium  3  (Bucker,  1970); 

r,3(K)=5i±ffi  ^-2ik2M^)d  (5.3) 

't'2-(K) 

This  equation  is  key  to  Winebrenner's  approach.  All  of  the  complexity  of  medium 
3  is  contained  in  r23.  Note  that  media  2  and  3  can  both  be  lossy,  with  a  complex 
wave  number,  but  the  sediment  in  region  3  is  not  constrained  to  be  a  fluid,  e.g.,  it 
could  be  viscoelastic  or  porous.  Furthermore,  arbitrary  gradients  in  density  and 
wave  number  are  allowed  in  medium  3  provided  variations  are  only  in  z. 

The  goal  of  this  chapter  is  to  find  the  scattered  field  in  terms  of  the  incident 
field.  Since 

'Fy(K)  =  I  At„(K,K')'F/(K')  (5.4) 

and 

'I'2-(K)  =  j  ATi2(K,K')'l^i(K')  (5.5) 

(Eq.  3.5  and  Eq.  3.6),  the  scattering  problem  is  solved  when  the  J-matrices  Tn, 
and  Ti2  are  found.  Perturbation  theory  solves  this  problem  by  expanding  the 
fields  and  each  T-matrix  in  terms  of  a  power  series  representation. 
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^(K)=y^^l^^4<W(K)  (5.6) 

IS, 

7-(K,K,)=f'^^r<">(K,Ki),  (5.7) 

n' 

1=0 

and  solving  one  order  at  a  time.  Both  the  factorial  and  the  factor  (kih)'^  do  not 
need  to  be  displayed,  but  this  representation  emphasizes  that  the  perturbation 
expansion  can  be  viewed  as  a  power  series  expansion. 

5.3  Derivation  of  Nth-order  recursion  relation 

Although  only  the  zero,  first,  and  second-order  fields  are  derived  in  this 
chapter,  the  results  of  this  section  can  be  used  as  a  starting  point  to  derive 
higher-order  results.  An  immediate  application  of  higher-order  results  would  be 
to  use  the  second-order  field  component  as  a  correction  to  the  first-order 
coherent  field.  Higher-order  results  serve  as  both  a  correction  and  measure  of 
accuracy  of  lower-order  results. 

5.3.1  Background 

Define  the  Fourier  transform  of  the  «th  power  of  the  surface, 

F„(K)  ^  J  ,  (5.8) 

so  that  F„(K)  is  obtained  from  a  convolution  of  F„_](K)  with  F(K), 

F„(K)  =  F„_i(K)®F(K).  (5.9) 

Also, 

Fo(K)  =  8(K),  (5.10) 


and 
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Fi(K)^F{K),  (5.11) 

is  the  Fourier  transform  of  the  surface  profile  function.  Taking  the  gradient  (Eq. 
3.4)  of 

=  J  AF„+i(K)e*'^  (5.12) 


and  then  dividing  by  «+l  yields 


n  +  1  n  +  1’' 


KF„+,(K)e*“. 


The  following  identity, 

r  (R)V/(R)  =  J  A  KF„+i(K)e'*^*^. 

n  +  n 

follows  from  Eq.  (5.13)  and 

V/"+Vr)=(«+1)/"(R)V/(R). 


(5.13) 


(5.14) 


(5.15) 


5.3.2  Continuity  of  pressure 


Invoking  the  Rayleigh  hypothesis,  we  use  Eqs.  (5.1)  and  (5.2)  in  Eq.  (2.3)  to 
obtain  the  equation  of  continuity  of  pressure: 

=  J  A['P2_(K)e“'^P2(K)¥(R)+  xj/2+(K)e'^p2(KW(R)je'*^R 


Expanding  the  exponential  terms  in  a  power  series. 


(5.16) 


(5.17) 


«=o 


and  using  Eqs.  (5.3),  (5.12),  and  (5.13)  yields 
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«=0  "• 

=  fd^K  X0'P2(K))”^^4^2-(K){[(-1)”  +  r23e'2^2p2(K)^ 


(5.18) 


After  using  the  operator  —^fd^Jie  the  above  expression  is  integrated 

(2Jif  J 

over  R,  and  over  K,  to  obtain 


J A  f;(,p,(K))"  ^^'P,(K)(-l)''  +  >P^(K)]f„(K"  -  K) 


=  jd^K  X(iP2(K))"^^'P2-(K)[(-1)"  +r23e‘"2fc<'‘)^]f„(K''-K)  (5  39) 


Expanding  fields  in  Eq.  (5.19)  in  a  power  series  (Eq.  5.6)  yields 
JA  f;(/pi(K)f^^^i(K)(-l)"F„(K"-K) 

n-0 


= K'f  S  (‘fccK))' 

m=0 


X  £^^'F^(.>(K)((-l)'”  +  r23e'^^2P2(KMj  F^(K"-K) 
J=0  J 


(5.20) 


Combining  powers  of  kih  results  in  the  following  expression  for  the  continuity  of 
pressure  on  the  interface: 
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=  J  A  (-if  (/Pi(K))^iVv(K"  -  KWi(K) 

+jd^K  X  (^5;)(iPi(K)f 

/n=l 

-jA  X  (;Vj(,„p2(K))’”F„(K"-K)[Y  +  (-ir]4'^"^”'(K) 

m=l 

+'I/W(K"), 

where 

Y  =  r23(K)e'^*2p2(KM_ 


5.3.3  Continuity  of  normal  velocity 

An  expression  for  the  equation  for  continuity  of  normal  velocii 
interface  is  obtained  by  substituting  Eq.  (3.3), 

\z-vhnR)\’ 

into  Eq.  (2.4), 

— n  •  V\|ri(r)  =  — n  •  V\i/2(r)|  , 

Pi  P2  ^  ’ 

yielding  an  expression  for  the  continuity  of  normal  velocity, 

1  1 


— n-Vxi/i 

Pi 


lz=/(R)  Pi  L 


^-Vx|/i-V/i/(R) 

dz 


'z=/(R) 


(5.21) 

(5.22) 

on  the 

(5.23) 

(5.24) 

(5.25) 


Using  Eqs.  (5.1)  and  (5.2)  in  Eq.  (5.25)  results  in  the  following  expression: 
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i-J  W(R)+ 

=  —ld^K  /:2P2(K)[-^2-(K)e“'*2p2(W(R)+  'F2+(K)e" 
-yjd^K  K  •  V(;i/(R))['F2-(K)e"‘^2p2(K)¥(R)+  ^2+(K)e'' 


jKR 


ik2h(mm) 


.;kr 


lk2h(K)hf(R) 


iKR 


which  becomes 


-  f  d^K  itiPi(K)  X(iPi(K)/(R))"  f-'FKK)(-l)"  +  ^'^(K) 

Pr  n=0  ”  „ 

-—  \d^KK-  V(;i/(R))  Y  (/pi(K)/(R))"  fT,(K)(-l)"  + 

Pi-"  n=0  ^ 

=  ±U^K  /t2p2(K)  f;(/P2(K)/(R))" ^^[-'I'2-(K)(-1)"  +  'F2+(K)]e*  *^ 

P2-’  n=0 

-—  \d^KK-  V(/i/(R))  y  (/p2(K)/(R)r  i^f'F2_(K)(-l)"  +  'P2+(K)]e*-® 
P2*’  n=0  ^ 


(5.26) 


(5.27) 


with  the  use  of  Eqs.  (5.6)  and  (5.17).  Using  Eqs.  (5.3),  (5.12),  and  (5.14)  in  the 
above  expression  yields 
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-  [Ai:,Pi(K)y(i3,(K))"i^ -(-l)">Pi(K)+'Py(K) 
Pi  „M) 

J  A/i  f;(/pi(K))"  ^^[^i(K)(-ir  +  ^jc(K)] 


n=0 


=  ±|AtjP2(K)i(ip2(K))''^^'H2_(K)[-(-l)"  +  r23e'“2Pi(K¥ 

-■^  J  d^Kh  f;(/p2(K))”  ^^^4^2-(K)[(-1)”  +  ^23  e'2^2p2(KM 

X—  f  A'K  •  . 

n  +  1*' 


1  f  9 

Equation  (5.28)  is  integrated  over  R  using  *',  and  over 

substitution  in  Eq.  (5.22)  is  made,  yielding  the  equivalent  expression, 


^  J  d^K  t,P,(K)  i(ipi<K))''  (^[-(-1)">P,(K)  +  'r/(K)]F„{K"  -  K) 

=  ^J/a:  taPiCK)  i(iP2(K))"  i^>F2_(K)[-(-l)”  +  y]f„(K"  -  K) 
-i-J<i2A:*X(‘P2(K))"^'f'2-(K)[(-l)”+Y]^^^^^f„+l(K''-K). 

n=0 


(5.28) 


K.  The 


(5.29) 


Dividing  the  above  expression  by  ki,  multiplying  by  pi,  and  expanding  in  a  power 
series  —  Eq.  (5.6)  —  yields 
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JA|3,(K)f;{i|3,(K))"fe^  -(-l)'’'P/(K)+i;^^'P}'>(K)  F„(K''-K) 


-iJA  S(‘Pi(K)) 


(n  +  1)! 


K  (K"-K) 


1=0  J  M 

= i  J  Ak|32(K)  f;(;|32(K))" 


+  Y]f»(K"-K) 

/-O 

f;(,p2(K)r 


(5.30) 


Separating  terms,  renumbering  indices  (setting  n  =  l+m),  and  equating  like 
powers  of  kih  yields 
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-jd^K  Pl(K)*  ■P,(K)F;^(K"-K) 

+|A|5,(K)  i  |^^'5^''-"’(K)f»(K''-K) 

m=f)  ■''  -' 


m=0 

.2,,  (-l)'^(iPl(K)f 
N! 

N 


JV!  p,(K)tf 

=  -!-j  AkP2(K)  i  +y]f„(K"-K) 


-—  f (f'K  "V  0*^P2(K))  vy(N-m) 

p-’ 


lK-(K"-K) 


y(^^)(K)[-(-l)"  +  y]  ^  Fm  -  K)- 


kP2(K)^i 

Corresponding  to  the  definition  in  Eq.  (5.22),  define 

7"  =  r23(K")e'^*2p2(K'V 

The  following  definitions, 

^^^(iPl(K))  y^(K---K), 


,p,(K)F(K"-K) 

M^F2(K"-K) 


A(K",K)  = 


(iPl(K)) 

N\ 


N 


-Fn(K"-K) 


(5.31) 


(5.32) 


(5.33) 


(5.34) 
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B+(K",K)  = 


jKp2(K)F(K"-K) 


(/kP2(K)) 


N 


N! 


■Fa,(K"-K) 


B_(K",K) 


-/Kp2(K)F(K"-K) 

H^^^pS-F2(K"-K) 


(-/kP2(K))^ 

A^! 


Fj^(K"-K) 


m(F-l) 
_*i _ 


(iV-1)! 


vi/(^  2) 

Z/ _ 


(A^-2)! 


(5.35) 


(5.36) 


(5.37) 


and 


Y^_(K)  = 


^(A^-1)  xj/(^-2) 

(A^-1)!  (A^-2)! 


...  'f: 


(0) 


2- 


result  in  more  concise  forms  of  Eq.  (5.21), 

wW/ 


N!  '  N\ 


+JA[-yB^(K",K)  +  BL(K",K)]Y2_(K) 
A(-1)^A;v(K",K)^,(K) 

r2, 


AA'(K",K)Yjr(K), 


(5.38) 


(5.39) 


and  Eq.  (5.31) 
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Pl(K") 


iv! 


=  J  A  Pi(K)(-l)^  A;v(K",K)'F,(K) 

-u(N - 1) f  A  ^ (-lfAi^{K",K)^i{K) 

-JAPl(K)A'(K",K)Yjr(K) 

i|  A  Kp2(K)[^5^(K",K)  -  BL(K",K)]Y2_(K) 


+-(y"-1)kP2(K") 
_ifAK.(K:-K) 
p-*  KP2(K)^f 


A^! 


■'^;(K",K)  -  bL(K",K)]Y2_(K) 


(5.40) 


Adding  Eq.  (5.39)  multiplied  by  (y"-1)kP2(K"),  and  Eq.  (5.40)  multiplied  by 
p(y"  +  1),  eliminates  and  results  in  an  expression  for  the  Nth-order 

scattered  field: 
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-[Ky" + I)|i, (k") - (t” -  i)k|32 (K”)] 

=  P(V'' +  l)j  A  p|(K)(-l)%(K",K)'F,(K) 

-p(7"  +  l)J<i^^:p,(K)A'(K",K)yy(K) 

-„(W-l)p(r  +  l)JA5^|||^(-l)%(K",K)>P,(K) 

+p(Y"  +  1)J  A  l^i!^A'(K”,K)YY(K) 

+(y"  +  1)J  A  kP2(K)['^UK",K)  -  b!.(K",K)]Y2_(K) 

-(y"  +  l)\d^K  - BL(K",K)1y2_(K) 

Kp2(K)itf  ^  ■' 

+(Y"-l)KP2(K")}^i^ii:(-l)^A;v(K",K)Y,(K) 

+(y"  -  1)kP2(K")}  d^K  A'(K",K)  Yjr(K) 

-(y"  -  1)kP2(K")J  d^K[iBiiK",K)  +  BL(K",K)]Y2_(K). 


(5.41) 


Define 


pPi(K")fi+Y"'i_A±nE^ 


kP2(K")Ii-y"J  i-r(K")’ 


(5.42) 


where  r  in  second  equality  in  Eq.  (5.42)  is  the  total  reflection  coefficient  of  the 
corresponding  flat  surface  problem  at  the  zero-mean  rough  interface,  obtained 
from  the  well  known  reflection  coefficient  expression  (see  for  example  Sec.  3.3  of 
Brekhovskikh  and  Lysanov). 


^  ri2(K)+r23(K)e^'"^^P^^*^)^ 

l  +  ri2(K)r23(K)e‘2^2p2(K)^’ 


(5.43) 


where 


pPi(K)  +  kP2(K) 


(5.44) 
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is  the  reflection  coefficient  corresponding  to  a  flat  surface  fluid-fluid  interface 
between  medium  1  and  medium  2  in  the  limit  as  the  width  of  medium  2,  d, 
becomes  arbitrarily  large,  or  medium  3  parameters  are  identical  to  medium  2 
parameters.  The  term  r23  is  a  general  reflection  coefficient  corresponding  to  the 
medium  2  -  medium  3  flat  interface.  Although  medium  2  is  constrained  to  be  a 
fluid,  medium  3  is  only  constrained  to  vary  in  the  z  direction.  Dividing  Eq.  (5.41) 
by  (y"-l)KP2  and  using  Eq.  (5.42)  yields 


■jW 
~N\ 
=  jd^K 


^"P,(K)  ^  ^-T^-(K"-K) 


jd^K 

jd^K 

JA 


|},(K")  p,(K")P,(K)^f 

rPl(K)  ,  ,  r'K  (K"-K) 
Pl(K”)  P,(K")Pi(K)t,^ 
^"kP2(K)  ^"K-(K"-K) 


.  pPi(K")  pPi(K")KP2(K)^f  J 

YBUK",K)  +  BL(K",K)]Y2_(K). 


\{-lf  Af^(K",K)%{K) 
A\K",K)Yf(K) 
['iBUK",K)-BL(K",K)]Y2_(K) 


(5.45) 


Dividing  the  above  expression  by  k^,  multiplying  by  pi,  substituting  Eq  (5.33) 
through  Eq.  (5.38)  into  the  above  expression,  using  the  binomial  coefficient 
notation. 


Nl 

^mj  (N-my.m\’ 


(5.46) 


and  collecting  terms  in  the  same  power  in  k\h  yields 
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X  (-l)^(/Pi(K))^F;v(K"-K)'Fi(K) 

+JA 


rPi(K)  ^  rK-(K--K) 


_Pl(K")  pi(K")Pi(K)^fJ 


^^^Pi(K)  ,  .  ^"K-(K"-K) 


_  1  _| - - 

n  /T7"//\0  /T7'\J_Z 


P,(K")Pi(K)tfJ 
X  i  (^)(i3l<K))"F„(K"-K)'p}''-”>(K) 


m=l 


'^^^Kp2(K)  ^^^K-(K^^-K) 

[  PPi(K")  pPi(K")Kp2(K)^f 


(5.47) 


Apply  Eq.  (5.42)  to  Eq.  (5.47) 


i+r(K")  Pi(K)  i-r(K")  i+r(K'')  k  (K"-k) 

2  2  p,(K")P,(K)tf 


Pl(K") 


x(-l) 


(ip,(K))"Fw(K"-K)'P,(K) 
i+r(K")  |3i(K)  i-r(K")  ^  i+r(K")  k-(K"-k) 

2  pi(K")  2  2  pi(K")Pi(K)A:fJ 

<  i  (jflWiCKjrF^CK-'  - 

m=l  ^ 
vjcP-K 


X  i  (ix)(''KP2<K)rF„(K"  -  k)[y  -  (-ir 

m=l 


1  +  r(K")  Kp2  (K)  1  +  r(K")  K  •  -  K) 

2  ppi(K")  2  pPi(K")Kp2(K)/:? 


(N-m), 


N 


+ 


1 

m=l 


A^V-r(K-),, 


2-  (K) 

(/Kp2(K))'”F^(K"  -  K)[y  +  . 


Equation  (5.48)  can  be  further  simplified  by  using  the  identities 


K  (K"-K)  K  K' 


kt 


-l+Pl(K) 


and 


(Pi(K)r-(KP2(K))"  = 


2_kl-K^  kl-K^ 


kf 


=  1-K^ 


(5.48) 


(5.49) 


(5.50) 


to  give 
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x(-I)''(iP,(K))''f;,(K"-K)'ri(K) 


(i+r(K")) 


ij/^(i+r(K")) 


2 

N 


Pl(K")p,(K) 


K  K 


-1-13, (K 


X I  (iX)(;|5l(K)rf„(K''-K)'P<.''-"'>(K) 

m=\ 

l|rfV(i+r(K"))i 


2 

N 


2  K-K" 
^ - r~ 

ki  j 


pPi(K")k132(K)| 

X  i  (^)(/KP2(K)f  F,(K" - K)[t  -  (-1)-]'f(^'”)(K) 

AT 

ij jV(l-r(K'0)S  (5()0KP2(K)rF^(K"-K)[Y  +  (-ir]^^^ 


(5.51) 


The  Nth-order  fields  scattered  away  from  the  interface  and  through  the  interface 
are  found  from  the  recursion  equations,  Eq.  (5.21)  and  Eq.  (5.51).  Each  order 
field  is  found  in  terms  of  all  lower-order  fields.  Although  reciprocity  is  not  obvious 
in  Eqs.  (5.21)  and  (5.51),  nevertheless,  it  should  be  possible  to  show  they  are 
reciprocal  when  each  order  is  calculated. 


5.4  Zero-order 


The  flat  surface  reflected  field  follows  from  Eq.  (5.51)  by  setting  N  =  0. 
=  i  J  A  |8(K''  -  K,T,(K)iLl^ 


■1^-h-Pi(K")Pi(K) 

kf 


Pi(K")Pi(K) 
(i-r(K^O)' 
(i  +  r(K")) 


(5.52) 


'Fj.^^K")  =  I  A6(K"  -  K)r(K")'F,(K) 


(5.53) 


or 
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'Fj°\K)  =  r(K)'P,(K).  (5.54) 

This  result  is  the  same  as  the  zero-order  result  obtained  in  chapter  3,  except  that 
the  reflection  coefficient  in  this  result,  given  by  Eq.  (5.43),  includes  medium  3  in 
the  partial  reflection  coefficient,  r23,  at  the  flat  interface  between  medium  2  and 

medium  3.  The  zero-order  T-matrix,  \  is  the  flat  surface  reflection  coefficient 
multiplied  by  a  delta  function; 

r/f^(K",K)  =  r(K)5(K"  -  K) ,  (5.55) 

where  r  is  given  in  Eq.  (5.43),  representing  the  total  reflection  coefficient  at  z  =  0 
when  the  rough  interface  is  replaced  by  its  mean.  The  zero-order  field 
transmitted  through  the  interface  follows  from  Eq.  (5.21)  by  setting  N  =  0, 

'f(?(K")(y"  + 1)  =  J  -  K)T',(K)  +  ¥p(K")  (5.56) 


or 


'F^?(K)  =  T(K)'F,-(K), 


(5.57) 


where 


T(K)  = 


i+r(K) 

1  + Y 


(5.58) 


is  the  flat  surface  transmission  coefficient.  Note  that  this  transmission  coefficient 
relates  the  down-going  transmitted  field  in  medium  2  to  the  incident  field.  For  the 
case  of  two  homogenous  fluids  —  medium  2  parameters  are  equal  to  medium  3 
parameters,  or  <»,  then  y  =  0,  and 


T(K)  =  l  +  r(K), 


(5.59) 


which  agrees  with  Eq.  (3.24). 
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5.5  First-order 

The  first-order  incoherent  field  scattered  back  into  medium  1  and  the 
downward  traveling  field  scattered  into  medium  2  are  derived  in  this  section  and 
compared  with  the  first-order  fluid-fluid  results  of  chapter  3. 


5.5.1  Derivation  of  T-matrix 


For  the  first-order  field,  set  jV  =  1  in  Eq.  (5.51): 

xj/(l)(K") 


K 


(i+r(K")) 

K  K" 

Pl(K") 

.  k? 

(l+HK”)) 

K  K" 

Pi(K") 

.  kf 

■l  +  Pi(K")Pi(K) 


(i-r(K^O) 


2 

2J  pPi(K")kP2(K)_ 

r(K"))KP2(K)^'(K"  -  K)[y  -  1]'f|?(K). 


(i+r(K")) 

(i-r(K-O)' 


F(K"-K)'P,(K) 


/(O), 


/r  o  I  (i  +  n^ 


^  - 1  -  Pi  (K")P,  -  K)^y>iK) 

k2  _  kP2(K)F(K"  -  K)[y  +  1]'f|?(K) 


k?  J 


The  first-order  F-matrix  then  follows  from  Eq.  (5.4),  the  zero-order  quantities 
given  in  Eqs.  (5.54)  and  (5.57),  and  from  the  above  equation: 
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rW(K",K) 

i/(i  +  r(K"))rK.K"  1  ,  R  rK"iR  ki 

l,r(K)(l  +  r(K''))rK.K"  «  ,K»,»  fK,(i-r(K"))l™”  ki 

1 1  /(i+r(K^o)  r  2  k-k" 

■^2ppi(K")Kp2(K)[’^ 

+|(1  -  r(K"))/Kp2(K)F(K"  -  K)|^(l  +  r(K)). 

By  applying  the  second  equality  in  Eq.  (5.42),  the  above  equation  simplifies  to 
rW  (K",  K)  =  Hi ,  (K",K)F(K"  -  K) ,  (5.62) 

where 

Hi  1  (K", K)  =  * --{Kl  -  nK"))(l  -  r(K))  +  a(l  +  r(K")Xl  +  r(K))},  (5.63) 

2zpi(K  ) 


Kp2(K)F(K"-K)(l  +  r(K)) 


K"  K 


(5.64) 


and 

^;  =  Pi(K")Pi(K)(p-l).  (5.65) 

Reciprocity  is  manifest  here.  Also  note  that  this  result  is  identical  to  the  fluid- 
fluid  first-order  result  in  chapter  3.  In  this  result,  the  complexity  of  medium  3  is 
included  in  the  reflection  coefficient  given  in  Eq.  (5.43).  Although  it  appears  as 
though  Eq.  (5.63)  follows  from  the  two-homogenous  fluid  result  in  Eq.  (3.37),  one 
would  not  typically  be  able  to  generalize  the  two  fluid  result  to  Eq.  (5.63)  simply 
by  substituting  the  general  expression  for  the  reflection  coefficient  given  in  Eq. 
(5.43)  for  the  two-fluid  reflection  coefficient  given  in  Eq.  (3.23).  Using  Eq.  (3.36) 
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in  Eq.  (3.35)  to  eliminate  the  Kp2  was  necessary  to  obtain  the  expression 
identical  to  Eq.  (5.63).  If  this  substitution  were  not  made,  then  the  first-order 
results  in  chapter  3  would  not  generalize  to  the  results  in  this  chapter  —  simply 
by  using  the  general  reflection  coefficient  given  in  Eq.  (5.43). 

The  first-order  field  scattered  through  the  surface  is  found  by  setting  iy  =  1  in 
Eq.(5.21): 

'fW(K")(y"  + 1)  =  -  J  (K)F(K"  -  Kf¥i(K) 

+l  d^Ki^i(K)FiK"  -  Ky¥f\K) 

-J  j2^iKp2(K)F(K"  -  K)[y  -  1]'P|?(K) 

+'p(^(K"). 

^  (5.66) 

Using  the  zero-order  results  in  Eqs.  (5.54),  (5.57),  (5.58),  and  (5.62)  in  the  above 
equation  results  in  the  following  expression  for  the  first-order  field  scattering 
through  the  rough  interface: 

^11.\k")(y"  + 1)  =  j  ^/^F/Pi(K)F(K"  -  K)r(K)T,(K) 

-J  A/pi(K)F(K"  -  K)'Fj(K) 

+J  t/2^/ppi(K)F(K"  -  K)(l  -  r{K))%iK) 

f  i/2/i:Hii(K",K)F(K"  -  K)'Pi(K). 

•'  (5.67) 

The  first-order  F-matrix  follows  from  the  above  equation  and  Eq.  (5.5): 

r/j^(K",K)(Y"  + 1)  =  i{p  -  l)Pi(K)(l  -  r(K))F(K"  -  K) 

-H/7„(K",K)F(K"-K).  (5.68) 


Similar  to  Eq.  (5.62),  define 

r/]^(K",K)  =  //i2(K",K)F(K"  -  K). 


(5.69) 
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If  we  use  the  representation  for  the  transmission  coefficient  in  Eq.  (5.58),  and  the 
definitions  of  a,  and  Z?  in  Eqs.  (5.64)  and  (5.65),  Hu  follows  from  Eq.  (5.68)  and 
its  definition  in  Eq.  (5.69): 

+ r(K))  -  i(i  -  r(K))}.  (5.70) 

This  result  is  a  generalization  of  the  result  in  Eq.  (3.45),  where  the  transmission 
coefficient  relates  the  downward  traveling  field  below  the  zero-mean  rough 
surface  to  the  incident  field  on  the  rough  surface.  For  scattering  back  into 
medium  1 ,  the  complexity  of  medium  3  is  contained  in  the  reflection  coefficient; 
for  scattering  through  the  sediment,  the  complexity  of  medium  3  is  contained  in 
both  the  transmission  and  reflection  coefficients. 

5.5.2  The  effect  of  shear  waves  on  the  differential  cross  section 

As  a  check  of  the  formalism  developed  in  this  chapter,  we  reproduce  a  result 
from  Essen  (1994)  that  includes  shear  supporting  sediment  below  a  rough  fluid- 
fluid  interface.  This  particular  example  is  an  appropriate  check,  because  the 
medium  immediately  below  the  rough  surface  is  a  fluid  and  Essen  uses 
perturbation  theory  to  calculate  the  differential  cross  section. 

The  differential  cross  section  obtained  from  first-order  perturbation  theory  for 
scattering  back  into  the  water  (medium  1)  follows  from  Eq.  (3.77)  by  replacing 
Hx2  from  Eq.  (5.70)  with  Hx\  from  Eq.  (5.63)  and  setting  the  speed  ratio,  v,  to 
one. 

o,l(d^,a,)  =  t,‘'sin2(e2)|i/„(Kj.,K,.)|V(K^  -K,),  (5.71) 

Here,  02  is  the  scattered  grazing  angle  measured  from  the  mean  surface,  and 
Kf  and  K,  are  the  transverse  components  of  the  scattered  field  and  incident 
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field  wave  vectors,  respectively.  This  example  is  concerned  with  the 
backscattering  strength,  lOlog(aii),  where  Kf  =-K,-. 

In  Eq.  (5.63),  the  complexity  of  medium  3  is  characterized  by  the  flat  surface 
reflection  coefficient,  r.  Eq.  (5.43)  is  used  to  find  the  reflection  coefficient  from 
the  partial  reflection  coefficients,  ri2  and  r23.  The  partial  reflection  coefficient 
Ti2  is  the  fluid-fluid  reflection  coefficient  at  the  interface  between  medium  2  and 
medium  3  (z  =  0).  The  partial  reflection  coefficient  r23  is  the  fluid-solid  reflection 

coefficient  at  the  interface  between  medium  2  and  medium  3  at  z  =  -d.  When 
medium  3  is  shear  supporting,  but  gradient  free,  the  partial  reflection  coefficient 
r23  is  given  by  (Brekhovskikh,  1980;  Brekhovskikh  and  Lysanov,  1991) 


r  2in-2i 

(5.72) 

where 

=  Z3^  cos2(2  Y3)  +  Z3,  sin^  (2Y3) , 

(5.73) 

7  P2 
^  KP2’ 

(5.74) 

II 

N 

(5.75) 

and 

II 

N 

(5.76) 

The  parameter  K/  is  the  ratio  of  the  compressional  wave  number  of  medium  3  to 
A:,,  and  the  parameter  is  the  ratio  of  the  shear  wave  number  of  medium  3  to 
A:,.  The  complex  angle  Ys  is  determined  from 
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cos(Y3)  =  p3/(K).  (5-77) 

where 

(5-78) 

In  Fig.  5.2,  the  parameters  given  in  Fig.  5c  of  Essen  (1994)  are  used  to  plot  the 
backscattered  field  strength,  101ogio(on(K^  ,K,)),  K^  =  -K,-,  as  a  function  of 

grazing  angle.  In  this  example,  medium  2  consists  of  a  very  fine  sand  and  has 
thickness  d  =  X.  The  wave  number  and  density  ratios  are  taken  to  be 
K  =  1/(1.12  -  0.004/)  and  p=1.85.  The  sedimentary  rock  in  medium  3  supports 
both  shear  and  compressional  waves  with  wave  number  ratios 
K^  =  1/(2.3-0.004/)  and  =1/(1 .3 -0.il/).  The  ratio  of  the  density  in  the 
sedimentary  rock  (medium  3)  to  that  of  the  fluid  in  medium  2  is  equal  to  2.5.  This 
example  also  uses  the  power-law  roughness  spectrum  given  in  Eq.  (3.55)  or  Eq. 
(3.66)  with  a  =  oo.  Using  a  finite  value  for  the  filter  parameter,  a,  would  only 
affect  the  result  at  larger  grazing  angles  than  those  plotted.  Since  y  =  4.0,  the 
backscattering  strength  does  not  depend  on  frequency.  The  spectral  strength  is 
given  to  be  W2  =0.04/(27t).  This  result  is  identical  to  the  result  obtained  by 
Essen. 

5.5.3  The  effect  of  sound-speed  gradient  on  incoherent  backscattered  field 

In  order  to  illustrate  that  gradients  in  sound  speed  can  have  a  significant 
effect  on  the  intensity  scattered  back  into  the  water,  the  backscattered  field 

strength,  101ogio(aii(%  =-Kp  is  computed  as  a  function  of  grazing 

angle  for  a  sediment  with  a  sound-speed  gradient  in  medium  3,  and  is  contrasted 
with  the  no  sound-speed  gradient  case  —  all  medium  parameters  in  medium  3 
are  identical  to  those  of  medium  2.  In  this  example,  we  have  followed  Mourad 
and  Jackson  (1993)  in  assuming  that  the  square  of  the  wave  number  in  the 
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medium  below  the  rough  seafloor  decreases  linearly  with  depth.  The  choice  of 
linear  decrease  in  squared  wave  number  with  depth, 


where 


^  =  -{z  +  d) 


(5.79) 


(5.80) 


yields  the  well  known  solution  for  the  field  in  medium  3  proportional  to  the  Airy 
function,  Az(t),  and  solution  for  the  reflection  coefficient  in  terms  of  Airy  functions 
(Brekhovskikh,  1980): 


r23(K) 


ik^^2i^)HAi{tQ)  +  AV{t^) 


(5.81) 


H  =  ikU 


(5.82) 


to  =  H\K^-kj).  (5.83) 

and 

X  =  -  ^2  )•  (5-^4) 

The  parameter  a  is  expressed  in  terms  of  the  sound-speed  gradient,  g,  the 
sound-speed  ratio,  v,  and  the  sound  velocity  in  medium  1,  q: 

a  =  2g/(vq),  (5.85) 

where 

g  =  ^£(51  .  (5.86) 

The  parameter  g  is  taken  to  be  the  sound-speed  gradient  at  the  interface 
between  medium  2  and  medium  3.  An  example  that  shows  the  effect  of  a  sound- 
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speed  gradient  on  scattering  back  into  the  water  is  given  below.  In  this  example, 
medium  2  is  taken  to  be  arbitrarily  thin;  the  reflection  coefficient  r  can  be  found 
from  substituting  Eq.  (5.81)  into  Eq.  (5.43),  and  setting  </  =  0. 

_  r^iPi(K)pHA/(roj- At  (/p)  /g  gyv 

ikM^)pHAi{to)  +  Ai'(tQy 

Since  the  medium  above  z  =  -d  must  be  medium  2,  z  =  -d  must  be  at  or  below 
the  lowest  point  on  the  rough  interface.  However  the  sound-speed  gradient  is 
small  enough  in  this  example  that  setting  d  =  0  results  in  a  negligible  difference 
in  scattering  strength.  Since  the  backscattered  field  strength  is  plotted  only  to  a 
grazing  angle  of  70°,  values  of  the  filter  parameter  a>lim  in  the  power-law 
roughness  spectrum  given  in  Eq.  (3.66)  result  in  the  same  scattered  field 
strength. 

Figure  (5.3)  compares  computations  of  seafloor  backscattered  field  strength 
made  with  and  without  a  sound-speed  gradient  for  a  frequency  of  300Hz  and  a 
water  sound  speed  of  1545  m-l.  with  parameters  appropriate  to  a  soft  silty  or 
clayey  bottom  (p2i  =  1.548,  V2i=056,  and  5  =  0.1445x10“^)  and  an  upward 
refracting  sound-speed  gradient  at  the  interface  of  1.3  s“^  This  example  shows 
that  the  sound-speed  gradient  can  cause  an  appreciable  change  in  the  seafloor 
scattering  strength.  The  gradient  in  this  example  is  high  enough,  and  the  loss  is 
low  enough,  so  that  much  of  the  energy  that  is  transmitted  through  the  interface 
is  refracted  upward  to  the  interface  to  strike  it  again,  resulting  in  constructive  or 
destructive  interference.  The  significance  of  this  interference  is  discussed  further 
in  chapter  6.  While  this  example  illustrates  the  effect  of  large-scale  gradients  on 
low  frequency  scattering,  very  steep  gradients  on  centimeter  scales  can  result 
from  biological  reworking  of  sediments  (Richardson,  1986)  and  may  influence 
high  frequency  scattering.  Gradients  in  porosity  are  equivalent  to  gradients  in 
density  (Hamilton,  1974),  which  can  have  strong  effects  on  the  reflection 
coefficient  (Rutherford  and  Hawker,  1978;  Holthusen  and  Vidmar,  1982). 
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Since  the  sound-speed  profile  at  the  experimental  site  is  not  known,  the 
effect  of  gradients  in  sound  speed  on  the  incoherent  and  coherent  intensity 
scattered  through  a  rough  interface  is  not  included  in  the  model  used  in  the 
simulations  (chapter  4). 

5.6  Second-order  field 

For  the  ocean  experiment  simulation  (Chotiros,  1995)  we  use  a  zero-order 
approximation  in  the  coherent  intensity  calculation.  There  is  a  slight 
inconsistency  in  this  approach,  in  that  the  incoherent  intensity  is  proportional  to 
the  second  power  of  kih,  while  the  coherent  intensity  is  only  computed  to  an 
accuracy  of  ifej/i  to  the  zeroth  power.  The  second-order  field  is  found  in  the  next 
section.  Using  this  result,  we  obtain  an  expression  for  the  second-order  coherent 
field. 

5.6.1  Derivation  of  T-matrix 

The  procedure  for  finding  the  second-order  field  is  similar  to  finding  the  first- 
order  field.  However,  both  first-order  and  zero-order  results  are  required  in  order 
to  obtain  the  second-order  field.  Setting  N=2  in  Eq.  (5.51),  substituting  K'  for  K 
in  the  terms,  and  using  Eqs.  (5.54),  (5.57),  (5.58),  and  (5.69)  and 

F2  (K"  -  K)  =  J  A'F(K"  -  K')F(K'  -  K)  (5.88) 


from  Eq.  (5.9)  yields 
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Tvf\K")  =  J  A'  I  d^KFiK"  -  K')F{K'  -  K)'F,(K) 
^(l.nK”))^[-5|:.l^P.(K'0P.(K)|^]r(K) 


,(1-y)kP2(K)  ,7  K  K 


.(i.r(K'0)(i.r(K))|^^^K- 

-(i-r(K"))(i+r(K))(KP2(K)f} 

+2i  j  d^K'ld^KFiK"  -  K')F(K'  -  K)'J^,(K) 


X  Hn(K',K)(l  +  r(K")) 


1  K'K' 


Pi(K")  k{ 


l-Pl(K 


+H,2(K',K)(l  +  r(K"))(l  +  r)^^[K^-^ 
-Hl2(K',K)(l-r(K"))KP2(K')(l-Y')}- 


(5.89) 


Define  Gii(K"K',K),  such  that 

=  J  A' J  Af(K"  -  K')F(K'  -  K)'F,(K)Gii(K",K',K).  (5.90) 

With  much  work  and  use  of  the  second  identity  in  Eq.  (5.42),  the  above 
expression  is  equivalent  to 
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Q 1  (K''',  K',K) 


(1 + r(K"))(i  -  r<K))|^  -  {1  -  r(K"))(i + r(K)) 

Pl(»-  ) 


+ 


iHn(K',K) 


Pi(K") 


K'  K" 


-1 


(1 + r(K"))  -  Pi(K")Pi  (K')(1  -  r(K")) 


+/Hl2(K',K)(l  +  r(K")) 


(l±ll 

Pl(K") 


p  p  A:] 


1 


(i+r(K"))i+r(K') 


In  order  to  find  the  second-order  field  scattering  through  the  interface,  set 
in  Eq.  (5.21). 

^(2)  (K")(l  +  r)  =  -jd^K  (Pi  (K)f  F2  (K"  -  K)'Pi  (K) 

+j  d^K  2/Pi  (K)F(K"  -  KwfkK) 

-]  A(Pi(K)f  F2(K"  -  KfvfkK) 

+1  AkP2(K)F(K"  -  K)[l  -  y]'f|1.^(K) 

+ J  d^K  (kP2  (K))^  F2  (K"  -  K)[1  +  (K) 

+¥p(K"). 

The  expressions 

=  j  A'H„(K,K')F(K-K')'FKK') 
and 

^i*)(K)  =  J  A'//i2(K,K')F(K  -  K')^,(K') 


(5.91) 

N  =  2 

(5.92) 

(5.93) 

(5.94) 


are  obtained  by  combining  Eqs.  (5.4),  (5.62)  and  Eqs.  (5.5),  (5.69).  Combining 
Eqs.  (5.57)  and  (5.58)  yields 
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>J,(0)(K)  =  l  +  r(K)y(K)^ 

(1  +  7) 

Using  Eqs  (5.88),  (5.93),  (5.94),  and  (5.95)  in  Eq.  (5.92)  yields 

'fP^(K")(i+y") 

=  J  d^K\<P-K'F{K''  -  K')F(K'  -  K)'F,(K) 

x{-p?(K)(l  +  r(K))  +  (KP2(K)f  (1  +  r(K))  +  2i//H(K',K)pi(K') 

+2iHi2(K',K)KP2(K')(l  -  Y)  +  i(K",K',K)}. 

Using  the  following  definition  for  Gi2(K",K',K), 

»j/(2)(K")  =  J  A  J  A'F(K"  -  K')F(K'  -  K)'F,(K')Gi2(K",K',K), 
in  Eq.  (5.96)  yields 

Gi2(K",K',K)(1  +  Y")  =  {-pf  (K)(l  +  r(K))  +  (Kp2(K)f  (1  +  r(K)) 

+2///n(K',K)Pi(K') 

+2///i2(K',K)kP2(K')(1  -  Y')  +  Gii(K",K',K)} 
By  noting  the  following  identities, 

(Pl(K)f  -  (KP2(K)f  =  =  1-k\ 

ki  ki 

Eq.  (5.98)  simplifies  to 

(1  +  y")Gi2(K",K',K)  =  (k2  -  l)(l  +  r(K))  +  2///n(K',K)Pi (K') 

+2///i2(K',K)kP2(K')(1  -  Y')  +  Gn(K",K',K). 


(5.95) 

(5.96) 

(5.97) 

(5.98) 

(5.99) 

(5.100) 
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5.6.2  Second-order  coherent  field 

In  order  to  check  the  zero-order  approximation  for  the  coherent  intensity,  the 
goal  in  this  section  is  to  find  the  coherent  component  of  the  second-order  T- 

matrix,  If  we  allow  G  to  represent  both  in  Eq.  (5.91)  and  G12 

in  Eq.  (5.98),  it  can  be  expressed  as  the  sum: 

G(K",K',K)  =  G^(K",K)  +  Gfo(K",K',K).  (5.101) 

The  T-matrix  is  found  from  Eq.  (5.90)  or  Eq.  (5.97): 

T^^\k",K)  =  J  d^K'F{K"  -  K')F{K'  -  K)G^(K",K) 


+1  A'F(K"  -  K')F(K'  -  K)Gb  (K",  K',  K) . 

(5.102) 

Substituting 

/,(R)  =  /(R) 

(5.103) 

and 

/2(R)  =  /(R)e*'‘' 

(5.104) 

into  the  identity 

,  '  ,  fd2R/,(R)/j(R)e-*»  =p2/i:'Fi(K-K')F2(K') 

{iKf  •'  ^ 

(5.105) 

yields 

/2(R)e"'(*^"“^^  ’^  -  J d^K'F{K"  -  K')F(K'  -  K). 

(5.106) 

Note  that  /  is  a  WSS  process  with  unit  variance.  Therefore,  taking  the  expected 
value  of  Eq.  (5.106)  over  the  ensemble  of  interfaces  results  in  the  following 
identity: 
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6(K"  -  K)  =  (J  d^K'Fx  (K"  -  K')F2(K'  -  K)).  (5. 1 07) 

Consequently,  the  expected  value  of  the  first  term  on  the  right  side  of  (5.81) 
results  in  the  following  expression  for  the  coherent  second-order  r-matrix: 

(I  d^K'FiK"  -  K')F(K'  -  K))g«(K",K)  =  5(K"  -  K)G^(K,K) .  (5.1 08) 

Define  Int2  to  be  the  second  term  on  the  right  side  of  Eq.  (102). 

Int2  =  J  A'F(K"  -  K')F(K'  -  K)Gf,(K",K',K).  (5. 1 09) 


Substituting 


m)= 


(5.110) 


for  both  F(K"-K')  and  F(K'-K)  in  Eq.  (5.109),  and  taking  the  expected  value, 
yields 


(Int2)  =  J ^/2F'Gfc(K",K',K)|j d^R 

j  jV  C(R'  -  R)//i^} , 


(5.111) 


where 

C(R'-R)sA^(/(R')/(R)) 


(5.112) 


is  the  surface  autocorrelation  function.  By  noting 


1 


d^R’e 


-i{K'-K'>R' 


C(R'-R)  =  e-'('‘"-'‘')* 

(2r)2J 

=  e-'(K"-K')Kvf(K"-K')/A2 


(5.113) 


and 
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f  d^R  e-'(K'-K)  Rg-i(K"-K')  R  _  g. A  ^  g/ _  K")iy(K"  -  K') ,  (5.11 4) 

Eq.  (5.111)  simplifies  to 

(Int2)  =  5(K  -  K") j d^K'Gt{K",K\K)W{K"  -  K')/h^ .  (5.1 1 5) 

By  using  Eqs.  (5.107)  and  (5.115),  the  expected  value  of  Eq.  (5.102)  simplifies  to 
^7’(2)(k",K)^ 

=  6(K"  -  K){g^(K,K)  +  J  dh'Gb{K,K',K)W{K  -  K')/h^].  (5.1 16) 

By  making  an  approximation  suggested  by  Thorsos  and  Jackson  (private 
communication),  Eq.  (5.116)  is  approximately 

^ j(2)(k"  K))  s  5(K"  - K){g^(K,K)  +  Gfc(K,K,K) j d^K'  W(K  -  K')/ h^} 


=  5(K"  -  K){GaiK,K)  +  G^(K,K,K)}.  (5.1 17) 

To  second-order,  the  expected  value  of  the  T-matrix  is  given  by 

(r(K^K))  =  ^r(®^(K",K)^  +  ^(A:i/if(7’^^^(K'',K)y  (5.118) 

The  intensity  reflection  or  transmission  coefficient,  ;5’/,  (Thorsos,  1990)  follows 
from  Eq.  (5.118). 

=|;?(K)p  =|;5o(K)P  +  (Mf  ^^o(K)^2(K)^  (5-119) 

where 

;?(K)  =  ;?o(K)  +  ^(^i/*f^2(K),  (5-120) 

(r(K",K))  =  5(K"-K);?(K),  (5.121) 


and  (see  also  Berman,  1992) 
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^7’(2)(k",K)^  =  5(K"  -  K);?2(K).  (5. 1 22) 

For  the  case  of  scattering  from  an  interface,  substitute  the  flat-surface  reflection 
coefficient  r(K)  for  ;?o(K)  and  ;?ii(K)  for  ;5'(K);  for  the  case  of  scattering 
through  an  interface,  substitute  the  flat  surface  transmission  coefficient  7(K)  for 
;?o(K),  and  ^i2(K)  for  ^(K). 

For  scattering  from  the  interface,  Eq.  (5.91)  along  with  the  definition  in  Eq. 
(5.101)  yields 


Gii«(K",K) 

=i(ic2 -i)  (i+r(K"))(i-r(K))|^-(i-r(K"))(i+r(K)) 


(5.123) 


and 


Gnf,(K",K',K) 


Pi(K") 


K'  K' 


-1 


V 


(1  -h  r(K"))  -  Pi(K")Pi  (K')(1  -  r(K")) 


+iHi2(K',K) 


(l+ll 


X 


Pl(K") 


p  p 


1  J 


(1 + r(K"))  -  pP,(K")Pi(K')(i  -  r(K")) 


i-r(K0 

i-t-r(K0 


(5.124) 


Likewise,  for  scattering  through  the  interface,  using  Eq.  (5.100)  and  the  definition 
in  Eq.  (5.101)  yields 


G,2,(K",K)  =  l)(l  +  r(K))  +  ^G,  ,„(K",K) 


(5.125) 


and 
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Gi2z,(K",K',K)(1  +  y") 

=  2/H„(K',K)pi(K')  +  2/i/i2(K',K)KP2(K')(l - Y)  ^26) 

+Gn,(K",K',K).  ^  ^ 

5.6.3  Second-order  coherent  field  calculations  for  incident  plane  wave 

In  Eq.  (5.119),  iRj  consists  of  the  zero-order  term  and  the  second-order  term; 
the  first-order  r-matrix  is  zero  mean.  In  order  to  compare  the  second-order  term 
to  that  of  the  zero-order  term,  the  second  term  on  the  right  side  of  Eq.  (5.119) 
divided  by  the  first  term  on  the  right  side  the  equation  is  plotted  in  Fig.  5.4  as  a 
function  of  the  RMS  (root  mean  square)  surface  height.  The  grazing  angle  is  20° 
in  Fig.  5.4a,  and  it  is  40°  in  Fig.  5.4b.  The  solid  curve  results  from  using  Eq. 
(5.116),  the  dashed  curve  is  due  to  Eq.  (5.117). 

5.6.4  Modification  of  coherent  field  calculation 

Taking  the  expected  value  of  both  sides  of  Eq.  (5.5)  and  using  Eq.  (5.121)  yields 

{'F2_(K)>  =  ;Sl2(K)>Fi(K).  (5.127) 

For  the  case  in  which  medium  3  parameters  are  identical  to  medium  2 
parameters, 

xj/2_(r)  =  j  (5.128) 

Taking  the  expected  value  of  the  above  equation  and  using  Eq.  (5.127)  results  in 
an  expression  for  the  coherent  field  penetrating  the  rough  interface. 

(x|/2_(r))  =  J  A;?i2(K)^i(K)e-'^>P>(*')^e**'.  (5.129) 

Likewise,  the  coherent  field  scattered  from  the  interface  is  given  by 


(5.130) 
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Note  that  ;5’i2(K)  in  Eq.  (5.129)  is  a  coherent  transmission  coefficient,  but  ;5’ii(K) 
in  Eq.  (5.130)  is  a  coherent  reflection  coefficient.  Referring  to  Eq.  (2.91),  the 
exact  expression  for  the  zero-order  field  is  easily  modified  by  substituting  ;5’i2(K) 
for  the  transmission  coefficient. 


0 


(K)-z*2p2  (K)) 

^lPl(K) 


JoiKR). 


(5.131) 


The  square  of  the  time-dependent  coherent  field  magnitude  scattered  through 
the  rough  interface  due  to  an  incident  point  source  with  baseband  source  signal 
f/(co)  is  therefore 


/2c  (**,0  = 


OO  OO  _  .  V 

J  Ja)f7(co)e“'®'  |  j^(^kR) 


(5.132) 


5.7  Summary 

Nth-order  perturbation  recursion  equations  are  derived  in  this  chapter  that 
allow  for  simple  computation  of  the  first  and  second-order  fields  for  scattering 
through  and  from  a  fluid-fluid  rough  interface.  An  expression  for  the  second- 
order  coherent  field  is  derived  that  illustrates  when  it  is  necessary  to  include  the 
second-order  field  in  the  coherent  intensity  calculations  for  scattering  through 
rough  interfaces.  In  the  ocean  experiment  simulations,  k\h  =  l,  and  including 
the  second-order  component  reduces  the  coherent  field  calculated  with  the  zero- 
order  result  by  almost  ten  percent.  Since  the  peak  incoherent  intensity  and 
coherent  intensity  typically  differ  by  over  an  order  of  magnitude  for  all  incident 
angles,  and  at  all  receivers,  the  difference  in  the  calculated  coherent  field  by 
excluding  the  second-order  correction  does  not  affect  the  intensity  velocity 
algorithm  output. 
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The  results  in  this  chapter  are  general  in  that  the  medium  below  the  lowest 
point  on  the  rough  interface  is  allowed  to  be  vertically  stratified,  viscoelastic,  or 
porous,  and  is  contained  in  the  reflection  coefficient  of  the  corresponding  fiat 
surface  problem.  For  the  case  of  increasing  sound  speed  with  depth,  these 
results  show  that  an  interference  effect  can  occur  when  enough  energy  is 
redirected  back  toward  the  interface.  Since  the  sediment  parameters  such  as 
sound  speed,  sound-speed  gradient,  and  density  are  not  everywhere  constant, 
the  effect  of  treating  the  sediment  parameters  as  random  is  discussed  in  the  next 
chapter. 


Figure  5.1  :  Definition  of  the  scattering  problem  for  an  arbitrary  two-dimensional  rough 
interface.  Medium  1  (z  >  hf{R))  and  medium  2  {-d  <z<  hf{R))  are  homogenous 
fluids:  the  medium  parameters  can  be  a  function  of  z  for  (z  <  -d).  This  diagram  can  be 
viewed  as  a  slice  through  a  two-dimensional  interface. 
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grazing  angle  (degrees) 


Figure  5.2 ;  Medium  2  is  very  fine  sand  with  k  =  1/(1. 12- 1  0.004),  and  p  =  1.85. 

Medium  3  is  sedimentary  rock  with  =  l/(2.3-/- 0.004),  Kj  =l/(1.3-t-0.11),  and 
p3/pl=2.5.  Surface:  W2  =0.04/(271:),  y  =  4.0.  The  solid  curve  uses  these 
parameters.  For  the  dashed  curve,  k,  =  1/1.13.  For  the  dotted  curve,  there  is  no 
layering,  \.e., 
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Upward  refracting  sound-speed  gradient 


grazing  angle  (degrees) 


Figure  5.3  :  Bottom  backscattering  strength  with  and  without  a  sound-speed  gradient 
for  soft,  silty,  or  clayey  bottom:  y  =  3.25,  wj  =  1.636  x  10“^  frequency  =  300  kHz, 
p  =  1.548,  5  =  1.445x10“^,  v=0.96,  q  =1545 m/s,  gradient  =  1 .3  s“^ 
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Coherent  intensity  ratio 


k^h 


Figure  5.4a  :  Ratio  of  second-order  coherent  term  to  zero-order  term, 

Re{;?o(K);?2(K)}/|j5o(K)f  with  frequency  =  20 kHz ,  p  =  2.0,  5  =  0.0163, 

y  =  3.0,  W2  =6-2xl0"^m,  v  =  1.13,  ci=  1536  m/s,  and  incident  grazing  angle, 

Gj  =20°. 


157 


CO 


k^h 


Figure  5.4b  :  Ratio  of  second-order  coherent  term  to  zero-order  term: 

(^i;i)^Re{;?o(K);?2(K)}/|^o(K)p  with  6,  =40“,  otherwise  same  parameters  as  Fig. 

5.4a  (frequency  =  20 kHz ,  p  =  2.0,  5  =  0.0163,  y=3.0,  ^2  =6.2x10  ^m,  v  =  1.13, 
cj  =  1536  m/s). 


CHAPTER  6 


First-order  perturbation  solution  for  scattering  from  a  rough  seafloor 
including  the  effects  of  random  subsurface  stratification 


6.1  Introduction  to  chapter  6 

The  example  in  Fig.  5.3  shows  that  upward  refraction  due  to  a  sound-speed 
gradient  in  the  sediment  results  in  constructive  and  destructive  interference  in  the 
field  scattered  back  into  the  water.  This  example  assumes  the  sound-speed 
gradient  and  other  seafloor  parameters  are  not  random.  Since  seafloor 
parameters  such  as  sound-speed  gradient  and  density  actually  vary  with 
position,  a  more  suitable  model  may  be  one  that  takes  into  account  the  inherent 
randomness  of  the  seafloor  parameters.  In  the  last  chapter,  a  first-order 
expression  for  the  scattering  cross  section  per  unit  area  per  unit  solid  angle 
(differential  cross  section)  is  derived  that  includes  the  complexity  of  the  medium 
below  the  lowest  point  on  the  rough  interface  (medium  3)  in  the  flat  surface 
reflection  coefficient.  In  this  chapter,  the  reflection  coefficient,  T23,  is  modeled  as 
a  random  variable.  The  expected  value  of  the  differential  cross  section  over  the 
random  sediment  parameters  is  then  expressed  analytically  in  terms  of  reflection 
coefficient  moments.  Analytical  expressions  for  the  reflection  coefficient 
moments  are  found  using  ray  theory.  These  analytical  results  are  compared  with 
modeling  the  sound-speed  gradient  as  a  random  variable  with  Gaussian 
probability  distribution  (Moe  and  Jackson,  1994b),  where  the  expected  value  of 
the  differential  cross  section  is  calculated  by  means  of  a  numerical  integration. 
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6.2  Averaged  scattering  strength  with  random  sound-speed  gradient 


The  first-order  perturbation  solution  for  the  differential  cross  section  for 
scattering  from  the  rough  interface  can  be  found  by  combining  Eqs.  (5.71)  and 
(5.63). 

x|a(l  +  nKf  ))(l  +  r(K,  ))  +  b(l  -  r(K^  ))(l  -  r(K,  )f  , 


where 


Ip  J  kf 


and 

Z;(k^,K,)  =  Pi(K^)Pi(K,)(p-1) 


(6.2) 


(6.3) 


are  given  by  Eqs.  (5.64)  and  (5.65).  The  vectors  and  K,  are  the  transverse 

components  of  the  scattered  field  and  incident  field  wave  vectors,  respectively. 
The  mean  (flat)  surface  reflection  coefficient,  T,  is  written  explicitly  in  terms  of  the 
reflection  coefficients  ri2  and  r23  in  Eq.  (5.43),  where  ri2  represents  the  mean 
surface  reflection  coefficient  at  z  =  0  as  d^oo  and  ^23  is  the  reflection 
coefficient  corresponding  to  the  interface  between  medium  2  and  medium  3 
{z  =  -d).  All  of  the  complexity  of  the  lower  medium  is  contained  in  r23.  As 
shown  in  Fig.  5.3,  the  lower  medium,  characterized  in  Eq.  (5.43)  by  Tja,  can 
have  a  notable  effect  on  the  scattered  field  strength  if  it  redirects  a  significant 
amount  of  the  downward  propagating  energy  in  the  sediment  back  up  towards 
the  rough  interface.  Two  first-order  effects  as  depicted  in  Fig.  6.1  can  result  in 
scattering  energy  away  from  the  interface,  resulting  in  destructive  or  constructive 
interference. 
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Since  sediment  parameters  are  typically  not  constant,  the  interference 
effects  in  Fig.  5.3  may  not  actually  be  observable.  By  treating  the  gradient  as  a 
random  variable,  the  expression  for  the  differential  cross  section  given  in  Eq. 
(6.1),  weighted  by  a  Gaussian  probability  distribution  for  the  sound-speed 
gradient  (linear  increase  in  squared  wave  number  with  depth)  with  variance 

(|g-(^)f)  =  0-lsec"^  and  mean  (g)  =  1.3sec"^  is  numerically  integrated.  The 

resulting  averaged  backscattering  strength  is  plotted  in  Fig.  6.2  as  a  function  of 
grazing  angle  and  compared  with  the  backscattered  field  strength  due  only  to  the 
rough  seafloor  surface  with  no  gradient,  r23  =  0.  In  this  example,  a  sound-speed 
gradient  in  the  sediment  results  in  a  substantial  increase  in  the  average 
backscattered  field  strength.  This  scattering  gain  occurs  at  all  grazing  angles, 
but  is  greater  at  smaller  grazing  angles.  Large  grazing  angles  result  in  long 
refracted  paths  in  the  lossy  sediment,  and  therefore  contribute  less  to  this  gain 
than  small  grazing  angles  with  shorter  paths. 

Using  the  same  parameters  as  in  the  above  example,  we  plot  the  bistatic 
scattering  strength  in  Fig.  6.3  as  a  function  of  scattered  grazing  angle  for  an 
incident  grazing  angle  9/  =  15°.  The  scattered  direction  in  all  plots  is  chosen 
such  that  the  incident  transverse  wave  vector,  K,-,  and  the  scattered  transverse 
wave  vector,  Kf  ,  are  in  the  same  line,  and  the  scattered  grazing  angle,  9y^,  is 
measured  from  the  horizontal  in  the  backward  direction.  In  this  case, 
backscattering  (K^  =  -Kj)  corresponds  to  an  angle  of  Qf  =  15°.  Scattering  gain 

due  to  the  upwardly  refracting  sediment  occurs  in  all  scattered  directions, 
including  the  backscattered  direction,  where  a  local  scattering  gain  maximum,  or 
backscattering  enhancement,  is  apparent.  The  reason  for  this  backscattered 
peak,  as  illustrated  in  Fig.  6.4,  is  that  the  coherent  sum  of  the  two  paths 
corresponding  to  the  incident  and  scattered  direction  occurs  only  in  the 
backscattered  direction,  because  only  in  this  direction  are  the  paths  ABCD  and 
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D'CBA'  of  equal  length.  This  effect  is  illustrated  for  similar  scattering  problems 
by  Lang  and  Khadr  (1992)  in  their  analysis  of  remote  sensing  of  a  vegetation 
canopy  over  a  flat  ground,  and  in  Hanson  and  Zavorotny  (1995)  in  their  analysis 
of  backscattering  from  an  ocean-like  surface.  The  width  of  the  peak  depends  on 

the  sound-speed  gradient  variance,  The  larger  the  sound-speed 

gradient  variance,  the  narrower  this  peak.  Another  oscillation  that  depends  on 
the  sound-speed  gradient  can  be  observed  just  below  a  scattered  grazing  angle 
of  10°.  At  these  angles,  the  oscillations  in  the  scattering  strength  are  not  quite  as 
rapid  as  those  at  larger  grazing  angles,  and  the  variance  chosen  for  the  gradient 
is  not  sufficient  to  completely  wash  out  the  dependency  of  scattering  strength  on 
scattered  grazing  angle. 

6.3  Formally  averaged  differential  cross  section  with  random  sediment 
parameters 

Rather  than  using  numerical  integration  as  in  the  last  example,  using  a 
general  approach  to  this  problem  leads  to  a  formally  averaged  analytic 
expression  for  the  differential  cross  section  per  unit  area  that  includes  the 
randomness  of  seafloor  parameters  below  /(R),  in  addition  to  the  randomness  of 
the  interface.  All  of  the  complexity  of  the  lower  medium  in  Eq.  (6.1)  is  contained 
in  the  reflection  coefficient,  r.  Consequently,  the  inherent  randomness  of  the 
sediment  is  equivalently  represented  in  the  random  reflection  coefficient,  r.  This 
point  of  view  implies  that  the  expression  in  Eq.  (6.1)  is  actually  a  conditional 
expected  value  over  the  rough  interface  given  deterministic  sediment  parameters 
below  the  interface,  or  a  deterministic  reflection  coefficient,  and  the  goal  is  to  find 
the  expected  value  of  Eq.  (6.1)  over  the  random  fluctuations  in  the  seafloor 
parameters  below  z  =  /(R), 
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WCKy.  -  K,. )(  a(l  +  Fy  ){1  +  Fi) + (.(l  -  Ty  )(1  -  r, ) 


(6.4) 


The  approach  is  to  express  Eq.  (6.4)  in  terms  of  moments  of  the  flat  surface 
reflection  coefficients.  The  notation  Tj-  =  r(K j-)  and  T,  =  r(K,)  is  used  for  the 
reflection  coefficient  in  Eq.  (6.1).  Expanding  the  square  in  Eq.  (6.4)  yields 


dii=^W(K^-K,) 


a(l  +  F^)(l  +  F,.)f)  +  fc2(|l-F^)(l-n)| 
+2fcRe(a(i + r/)(i + r,)(i  -  r})(i  -  r*)' 


(6.5) 


For  monostatic  scattering,  the  random  variables  and  F,  are  equal 
(Ky  =  -K/),  and  Eq.  (6.5)  becomes 


du=^W(K/-K;) 


x{lap[l  +  4(|Fp)  +  (iFp)  +  2Re(F2)  +  4Re(F)  +  4Re(F|Fp)] 
+b^\l  +  4(|Fp)  +  (|Fp)  +  2  Re(F2)  -  4Re(F)  -  4Re(F|Fp)] 
+bRt(a)[2  -  8(|Fp)  +  2(|Fp)  +  4Re(F2)] 
+8i>Im(a)[lm(F|Fp)  -  Im(F)]|. 


(6.6) 


Small  changes  in  medium  parameters  result  in  rapid  changes  in  phase  of  the 
reflection  coefficient.  The  reflection  coefficients,  Tf  and  F,,  are  therefore 
assumed  to  be  independent  random  variables  for  bistatic  scattering.  However, 
for  scattering  in  directions  close  to  the  backscattered  direction,  the  reflection 
coefficients  corresponding  to  the  incident  and  scattered  directions  are  neither 
equal  or  independent  (see  Fig.  6.3).  For  independent  reflection  coefficients,  Eq. 
(6.5)  becomes 
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dn=^W(K/-K,-) 

xjlap^i + 2  Re(r^) + ^|r^f  )][i + 2  Re(r,> + (|r,-p)_ 

+b^\^i  -2Rc{Tf) + (|r/f  )][i  -  2  Re(r,.) + (|r,.p)] 

+2bRc{a)^[l  -  -  (|r,.p))  -  4Im(r^)lm(r,-) 

In  order  to  find  expressions  for  the  moments  (F),  (r^).  (in^)>  (r|rp^,  and  (|rp). 

an  equivalent  form  of  Eq.  (5.43)  is  used  to  express  the  reflection  coefficient,  F,  in 
terms  of  the  partial  reflection  coefficients  F12,  and  F23  as  a  geometric  sum, 

r=r,2+(i-rf2)r23e^‘'‘*''  i(-ri2r23e'**‘')".  (6.8) 

n=0 

Small  variations  in  medium  parameters  result  in  small  deviations  of  the 
magnitude  of  the  reflection  coefficient  F23  from  its  mean,  but  result  in  large 
fluctuations  of  its  phase.  The  partial  reflection  coefficient  F23  is  assumed  to  have 
the  form 

F23(K)e^‘‘'2p24  ^  (6.9) 

where  (|)  is  a  uniformly  distributed  random  variable  between  0  and  2jt,  and 
ti^(K)  is  approximated  as  the  mean  square  reflection  coefficient  ^|F23(K)py 
With  much  work,  Eqs.  (6.8)  and  (6.9)  yield 

(F)  =  Fi2,  (6.10) 

r^)=r?2 


(6.11) 
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(6.12) 


(6.13) 


=  1^121  +^1^12]  l-r'12 


i-|ri2rTi^ 


+i-ri2  Ti 


2  "^.^4  1  +  1^21 


-4Re  rf2  i-ri2  i-rf2 


.2  |2  Tj^ 


i-|ri2ni' 


(6.14) 


The  derivations  for  Eqs.  (6.12),  (6.13),  and  (6.14)  are  given  in  Appendix  E;  Eqs. 
(6.10)  and  (6.11)  follow  immediately  from  Eqs.  (6.8)  and  (6.9). 


6.4  Effect  of  sound-speed  gradient 

The  above  formalism  is  now  applied  to  the  case  of  a  linear  decrease  in 
squared  wave  number  with  depth  introduced  earlier  (Figs.  5.3,  6.2,  and  6.3).  In 
order  to  use  the  differential  cross  section  expressions  in  Eqs.  (6.6)  and  (6.7),  the 

value  for  Ti^(K)  =  ^|r23(K)f^  {d2=0)  is  computed  and  used  in  Eqs.  (6.10) 

through  (6.14).  These  equations  are  then  substituted  in  Eq.  (6.6)  to  obtain  the 
monostatic  differential  cross  section,  or  in  Eq.  (6.7)  to  obtain  the  bistatic 
differential  cross  section  for  scattered  directions  not  close  to  the  backscattered 
and  specular  directions.  Rather  than  use  Eq.  (5.81)  to  determine  ti^(K)  exactly, 
a  computationally  efficient,  yet  accurate,  ray  theory  expression  (Appendix  F)  is 

used  to  obtain  the  following  analytical  expression  for  ^|r23(K)p^: 
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|r23(Kf)  =  Tl^(K)  =  exp 


— sin^(02)  +  2cos^(92)sin(92) 
V  ^  L3 


(6.15) 


where  02  is  the  grazing  angle  just  below  the  mean  surface, 


cos(02)  = 


K 

Re(A:2)' 


(6.16) 


The  gain  parameter,  q,  is  the  sound-speed  gradient,  g  =  dcldz  evaluated  at 
2  =  0,  divided  by  the  product  of  the  frequency,  /,  and  the  loss  parameter,  6; 


g  _l  dX 
fb  b  dz 


(6.17) 


This  parameter  will  be  shown  to  determine  the  average  gain  of  the  scattered  field 
strength  due  to  gradients  and  discontinuities  below  the  rough  water-sediment 
interface  as  compared  to  the  scattered  field  strength  due  solely  to  the  rough 
water-sediment  interface.  For  example,  a  small  value  of  q  due  to  a  relatively 
large  loss  at  a  high  frequency  will  result  in  the  attenuation  of  much  of  the  signal 
energy  before  it  strikes  the  rough  interface  again  from  below,  resulting  in  little 
scattering  gain  over  the  homogenous  fluid-fluid  case.  Also,  a  small  sound-speed 
gradient  results  in  large  path  lengths,  and  attenuation  of  much  of  the  signal 
energy  before  it  strikes  the  rough  interface  for  the  second  time.  When  q  is  small, 
only  the  sediment  properties  at  the  water-sediment  interface  are  of  importance  in 
the  scattered  field  strength  calculation  . 

An  expression  for  ti^(K)  for  the  linear  increase  in  sound  speed  with  depth 
case  is  also  derived  in  Appendix  F. 


166 


6.4.1  Effect  of  sound-speed  gradient  on  backscattered  field  strength 

As  a  check  of  the  analytical  results,  the  backscattering  strength  using  Eq. 
(6.6)  is  plotted  in  Fig.  6.5  along  with  the  numerically  obtained  result  from  Fig.  6.2, 
and  the  bistatic  scattering  strength  obtained  from  Eq.  (6.7)  is  plotted  in  Fig.  6.6. 
It  is  evident  that  the  scattered  field  strength  obtained  by  numerically  averaging 
the  expression  in  Eq.  (6.1)  over  a  Gaussian  probability  distribution  function  for 
the  gradient  is  indistinguishable  from  the  result  due  to  the  formally  averaged 
differential  cross  section  in  the  backscattered  direction  and  for  bistatic  scattering 
at  least  a  few  degrees  from  the  backscattered  direction.  In  this  example,  the 
gradient  gain  parameter  has  the  relatively  large  value  ^  =  30.  Equation  (6.6)  is 
used  to  find  the  point  on  the  plot  corresponding  to  backscattering. 

When  =-Kj,  then  r(Ky)  =  r(Kj),  resulting  in  a  peak  in  the 

backscattered  direction.  Depending  on  the  value  of  the  sound-speed  gradient 
variance  in  the  numerical  calculation,  the  reflection  coefficients  r(Ky)  and  r(Kj) 

are  correlated  when  Kf  is  close  to  Ki,  and  therefore  the  bistatic  scattered  field 

strength  is  dependent  on  the  gradient  variance  in  the  region  close  to  the 
backscattered  direction.  However,  even  for  a  small  variances  in  seafloor 
parameters,  this  region  is  very  small. 

The  moments  obtained  in  Eqs.  (6.10)  through  (6.14)  using  the  representation 
for  the  random  reflection  coefficient  in  Eq.  (6.9)  also  agree  very  closely  with  the 
corresponding  reflection  coefficient  moments  obtained  numerically  with  a 
Gaussian  probability  distribution  function  for  the  sound-speed  gradient. 

Figure  6.7a  illustrates  the  dependency  of  the  scattering  gain, 
101ogio(an/aii),  (K/  =  -K/),  where  On  is  given  by  Eq.  (5.71),  on  the  mean 
density  ratio  and  mean  sound-speed  ratio  for  q  =  30.  This  value  for  q  is  the 
same  as  that  of  Figs.  6.2,  6.3,  6.5,  and  6.6.  Fig.  6.7b  is  an  identical  plot,  but  with 
a  lower  gain  parameter,  ^  =  10,  resulting  in  correspondingly  lower  values  for  the 
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scattering  gain.  For  a  density  ratio  close  unity,  the  scattering  gain  can  be  very 
large.  However,  the  density  ratio  must  always  be  greater  than  one  for  the 
seafloor. 

The  dependence  of  the  scattering  gain  on  q  and  grazing  angle  is  shown  in 
Fig.  6.7c.  As  can  be  seen  in  Fig.  6.3  and  Fig.  6.7c,  the  scattering  gain  is  strongly 
dependent  on  the  grazing  angle.  A  low  grazing  angle  results  in  a  significantly 
higher  scattering  gain  because  the  resulting  ray  path  length  in  the  lower  sediment 
is  much  shorter  and  therefore  attenuated  much  less  than  the  large  grazing  angle 
case.  The  dependency  of  the  scattering  gain  on  q  is  greater  for  larger  grazing 

angles. 

6.4.2  Bistatic  scattering  and  backscattering  enhancement 

Figures  6.8a,  6.8b,  6.8c,  and  6.8d  use  Eqs.  (6.6)  and  (6.7)  to  show 
backscattering  enhancement  corresponding  to  the  same  sediment  parameters  in 
Figs.  6.7a,  6.7b,  6.7c,  and  6.7d.  Note  that  these  results  do  not  depend  on 
interface  roughness.  Since  the  two  first-order  paths  through  the  sediment  add 
coherently  only  for  monostatic  scattering  (see  Fig.  6.4),  the  field  strength 
scattered  back  in  the  direction  of  the  transmitter  is  greater  than  the  field  strength 
scattered  in  other  nearby  directions.  As  is  clear  from  Fig  6.8c,  the  dependency  of 
backscattering  enhancement  on  q  is  negligible. 


6.5  Conclusion 

A  formally  averaged  expression  for  the  scattered  field  that  includes  the 
random  stratification  below  the  rough  seafloor  surface  as  well  as  the  randomness 
of  the  rough  seafloor  is  derived  in  this  chapter.  Although  this  formalism  can  be 
applied  to  a  rough  interface  above  any  arbitrary  vertically  stratified  sediment,  ray 
theory  is  put  to  use  to  obtain  formally  averaged  closed  form  analytical 
expressions  for  the  special  case  of  linear  decrease  of  squared  wave  number  with 
depth  sound-speed  gradient.  This  result  shows  that  the  scattering  gain  due  to 
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upward  refracting  sound-speed  gradients  is  determined  by  the  quantity 
q  =  gj[f^),  the  sediment  parameters  (speed  ratio,  and  density  ratio),  and  also  the 

incident  angle. 

Gradients  and  layering  below  the  seafloor  can  strongly  affect  the  scattering 
from  the  interface  if  a  significant  amount  of  the  energy  that  penetrates  the 
interface  is  redirected  back  up  towards  the  interface.  The  effect  of  gradients  in 
sound-speed  on  scattering  can  be  significant  even  for  incident  energy  at  large 
grazing  angles  (for  situations  with  large  q  —  large  gradient,  low  frequency,  low 
loss),  but  the  effect  of  gradients  on  the  scattered  field  strength  is  most  significant 
at  small  grazing  angles.  Backscattering  enhancement  is  also  shown  to  occur  in 
sediments  that  exhibit  a  relatively  large  scattering  gain.  By  using  numerical 
solutions  for  the  flat  surface  reflection  coefficient,  the  monostatic  and  bistatic 
differential  cross  sections  of  a  rough  interface  above  an  arbitrarily  complicated 
vertically  stratified  sediment  that  includes  the  randomness  of  the  sediment 
parameters  below  the  interface,  in  addition  to  the  randomness  of  the  interface, 
can  be  solved  using  the  formalism  in  this  chapter. 

The  results  of  this  chapter  can  also  be  used  to  obtain  expressions  for  the 
field  scattered  through  the  rough  interface  when  the  sound  speed  is  depth 
dependent.  These  results  could  then  be  incorporated  in  the  model  used  in  the 
ocean  experiment  simulation  discussed  in  chapter  4.  Except  for  downward 
propagating  energy,  sound-speed  gradients  would  increase  the  path  length  of 
energy  scattering  into  the  sediment.  Consequently,  sound-speed  gradients 
would  seem  to  favor  downward  propagating  energy,  resulting  in  a  slower 
apparent  slow  wave  propagation  than  that  predicted  by  the  model  used  in 
chapter  4.  The  effect  of  sound-speed  gradients  on  the  scattered  field  through  the 
interface  is,  however,  beyond  the  scope  of  this  dissertation,  and  a  subject  for 
future  investigation. 
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Constructive  or  destructive  interference 


Figure  6.1  :  Two  first-order  effects  that  can  result  in  destructive  or  constructive 


interference. 
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Average  effect  of  gradients 


grazing  angle  (degrees) 


Figure  6.2  :  Averaged  backscattering  strength.  Linear  decrease  in  squared  wave 
number  with  depth  sound-speed  gradient,  g,  with  Gaussian  probability  function 

((g)  =  1.3  s'^  /|g|^  -  {g)\  =  0.1  s'^).  All  other  parameters  same  as  Fig.  5.3. 


bistatic  scattering  strength  (dB) 
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Figure  6.3  :  Numerically  averaged  bistatic  scattering  strength;  same  parameters  as  Fig. 
6.2  and  Fig.  5.3. 
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Backscattering  enhancement 


Figure  6.4  :  Backscattering  enhancement.  If  paths  AJBCD  and  D'CBA'  are  equal 
(A  =  A'  and  D  =  D'),  then  backscattering  enhancement  due  to  the  sound-speed 
gradient  occurs. 


backscattering  strength  (dB) 
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Comparison  of  analytical  and  numerical  results 


Figure  6.5  :  Average  backscattering  strength:  comparison  of  numerical  average 
obtained  with  random  sound-speed  gradient  (Fig.  6.2)  with  analytical  formal  average 
approximation.  Same  parameters  as  Fig.  6.3. 
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Comparison  of  analytical  and  numerical  results 


scattered  grazing  angle  (degrees) 


Figure  6.6  :  Average  bistatic  scattering  strength:  comparison  of  numerical  average 
obtained  with  random  sound-speed  gradient  (Fig.  6.3)  v>/ith  analytical  formal  average 
approximation.  Same  parameters  as  Fig.  6.3. 
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Figure  6.7a  :  Dependency  of  scattering  gain,  101ogio(a/a),(Kj  =-Ky),  on  the  mean 

density  ratio  and  mean  sound  speed  ratio  for  q  =  30  (frequency  =  1  kHz,  g  =  30  s  ^ 
6  =  0.001,  incident  angle  =  45°) 
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Backscattering  Enhancement  (dB) 


q 


Figure  6.8c  :  Backscattering  enhancement  as  a  function  of  grazing  angle  and  q 
p  =  1.5,  f  =  lkHz,  5  =  0.001,  v  =  1.0.  q  =1500m/s. 


CHAPTER  7 


Dissertation  Summary  and  Proposed  Future  Work 


7.1  Introduction  to  chapter  7 

By  means  of  analytical  models  for  seafloor  scattering  derived  in  this 
dissertation,  roughness  is  shown  to  be  an  important  mechanism  for  penetration 
of  acoustical  energy  into  the  sediment  at  subcritical  grazing  angles. 
Experimental  results  (Boyle  and  Chotiros,  1992;  Chotiros,  1995)  showing 
acoustic  penetration  from  water  into  sandy  sediments  at  grazing  angles  below 
the  compressional  critical  angle  in  relation  to  the  mean  surface  are  explained  in 
terms  of  diffraction  of  an  ordinary  longitudinal  wave,  rather  than  refraction  of  a 
slow  wave. 

Since  the  transmitted  signal  in  the  experiments  is  a  pulse,  formalism  for 
scattering  of  a  pulse  through  a  rough  interface  as  well  as  from  a  rough  interface 
is  derived  in  terms  of  the  scattering  cross  section  per  unit  area  per  unit  solid 
angle  (differential  cross  section).  This  analytical  result  for  the  time-dependent 
scattered  field  intensity  is  used  to  show  that  the  CW  approximation  is  appropriate 
for  simulations  of  both  experiments.  First-order  perturbation  theory  is  used  to 
obtain  an  approximate  expression  for  the  differential  cross  section.  Since  using 
the  zero-order  field  to  represent  the  coherent  intensity  is  only  valid  for  small 
amounts  of  interface  roughness,  the  second-order  component  of  the  coherent 
field  intensity  is  derived,  and  compared  with  the  zero-order  component  of  the 
coherent  intensity  for  the  sediment  roughness  parameters  used  in  the  ocean 
experiment  simulation.  These  analytical  results  justify  using  the  flat  surface  field 
to  represent  the  coherent  component  of  the  field.  Although  the  main  problem  is 
acoustic  scattering  through  the  water-sediment  interface  for  an  incident  field  in 
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the  water,  analytical  results  derived  in  this  work  also  include  scattering  from  the 
interface.  Furthermore,  these  analytical  results  are  general,  and  can  be  applied 
to  other  rough  interface  scattering  problems. 

Gradients  in  sediment  sound  speed  are  also  allowed  in  the  perturbation 
results.  A  novel  approach  for  finding  the  differential  cross  section,  given 
vertically  stratified  random  sediment  parameters,  yields  analytical  results  that 
show  the  average  effect  of  sound  speed  gradients  in  the  sediment  on  the  field 
scattered  back  into  the  water.  Although  gradients  in  sound  speed  may  affect  the 
apparent  slow  wave  speed,  this  effect  is  not  calculated  in  this  dissertation. 
However,  the  formalism  in  this  dissertation  can  be  generalized  to  include  the 
effect  of  sound-speed  gradients  on  the  field  penetrating  the  surface. 

7.2  Time  dependent  scattered  intensity  —  chapter  2 

An  original  expression  for  the  two-frequency  mutual  coherence  function  for 
scattering  from  and  through  a  2-D  fluid-fluid  rough  interface  due  to  a  narrowband 
incident  plane  wave  is  derived  and  expressed  in  terms  of  the  second  moment  of 
the  r-matrix,  resulting  in  a  general  analytical  expression  for  the  time  dependent 
incoherent  intensity  due  to  a  narrowband  incident  plane  wave. 

Two  approximations  are  made  to  the  two-frequency  mutual  coherence 
function  when  obtaining  the  results  of  chapter  2. 

•  The  r-matrix  correlation  function  is  assumed  to  vary  slowly  over  the  source 
frequency  range.  This  assumption  is  shown  to  be  valid  in  chapter  3  for  a 
narrowband  source  signal. 

•  The  phase  terms  are  expanded  in  a  power  series  in  frequency. 

An  equivalent  expression,  omitting  evanescent  waves,  is  obtained  in  terms  of  the 
differential  cross  section.  It  has  been  demonstrated  that  the  evanescent 
component  of  the  scattered  incoherent  field  is  insignificant  compared  to  the 
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propagating  component  of  the  incoherent  field  for  distances  as  small  as  one  or 
two  wavelengths  from  the  interface,  depending  on  the  surface  roughness 
spectrum.  In  the  limit  as  the  narrowband  signal  approaches  the  CW  (continuous 
wave)  case,  the  result  in  terms  of  the  T-matrix  correlation  function  is  exact,  and 
the  result  in  terms  of  the  differential  cross  section  approaches  the  exact  solution 
for  distances  sufficiently  far  from  the  interface.  This  result  is  interesting,  because 
the  differential  cross  section  is  defined  as  a  far-field  entity. 

This  expression  for  the  incoherent  time  dependent  scattered  field  intensity  in 
terms  of  the  differential  cross  section  is  used  to  heuristically  obtain  the 
corresponding  expression  for  the  scattered  field  due  to  a  narrowband  point 
source.  A  simple  analytical  expression  is  obtained  for  the  case  of  a  Gaussian 
shaped  source  pulse,  also  in  terms  of  the  differential  cross  section.  The 
analytical  results  and  calculations  in  this  chapter  show  that  the  sonar  equation  is 
appropriate  for  the  simulations  of  the  experiments. 

7.3  First-order  perturbation  theory  —  chapter  3 

Derivations  for  the  first-order  perturbation  expressions  used  for  the  numerical 
computations  in  chapter  2  and  in  the  computer  simulations  (chapter  4)  for 
scattering  both  from  and  through  a  rough  fluid-fluid  interface  separating  two 
homogenous  fluids  are  presented  in  chapter  3.  These  first-order  perturbation 
theory  results  are  also  used  to  show  that  the  T-matrix  correlation  function  varies 
sufficiently  slowly  over  the  source  frequency  range  so  that  its  dependence  on 
frequency  can  be  ignored  in  the  time  dependent  scattered  field  derivation. 

7.4  Computer  simulations  of  recent  experiments  —  chapter  4 

The  results  of  computer  simulations  of  recent  experiments,  presented  in 
chapter  4,  show  the  effect  of  roughness  on  acoustic  penetration  of  the  seafloor 
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and  are  compared  with  recent  sediment  penetration  experiments.  These  results 
show  that  the  incoherent  intensity  scattered  through  the  rough  seafloor  can 
mimic  a  slow  refracted  wave. 

By  assuming  a  small  amount  of  roughness  at  the  water-sediment  interface, 
and  using  the  3-D  scattering  model  developed  in  chapter  2  and  the  first-order 
perturbation  theory  results  for  2-D  fluid-fluid  interfaces  derived  in  chapter  3,  it  is 
possible  to  reproduce  the  results  of  the  tank  experiments  of  Boyle  and  Chotiros 
(1992),  both  in  magnitude  and  arrival  time  of  the  incoherent  intensity  pulse 
(apparent  slow  wave).  These  results  indicate  that  the  acoustic  penetration  of  the 
surface  may  be  due  to  scattering  (diffraction)  from  low  levels  of  roughness  rather 
than  slow-wave  refraction. 

Simulations  of  the  ocean  sediment  penetration  experiments  (Chotiros,  1995) 
are  presented  and  compared  with  Chotiros'  experimental  results.  Assuming 
roughness  parameters  appropriate  to  the  experiment  location  and  using  the 
intensity  propagation  speed  and  direction  finding  algorithm  used  by  Chotiros, 
these  simulations  also  show  propagation  of  energy  at  speeds  slower  than  the 
compressional  speed  of  the  sediment,  in  a  direction  concurring  with  Snell  s  law. 
Although  the  simulations  yield  an  apparent  slow  wave  that  is  between  50m/s  and 
300  m/s  faster  than  the  results  presented  by  Chotiros,  the  apparent  slow  wave 
speed  is  dependent  on  the  roughness  spectrum.  Because  of  the  potential  for 
disturbance  of  the  sediment  above  the  receivers,  the  roughness  spectrum 
measured  in  the  area  of  the  experiment  is  not  the  same  as  that  of  the  actual 
roughness  spectrum  at  the  experiment  site. 

7.5  Perturbation  theory  results  for  a  fluid-fluid  interface  including 
gradients  —  chapter  5 

Using  Rayleigh-Rice  perturbation  theory.  Nth  order  recursion  equations  are 
derived  in  chapter  5  that  provide  a  starting  point  for  evaluating  the  Nth  order 
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scattered  field  in  terms  of  lower  orders.  The  sediment  below  the  lowest  point  on 
the  interface  can  be  viscoelastic,  or  porous,  supporting  shear  or  Biot  slow  waves 
as  well  as  gradients  in  physical  properties  (e.g.,  sound  speed  can  vary  with 
depth),  and  is  represented  by  the  corresponding  flat  surface  reflection  coefficient. 
For  the  case  of  increasing  sound  speed  with  depth,  these  results  (first-order) 
show  that  constructive  and  destructive  interference  in  the  field  strength  can  occur 
for  scattering  back  into  the  water  when  a  significant  amount  of  energy  is 
redirected  back  toward  the  surface. 

An  expression  for  the  second-order  coherent  field  is  derived  from  the  Nth 
order  recursion  equations  that  illustrates  the  necessity  of  including  the  second- 
order  component  in  the  coherent  intensity  calculations  for  k\h>\. 

7.6  Formal  average  of  sediment  parameters  —  chapter  6 

Since  the  sediment  parameters  such  as  sound  speed,  sound  speed  gradient, 
and  density  generally  show  variability  in  lateral  direction,  the  effect  of  treating  the 
sediment  parameters  as  random  is  shown  in  chapter  6.  The  first-order 
expression  for  the  differential  cross  section  derived  in  chapters  2  and  5  is 
generalized  to  account  for  randomness  in  the  medium  below  the  lowest  point  on 
the  rough  interface  by  expressing  the  differential  cross  section  in  terms  of 
moments  of  the  stochastic  flat  surface  reflection  coefficient.  This  approach 
results  in  a  formally  averaged  expression  for  the  scattered  field  that  includes  the 
random  stratification  in  the  sediment,  in  addition  to  the  randomness  of  the  rough 
water-sediment  interface. 

Although  this  formalism  can  be  applied  to  a  rough  surface  above  any 
arbitrary  vertically  stratified  sediment,  ray  theory  is  used  to  obtain  formally 
averaged  closed  form  analytical  expressions  for  a  few  special  cases  that  include 
sound  speed  gradients.  This  result  shows  that  the  scattering  gain  due  to  upward 
refracting  sound  speed  gradients  is  determined  by  the  quantity  q  =  gl{0),  the 
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sediment  parameters  (speed  ratio  and  density  ratio),  and  also  the  incident  angle. 
Upward  refraction  or  reflection  due  to  large  sound  speed  gradients  in  the 
sediment  results  in  increased  scattering  strength  over  the  non-gradient  case. 
Backscattering  enhancement  is  demonstrated  in  this  first-order  result,  since  the 
scattering  gain  has  a  relative  maximum  when  the  incident  and  scattered  grazing 
angles  are  equal. 

Summary  and  Ideas  for  future  work 

The  tank  and  ocean  experiment  simulation  results  presented  in  this 
experiment  show  that  surface  roughness  is  a  mechanism  for  subcritical 
penetration  of  sediment.  Although  the  issue  of  the  dominant  mechanism  for 
subcritical  penetration  of  the  seafloor  can  only  be  resolved  by  experiment,  further 
improvements  on  the  scattering  model  can  be  made.  For  example,  the  effect  of 
sound  speed  gradients  on  the  scattered  field  through  the  interface  is  a  subject  for 
future  investigation.  The  results  of  chapters  5  and  6  can  be  used  to  obtain 
analytical  expressions  for  the  field  scattered  through  the  rough  interface  when 
the  sound  speed  is  depth  dependent.  These  results  could  then  be  incorporated 
in  the  model  used  in  the  ocean  experiment  simulation. 

The  Nth  order  recursion  equations  provide  a  starting  point  to  solve  higher 
orders.  Solving  for  the  third-order  scattered  field  would  provide  a  means  of 

obtaining  the  differential  cross  section  to  This  result  could  then  be  used  to 

determine  the  accuracy  of  the  first-order  results  for  scattering  through  and  from 
2-D  surfaces. 

The  results  in  chapter  5  provide  a  straightfonvard  approach  for  obtaining  the 
first-order  scattered  field  from  multiple  rough  interfaces  in  terms  of  partial 
reflection  coefficients.  The  effect  of  multiple  rough  interfaces  could  then  be 
compared  with  the  results  of  chapter  6. 
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APPENDIX  A 


Derivation  of  narrowband  dispersion  coefficients 

Equations  (2.38)  -  (2.41)  (expressions  for  rj,  t2,  P,  and  Q)  are  derived  in  this 
section.  For  convenience,  define 
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Equation  (2.42)  is  included  below  for  convenience: 
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i2 


2B2c 


2k^  2  f  cos  9; 

2 


JX  J 

C\  C\ 


rB2 


9co 


//2 


l(0j=(0"=(0^ 


4^ 


1  lOi-K  „ 
— - COS0,- 


n2 


1  f  COS0A^ 


q^^A  2 


and 


d^B2 


9c0;a), 


co;=co,=cOc 

=  _-5^iJ^cos0j 
4q5|c 


^2k^G)c  PLj-K 


V 


q 


q 


COS0,- 


q^^c 


'^C0S9;^^ 


(A13) 

(A14) 


(A15) 

(A16) 


(A17) 


Substituting  (A12)  -  (A17)  into  (A4)  yields  an  approximate  expression  for  B2. 
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B2(i^,C0c  +o)',a)^  +0)")  = 

^^-^cose, 


+0) 


,  1  flK^Gi,  -•  ^ 


+(0' 


2B2c 

d,K 


COS0j 


-la)'2 

2 


-co'0"' 


^2k^o)^  d,'K 


4^ 


^2 


a,-  •  K  . 
— ‘-^COSS; 

4qS2c 


cf  q 

(^Jl 


cos9,- 


1 


J 


252c 

^ 


''2k^  2 


V 


2k  co^  a;  *  K  . 
- ^ - i - cos  0; 

I.  cl  / 


q  q 

1  f  COS0J 


q^52c 


--co"^ 

2 


1 


44 


fdrK  1  rcos0,f 

V,  q  J  q52c^  2  y 


Similarly, 

S2(^,Wc  +m"5®c  +“')  = 


+0) 


1  (  2k^G)c  d;  -K 


+© 


252c 
.  «rK 

2q^ 


I 

cos0j 


COS0J 


2 


1 


452^ 


/■  9  ^  \2 

2k  CO.  a.-  K  -  . 

- - - - COS0/  cos(p 

q 


2^c 


-©'(0"' 


-i-(0'2 

2 


^COS0, 

4q52\ 


^ 


1 


45ic 


f  0 

2k  (0^  a.  -K  _ 

— ^  — * —  cos  0,-  cos  (]) 

q  ^1  y 

^2 

I  /  f  iw  n 

i 


^2k^  2  f  COS0,- 

1  ^  COS0,- 

^  2  y 


2^ 


q^52c 


1 


COS0;A^ 


y 


q^52c 


^  2 


(AI8) 


(A19) 


Subtracting  Eq.  (A18)  from  the  conjugate  of  (A19)  yields  an  approximation  for 
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:^{K,(Sic  +  co",cOc  +  co')  -  B2(K,(Sic  +  co'.co^  +  co")  = 


*2 

K  CO. 


d,-K  „  a,K  „ 

- ^  '  cosO. 

-  2ciB2c 


cMc  ^C\Bic 


+CO 


.  '  «rK 
l2ci52*. 


COS0/ 


k'^cOc  di  -K 


c\B2c  2qB2c 


cosOj 


,  1  ..//2j  1 

H — CO  < - ^ 

2  UbI, 


arK 

I 


cosOj 


^  1  /^cosO,-^^ 


C\B2c 


n  * 

44 


:t=2 

2k  co 


£-^cose,. 


a2 


1 


24 


2k*^  2  f  COS0,- 


+ic0'2 

2 


'^2k^co 


44 


I 


L-^il^cosQi 


^2 


q 


2^c 


c  <^1  ^1  ^ 

^2k^  2  f  COS0; 


2  y 
2\ 
y 


Q  * 

44 


^cos0, 


1  f COS0, y 


2 


+co'co"2/Im 


^oose; 

^ClHc 


^2kV  6i,K 


V 


COS0, 


1  1 

^ COS0J 

^  2  J  ^ 

(A20) 


Combining  the  co'  terms  in  (A20), 


^cos0;-  Re(l/£’^ )  -  ^  ^2 


K 

^  r>^ 

^  1  ^ 

K  _1 

COS\7j  iv6 

®c 

^Kp2(is:,o)c)j 

q32(^,®c)^ 

HCO  Z, 


(A21) 


and  making  the  definitions 


t2  = 


P2(^-®c) 


(A22) 


R  Q  .  •  D  .  «rKcos0, 

= - ^-cos0, +— sin0, +z— ^ - ; - ^Re 

q  q  ^1^1 


k|52(K,Mc)J  PziK,®.) 


(A23) 
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Equation  (2.47)  is  obtained  by  rearranging  terms  in  (A24),  and  Eq.  (2.48)  is 
obtained  from  the  cross  term,  co"©'. 
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APPENDIX  B 


Dispersion  of  a  Gaussian  pulse 


_t2  L2 

Beginning  with  the  baseband  source  term  in  Eq.  (2.73),  u{t)  =  e  '  " ,  and  the 
definition  of  g^{t)  in  Eq.  (2.42),  Eq.  (2.68)  is  derived.  Using  the  well  known 


identity 


oo 


c^/2a^gi(ax  _g-<r^®V2 


in  the  second  integral  of  Eq.  (2.49)  yields 


where 


<1 


q  =  -\-iP . 

Using  (B2),  the  expression  for  g^{t)  becomes 


^2/^)  ^  f  j(o"[/*(co")e 


where 


P"  =  P*  -i{2Plq'f. 

Using  the  identify  in  (B1)  to  evaluate  (B5)  yields  Eq.  (2.75). 
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APPENDIX  C 

RMS  height :  Subtracted  Gaussian  weighted  moving  average 

The  RMS  height  of  the  surface  is  the  surface  correlation  function  evaluated  at 
R  =  0. 


h^  =  \d^KW{K). 

Substituting  Eq.  (3.66)  for  the  roughness  spectrum  into  (C1)  yields 


(C1) 


=lnw\dKK\-Q 

0 

oo 

=  lTiw\dKK^~"‘ 

0 


-¥}a^  !  2 


(C2) 


Evaluate 


5(e)  =  lim  ? dKK^~A  1  - '^  +  t 

e-40''  L 

e 

Let 

oo  2— Y 

/i  = 

J  Y-2 


=  \dKK^ 


-Y^-A’V/2 
'e  , 


and 


=  \dKK^ 


-Y, 


2 

—K  a 


e 


(C3) 


(C4) 

(C5) 


(C6) 
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The  integral  I2  is  solved  first.  Using  integration  by  parts  (let  M  =  e  ^  ^  1'^)^  this 


integral  simplifies  to 


-ih 


Y-2-'e 


Next,  73  from  (C6)  is  evaluated  using  the  solution  for  I2 .  Let  =2a^\  then 


fp. 


u  =  - i—^rdKK 

^  y_2  y -2^^ 


3-y-Kh^  12 

&  I 


^  y-2  Y-2-'e 


3-y-K^a^ 

c 


Therefore 


S(e)  =  h~  272  +  h 


o2-Y  -t^a^l2.2-y  „-e^a^o2-Y 

e ^  I  ^  ^ 

y-2  y-2 


y-2  Y“2 

+2-2Lr,iA:A;3-V*''»''^ 

Y-2’'e 

Y-2-'e 


(C10) 


The  first  three  terms  in  (C10)  sum  to  zero  in  the  limit;  the  last  two  terms  can  be 
evaluated  using  the  identity 


fOO 

f  dKK 

JO 


3-y-K^a^  1 


-r(2-Y/2), 


(C11) 


or  equivalently, 


#oo 

I  dKK 

Jo 


3-y-K^a^  12  _  2 


I-Y/2 


-r(2-Y/2). 


(C12) 


Using  Eqs.  (C11)  and  (C12)  in  (C10)  and  (C2)  yields 


^2  ^  _  y/2)[2^-V/2  _  ij 

The  RMS  surface  height  is  therefore 
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(C13) 


(C14) 
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APPENDIX  D 

Number  and  placement  of  array  elements 

Using  the  intensity  front  velocity  algorithm  (Chotiros,  1995),  determine  the 
number  of  receivers  that  are  required  to  obtain  a  point  solution  in  speed- 
depression  angle  space. 

Define  the  reference  sensor  tq  to  be  at  the  origin 

ro=0.  (D1) 

The  signal  arrives  at  the  yth  receiver  located  at  rj  =  (xj,yj,Zj)  at  time  tj. 

(D2) 

c'  ’ 

where  a  =  +  a^y  +  a^-z  is  the  unit  magnitude  direction  vector  of  the  plane 

wave,  and  c  is  the  propagation  speed.  Equation  (D2)  for  receivers  ro,ri,...,r„ 
can  be  written  as 

’-^1  zi~\  pi 

X2  y2  Z2  fajc]  h 

.  .  .  ay  =c'  .  .  (D3) 

•  •  •  l^zi 

Xn  yn  ^n\  W. 

There  are  n  + 1  receivers  including  the  reference  receiver  tq.  Note  that  all  of  the 
receiver  locations  ri,r2,...r„  and  the  receiver  times  ri,r2.  are  assumed  to  be 
known.  The  vector  a  and  the  propagation  speed  c'  are  unknown  —  three 
unknowns.  In  order  to  have  a  unique  solution,  the  receiver  position  matrix  r  is 
required  to  have  rank  3,  requiring  a  minimum  of  three  receivers  (ri,r2,r3)  that 
span  3-dimensional  space.  Including  the  reference  receiver,  a  minimum  of  four 
receivers  are  required  that  are  not  in  the  same  plane. 
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However,  Chotiros’  speed-depression  angle  solution  space  need  only  be  a 
slice  of  the  3-D  solution  space  since  the  plane  wave  direction  has  no  y 
component  (a^  =0).  Therefore,  the  receiver  matrix  r  is  required  only  to  have 

rank  2,  spanning  only  the  x,z  plane.  Including  the  reference  receiver,  a  minimum 
of  three  receivers  (w  =  2)  are  required. 

The  following  example  shows  that  a  linear  array  results  on  a  curve  in  the 
speed-depression  angle  solution  space.  Assume  a  four  element  linear  array, 
such  that 


ri  =(0,0,zi) 

*•2  =  (0.0.Z2) 

r3  =  (0,0,Z3) 


The  reference  receiver  is  fq  =  (0,0,0).  Equation  (D3)  results  in 

/  ^1  /  ^2  f  ^3 

=  c  —  =  c  —  =  c  — 

^2  ^3 


«5=i- 


-ay  =  1- 


t  f 
^iJ 


(D4) 


(D5) 

(D6) 


because  of  the  assumption  a^  =  0.  A  linear  array  therefore  results  in  a  curve 
(2-D  solution)  in  speed-depression  angle  (velocity)  solution  space. 
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APPENDIX  E 

Calculation  of  reflection  coefficient  moments 


The  moments  of  the  reflection  coefficient  (T)  are  obtained  by  multiplying  Eq. 
(6.8)  by  the  conjugate  of  itself.  To  simplify  Eq.  (6.8),  let 


a  =  ri2. 

(El) 

II 

1 

(E2) 

C  =  r23e'^^2M_ 

(E3) 

so  that 

r  =  a+YCX(-«cr- 

n=0 

Taking  the  expected  value  of  (E4)  multiplied  by  its  conjugate  yields 

(|rp)  =  (|ap)+/w^iieaOV*^'r)+(2R={“*->^i(-“0"^  (E5) 


Noting  that 

and  summing  the  geometric  series  results  in 


(E6) 


(E7) 


Equation  (6.12)  is  obtained  by  substituting  (El)  -  (E3)  into  (E7).  Multiply  (E4)  by 
(E5)  and  take  the  expected  value  to  obtain  the  moment 


moment 
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I  OO  OO  ] 

(|rr)=|ap(|rp}+|o.pM^  ICf  S  K-aQ'hr) 


=0m=0 


+2\af{  Re- 
f 


a*YCE("“^)” 


n=0 


+1/  icr 


\( 


S  I(-“0'KC)"  s  X(-a9”(-aT)' 

V,^=0/n=0  y\^n=0/?=0 


(E12) 


-Kl 


Re 


a*YCS(““0” 


n=0 


a^CK-aCr  +«yT  iKCf  I  S(-«C/hTr  • 


n=0 


«=0 


je=0m=0 


We  use  (E7)  for  the  first  term  on  the  right  side  and  use  the  relation  (E6)  for  the 
second,  fourth,  and  sixth  terms.  The  third  term  evaluates  to  zero,  and  the  fifth 
term  uses  the  identity 


Re 


SK)"  4  IK)" 

„=0  ]j  ^n=0 


(E13) 


With  much  work,  the  fourth  order  moment  of  the  reflection  coefficient  can  be 
shown  to  be 


(rr ) = + 4|aPlYp 


■laro 


^Rel 


(a^Y*)|Yr 


(E14) 


(l-lapTi^r 


Using  the  relations  in  (E1)-(E3),  the  above  equation  is  equivalent  to  Eq.  (6.14). 
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APPENDIX  F 

Ray  approximation  calculation  of  reflection  coefficient  magnitude 


F.1  Linear  decrease  in  squared  wave  number  with  depth 

The  ray  paths  for  propagation  into  the  lower  medium  are  parabolic.  If  the 
lower  medium  were  lossless,  then  all  of  the  energy  propagating  into  the  lower 
medium  would  be  redirected  upward,  and  the  magnitude  of  the  reflection 
coefficient,  Tjs,  would  be  unity.  However,  loss  in  the  lower  medium  results  in  a 
magnitude  squared  reflection  coefficient  equal  to  the  energy  that  reaches  the 
surface.  The  square  of  the  wave  number  in  medium  3  is  given  by 

where  k2  is  the  constant  wave  number  in  the  arbitrarily  thin  layer  below  the  mean 
surface,  and  the  depth  coordinate 

C  =  -z  (F2) 

increases  with  depth.  By  defining  the  sound-speed  gradient  as  the  change  in 
sound  speed  with  respect  to  depth  at  a  point  just  below  the  mean  surface,  a  can 
be  expressed  in  terms  of  the  sound-speed  gradient  g  at  ^  =  0 ,  the  speed  ratio  v, 
and  the  sound  speed  in  medium  1,  cj. 

0  =  ^.  (F3) 

VC, 

Since  energy  propagating  through  the  sediment  will  be  attenuated  by  an  amount 
determined  by  the  line  integral  over  the  ray  path,  we  can  approximate  the 
expected  value  of  the  magnitude  squared  of  the  partial  reflection  coefficient  r23 

as 
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( 


=  exp 


-2Im  j 

;=0  sin(e(C)) 


(F4) 


where  the  path  of  integration  in  over  the  ray  path,  0(0  is  the  grazing  angle  of  the 
ray  as  it  propagates  through  the  lower  medium,  and  k2  is  the  wave  number  just 
below  the  surface  (medium  2).  In  this  example,  the  mean  depth  of  medium  2,  di, 
is  allowed  to  be  arbitrarily  small.  The  turning  point  of  the  ray,  's  found  from 
Snell’s  law  noting  that  the  grazing  angle  is  zero  at  this  point: 


COS(6)  _  €05(02)  (P5) 

C3(C)  C2 

where  02  is  the  grazing  angle  of  the  field  just  below  the  mean  surface,  related  to 
the  incident  field  grazing  angle  0i,  and  speed  ratio  V21  is  by  Snell’s  law 


cos(02)^  (F6) 

cos(0i) 

Evaluating  the  integral  in  (F4)  yields  Eq.  (6.15). 


F.2  Linear  increase  in  sound  speed  velocity  with  depth 


The  sound  speed  in  medium  3  is  set  equal  to  the  sound  speed  in  medium  2 
at  the  medium  2  -  medium  3  interface;  it  increases  linearly  with  depth,  and  is 
given  by 

C3g)=(l+<lOc2.  (F7) 


The  lower  medium  reflection  coefficient  follows  from  the  following  line  integral: 


(lr23p)  - 


exp 


-2Im^ 


C“^max  1 

h  i 


^=0  l  +  a^sin(0(Q) 


(F8) 


By  using  Snell’s  law,  the  ray  turning  point,  for  th's  constant  sound-speed 
gradient  profile  is  easily  shown  to  be 
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( 


Z  =-“ 
^max 

a 


— — 1 
005(62)  y 


(F9) 


The  parameter  a  can  be  expressed  in  terms  of  the  sound-speed  gradient  and  the 
medium  1  sound  speed. 


a  = 


g 


(F10) 


Evaluating  the  integral  in  (F8)  yields  an  expression  similar  to  Eq.  (6.15)  for  the 
magnitude  squared  reflection  coefficient: 


f 

=  exp 

V 


8ji, 
- In 


1-Hsin(e2)y 

cos(02)  Jj 


(F11) 


where 


is  the  sediment  gain  parameter. 


(F12) 


F.3  Linear  increase  in  sound  speed  over  discontinuous  interface 


in  the  next  case,  the  sediment  sound  speed  is  assumed  to  increase  linearly 
down  to  a  depth  d^,  where  the  sediment  parameters  are  discontinuous.  When 
the  maximum  depth  of  the  ray  in  medium  3,  Zn^gx-  's  less  than  the  depth  of  the 
surface  separating  medium  3  and  medium  4,  the  loss  due  to  the  path  in  the 
sediment  above  the  discontinuous  interface  is  given  by  (F8).  However,  when 
Zmax  >^3.  the  ray  intersects  and  reflects  off  of  the  medium  3  -  medium  4 
interface,  resulting  in 


Ti^=lr34pexp 


-2W 


V 


C=rf3 

h  1 


1  ll 

l  +  aC  sin(e(0)f , 


Zmax  >‘^3- 


(F13) 
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The  above  equation  is  the  same  as  (F8)  except  the  upper  limit  of  the  integral  is 
c/3,  rather  than  Z^^ax.  arid  the  magnitude  of  the  partial  reflection  coefficient  |r34| 
is  included.  Setting  a  =  gl{v2\Ci)  and  evaluating  the  above  integral,  we  get 

2  r  8jt,  fl  +  sin(e2)^^  rj 

3^^  =  exp - In  - — ^  Zmax  ^  «3 

[  q  {  008(02)  J)  _ ■ 

fl  +  sin(92)'|  {\  +  ^\-(\  +  ad2f 

-exp  ^  nj^  cos(02)  J  ”  (l  +  ac/2)cos(02) 

V  L  ^ 

^max  ^  ^3  * 
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